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ABSTRACT 


This publication contains summaries of the papers presented at the Third Spaceborne 
Imaging Radar Symposium held at the Jet Propulsion Laboratory (JPL) California Institute of 
Technology, in Pasadena, California, on January 18-21, 1993. The purpose of the s^posi^wasto 
present an overview of recent developments in the different scientific and technological fields related 
to spaceborne imaging radars and to present future international plans. 

This symposium is the third in a series of "Spaceborne Imaging Radar" symposia held at JPL. 
The first symposium was held in January 1983 and the second in 1986. 
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INTRODUCTION 


This publication contains summaries of the papers presented at the Third Spaceborne 
Imaging Radar Symposium held at the Jet Propulsion Laboratory (JPL), California Institute of 
Technology, Pasadena, California, in January 1993. The purpose of the symposium was to present 
an overview of recent developments in the different scientific and technological fields related to 
spaceborne imaging radar and to present future international plans. 

The last five years have seen major advances in the development and utilization of 
spaceborne imaging radar for Earth and planetary studies. The National Aeronautics and Space 
Administration/JPL Magellan spacecraft radar mapped the totality of the surface of Venus, which is 
hidden by the dense global cloud cover, giving us a detailed view of a very dynamic (historically) 
surface covered with volcanic, tectonic, eolian, and impact features. The European Space Agency 
(ESA) and Japan launched semioperational Earth-orbiting systems that are providing a wealth of 
data for ocean study, geologic mapping, and environmental monitoring. 

The next five years are expected to be equally exciting, opening a new frontier in the 
utilization of spaceborne radar systems. The joint U.S. /Germ any/Italy SIR-C/X-SAR will provide the 
first multispectral ("color"), multipolarization spaceborne capability using active arrays. ESA will 
launch the Earth Resource Satellite 2 (ERS-2) (then followed by Envisat), which will ensure long- 
term continuity of the ERS-1 data. This continuity is crucial for long-term monitoring of 
environmental changes. Canada will launch RADARSAT-1 (followed by RADARSAT-2), which will 
allow wide-scale repeat coverage of the world's polar regions. Japan and the U.S. will be launching 
the first rain-mapping radar, and the U.S. will launch the Titan radar on the Cassini spacecraft, 
which will orbit Saturn and map the cloud-covered world of Titan. 

The next five years also hold the potential for major new exciting missions. Radar 
interferometry is maturing as the preferred technique for global digital topography mapping, and the 
U.S. and Italy are planning a TOPSAT Mission. SIR-C/X-SAR could be transformed into a free flyer. 
Japan, France, and the U.S. each are considering advanced spaceborne systems, and the U.S. and 
United Kingdom are studying cloud-mapping radar systems. 

Collectively these missions, if coordinated, could provide a very powerful element of the 
Mission to Planet Earth to understand and monitor the global and local changes on our planet. 

We all look forward to these exciting developments and to the Fourth Spaceborne Imaging 
Radar Symposium, which will be held in Spring 1998, the twentieth anniversary of Seasat. 


Charles Elachi 

Assistant Laboratory Director 
Office of Space Science and 
Instruments 
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Magellan: Principal Venus Science Findings 

R. Stephen Saunders 
Jet Propulsion Laboratory 
California Institute of Technology 
Pasadena, California 91109 

This is a brief summary of the science findings of the Magellan 
mission, principally based on data from the radar system. Future 
plans for Magellan include acquisition of high resolution gravity 
data from a nearly circular orbit and atmospheric drag and 
occultation experiments. The Magellan science results represent the 
combined effort of more than 100 Magellan investigators and their 
students and colleagues. More extensive discussions can be found in 
the August and October, 1992 issues of the Journal of Geophysica l 
Research. Planets [1] . The Magellan mission's scientific objectives 
were (1) to provide a global characterization of landforms and 
tectonic features; (2) to distinguish and understand impact 
processes; (3) to define and explain erosion, deposition, and 
chemical processes; (4) to model the interior density distribution* 
All but the last objective, which requires new global gravity data, 
have been accomplished, or we have acquired the data that are 
required to accomplish them. 

Synthetic aperture radar imaging and altimetry were acquired over 
nearly 99% of the planet with resolution between 120 m and 300 m and 
at least four looks. Several image geometries were obtained in 
order to provide the best interpretation of the landforms. For the 
first 243 day cycle, one Venus rotation, an incidence -angle, or 
look-angle, profile was used that maximized the image resolution 
and overall quality everywhere along the orbit. This profile caused 
the incidence angle to vary from about 15" over the north pole to 
45" at the equator. In the first cycle we mapped 83% of Venus, more 
than meeting the primary mission objectives. In the second mapping 
cycle, mapping was restricted in order to control spacecraft 
temperature . Image data were obtained at a constant incidence 
angle, and looking to the right (toward the west) in the opposite 
direction from cycle 1. Also, in cycle 2, some of the major gaps 
were filled with the same incidence angle profile as used in cycle 
1. In cycle 2, we also conducted a successful test of a stereo mode 
in which we imaged at a slightly different angle than in cycle 1. 
The stereo was so useful that it was decided to devote much of the 
third mapping cycle to acquiring stereo images. All of the radar 
image data were processed at JPL in a complex flow that begins at 
the DSN stations at Goldstone, Madrid, and Canberra. In addition to 
images and altimetry, Magellan also acquired radiometer data 
whenever images were obtained. [2] • The radiometry samples the radio 
emission of the surface at the radar wavelength. Emissivity vanes 
from place to place because of variations in surface properties. 

Imaging was terminated at the end of the third cycle and gravity 
data are being acquired during cycle 4. This is accomplished by 
pointing the 3.8 m high— gain antenna toward Earth during, the 
periapsis part of the orbit and recording the returned radio signal 
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at the DSN. From this signal we extract the slight accelerations of 
the spacecraft as it orbits Venus and convert these accelerations 
into gravity maps that tell us about density variations in the 
interior. 

Magellan image data have provided several improvements in knowledge 
of the fundamental planetary constants for Venus. The rotation 
period of Venus was refined to 243 . 0185 ±0. 0001 days and the north 
pole direction, in J2000 coordinates, has been refined to right 
ascension 272.76° ±0.02° and declination 67.16° ±0.01° . The mean 
radius was refined to 6051.84 km, with the lowest point 6048.0 km 
and the highest point 6062.57 km [3]. 

Magellan has established volcanism as the dominant surface process 
on Venus [4] . Volcanism is broadly distributed, not completely 
random, but does not form linear patterns as on Earth where major 
volcanic activity tends to occur along plate boundaries. Image 
analysis reveals thus far 556 shield fields, 274 volcanoes 20-100 
km, 156 volcanoes 100 km, 86 calderas (not on shields) , 259 
arachnoids, 53 lava flow fields, 200 sinuous lava channels, 145 
steep-sided domes (pancakes) [4]. Over 360 coronae and corona-like 
features have been identified [5] . 

Tectonics is a major process [6], with evidence for extension and 
compression. Steep slopes (20“-30°), up to tens of km in extent, 
provide evidence of active tectonics. Deformation is more 
distributed than on Earth. Shear zones are seen in complex ridged 
terrain. Trench topography resembles terrestrial subduction. An 
extensive equatorial zone of fractures is among the most recent 
tectonic features. 

More than 900 impact craters 1.5 km to 280 km have been identified 
[7,8]. There appears to be a globally random distribution yielding 
an average surface age of about 500 Myr. Both bright and dark 
splotches appear to be shock signatures. Most craters are 
unmodified. Bright and dark E-W oriented parabolic halos are 
associated with about 20% of craters [9]. 

Surface processes and surface properties [10,11,12] investigations 
yield more than 8000 mapped wind streaks with directions consistent 
with Hadley circulation. Possible dune fields have been identified 
and there is widespread evidence of landslides [13]. Anomalous 
left-right reflectivity behavior indicates unusual surface 
reflectivity behavior, possibly caused by asymmetric shapes. 
Anomalous low emissivity in elevated regions has been confirmed 
[ 2 ]. 

Major questions about Venus remain unresolved, pending acquisition 
of new data and further analysis of existing data. Interpretation 
of the impact crater population suggests a major secular change in 
the rate or style of resurfacing [7], but the details or even the 
reality of this change, whether catastrophic, cyclical, local or 
global scale is unknown. High— resolution global gravity data will 
help address some of the unresolved issues concerning the 
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ERS-1 MISSION OVERVIEW AND STATUS 
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INCLINATION CORRECTION 

INCLINATION CORRECTION PERFORMED 16 DECEMBER 1992, . 

satellite brought back from +1 Km to -1 Km boarder at argument of latitude 90 
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ERS-1 

ESOC Mission Management and Planning Office - Pleval 
Reference Measurement System 
1 6: Long term stability rest, bit clock link 2 PB (USO) 
Nominal Frequency: 15000000.000 

In ffreq— nom) : Mean = 0.181 Hz, Sigma = 0.051 Hz 

In (freq-nomj/nom: Mean = 1 .206e-08, Sigma = 2.068e-09 


4.0x10 


-9 


.’8 


0.0 

30.0 

60.0 
90.0 =- 

120.0 ^ 

150.0 |- 

180.0 |- 

210.0 e- 

240.0 §- 

I 270*0 §- 

3 - 

O = 

| 300.0 E- 

'o E 

§* 330.0 §- 

o E 

360*0 |* 

390.0 E- 

420.0 E- 

450.0 §- 

480.0 |- 

510.0 |- 

540.0 E- 

570.0 E- 


H I I i I I I M j i i i II l I I l | M I I II I IT 


600.0 


Relotive: (compensoted_frequency - nominal) / nominal 
8-OxlO -9 1.2x10“ 8 1.6x10 -8 2.0x10~ 8 2.4x10~ 8 2.8x10 w 3. 

1 1 1 rnrrn 1 1 1 m 1 1 1 1 i j n n i n n | t i i 1 1 i rn immmi 


M t I 11 I I 1 1 II 1 i 1 t 1 1 i 1 1 1 1 1J-.11. 


2x10 

TT3 


0 143 Hz 06- AUG -91(2 18) 09:40:37, 291. L 
.179 Hz, 21 -AUC-9 1(233) 09:40*6, 506, L 



0.141 Hz, 02— JUH— 92(154) 11:09*6,4605, L 


0.120 Hz, 1 4— JUL-92( 1 96) 10*0:10*5206, L 


7.00000*- 09. 0.105 Hz, 25- AUG -92(238) 10*0^1*5807, L 


6.26667*— 09, 0.094 Hz* 29- SEP -92 (273) 10*0*0,6308, L 


5.40000*— 09, 0.081 Hz, 03-NCV-92(308) 10*0:17,6809, L 


4.46667*— 09, 0*67 Hz, 08-0EC-92(343) 10*0:07,7310, L -3 


lllllllll lllllMmllllll LLi-L 


[ 1 I : t 1 i I I I I 1 I I 1 t Eh 


0.00 


0.06 


0.12 


0.36 


0.18 0J24 0.30 

(comp«nsated_frequency — nominal) Hz. 

Legenda: relative, freq— nom, date, year^day, hh:mm:ss, orbit, flag 

Flog V: Valid measurement 

Flag L: Volid measurement: Limits exceeded 

32 


0.42 


0.48 


ORIGINAL PAGE IS 
OF POOR QUALITY 




Internal Gain 



8P 


34 


28-JUL- 1 991 28-NOV- 1 991 29-MAR- 1992 30-JUL-1992 30-N0V-1 t m. 




System Noise 



35 


25-JUL- 1 99 1 25-NOV-1991 28-MAR-1992 29-JUL-1992 30-NOV-1992 





INSTRUMENTATION (5) 



ERS-1 MISSION OVERVIEW AND STATUS 


INSTRUMENT OPERATIONS 


o 



37 


ERS-1 MISSION OVERVIEW AND STATUS 


PLATFORM AND INSTRUMENT AVAILABILITY 
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ERS-1 MISSION OVERVIEW AND STATUS 


ERS-1 OVERALL GROUND SEGMENT 
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ERS-1 MISSION OVERVIEW AND STATUS 






LBR DATA ACQUISITION 
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ERS-1 MISSION OVERVIEW AND STATUS 


LBR FD PRODUCTS 
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-1 MISSION OVERVIEW AND STATUS 



LBR FD DISSEMINATION 



ERS-1 MISSION OVERVIEW AND STATUS 


SAR FAST DELIVERY PROCESSING AND DISSEMINATION 
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ERS-1 MISSION OVERVIEW AND STATUS 



NEAR REAL-TIME DISSEMINATION 
OF ERS-1 SAR IMAGES OVER EUROPE 
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ACQUIRED SAR FRAMES FROjM BEGINNING OF MISSION 



ERS-1 MISSION OVERVIEW AND STATUS 



SAR FRAMES PLANNED FOR ACQUISITION 

(STATUS AS OF 30 OCT. 92, PLANNING UP TO 28 FEB.93) 
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AN EXPANDING USER COMMUNITY 
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ERS-1 MISSION OVERVIEW AND STATUS 


IMPROVEMENTS TO SYSTEM AND SERVICES 
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ERS-1 MISSION OVERVIEW AND STATUS 


OBJECTIVES STILL TO BE ACHIEVED 
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ALMAZ 
V. Viter 

NPO Machinostroyenia 
33, Gagarin Street 
Reutov, Moscow District 143952 


N94- 15 889 


/ 

ir? 


Russia 


R AST r n AT A OF THE AUTOMATIC SPACE STATION. 

"ALMAZ- IB" 

Launch vehicle 

Proton 

Operational orbit parameters: 


• Altitude 

• Inclination 

350-400 km 
73 grad 

On-orbit mass 

18,550 kg 

On-board scientific equipment mass 

up to 4,500 kg 

Active operation time 

up to 3 years 

Orientation & stabilization accuracy: 


• Orientation 

• Stabilization 

15-30 min. of arc 
4-6 min. of arc 

Power-generating system characteristics: 


• Available average power per turn 

• Maximum power (at peak load up 

to 15-20 min.) 

2,300-3,300 W 
8,600 W 


AT MA7-1B AUTOMATIC S PACE STATION 


1 . SAR - 70 
SAR - 10 r 

2. SAR - 10 1 

3. Multizonal Scanner MSU-SK 

4. Spectroradiometer 

5. Multizonal Scanner MSU-E 

6. Side-looking Radar 

7. Optoelectronic Equipment 
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PRTNrTPAT . TECHNICAL CHARACTERISTICS OF REMOTE SENSING 
F.OT TTPMENT ON THE ALMA Z- IB SPACE STATION 

A Set of Variant 1 


1 . 

1 . 1 . 


Principal technical characteristics of radar equipment 


Synthetic aperture radar 


SAR-10 


Operating wavelength range 
Polarization 


Scanning direction 
Angle of sight from nadir 
while turning 

Scanning field for re-shooting 
while turning 
Time of continuous survey 
Mode of operation detailed 

scanning field 25 - 50 km 

resolution 5 - 7 m 


9, 6 cm 

- horizontal (port) 

horizontal or vertical 
(starboard) 

port or starboard at will 

25 - 50 grad 

320 km 
up to 1,200 s 

general view 
100 - 150 km 
15 m 


1.2. Synthetic aperture radar 


SAR-70 


Operating wavelength range 
Polarization 
Scanning direction 
Angle of sight from nadir 
Scanning field for re-shooting 
while turning 

Time of continuous survey scanning 

Scanning field 

Resolution 


70 cm 

horizontal or vertical 

starboard 

25 - 50 grad 

320 km 
up to 1,200 s 
100 - 150 km 
15 - 25 m 


1.3. 


2 . 


Side-looking radar 


Angle of sight from nadir 
while turning 
Operating wavelength 
Swath width 
Resolution 

Optoelectronic equipment 


25 - 55 grad 
3, 6 cm 
450 km 
400 - 1,200 m 

for stereophotography 


Number of operating ranges 5 

Resolution on the ground 2 - 4 n 

Swath width 70 km 
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3. Characteristics of multizone high-resolution 
electronic scanner MSU-E 

Number of operating ranges 
Resolution on the ground 
Swath width 
Potential angles of view 
Number of on-board sets 


4. Characteristics of multizone middle-resolution 

optomechanical scanner MSU-SK 

Number of operating ranges 
Resolution on the ground 
in visual range 
- in IR range 
Swath width 

5. Characteristics of” spectroradiometer for 

ocean satellite monitoring 

Scanning field 
Number of spectral zones 
Boundaries of spectral zones 
Linear resolution at nadir 
Temperature resolution 

6. Information transmission 

Direct transmission to available receiving stations located on the 
Russian territory and abroad 

Data collection on board and transmission to these stations when 
flown over 

Data collection on board and transmission via relay-satellite to 
central receiving station 

7. Effectiveness of information reception 

Directly to receiving stations no more 12 h 

via relay-satellite no more 2.5 h 

8. Rate of information transmission 

Directly to receiving stations. 

Via relay-satellite 
Direct transmission to mobile and 
small receiving stations 


122.8 Mbit/s 
10 Mbit/s 
3 Mbit/s or 
960 Kbit/s 


2 x 1,100 km 
11 

0.4 - 12.5 mkm 
0.6 km 
0.1 K 


5 

80 

300 m 
300 km 


3 

10 m 

2 x 24 km 
32 grad 
2 
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MATN FIELDS AND DIRECTIONS O F ALMAZ- IB 
SPACE INFORMATION APPLICATION 

CARTOGRAPHY 

Updating and completing topographic and thematic maps 
Cities development planning 
Transport networks planning 
Mapping of shelf shallow sections 

LAND MONITORING 

Information supply of earth catastrophe, lands, forests, pastures and 
water resources 

Compiling soil maps and determining humus content 
Taking inventory and forecasting harvest yield 
. Determining soil temperature and moisture content 
> Detecting underground water sources 

GEOLOGY 

Revealing perspective areas for mineral prospecting 
Geologic mapping 

Planning and controlling mineral/raw-mineral base 
Specifying data of mining and oil-bearing regions 
Monitoring large open-cast mining development 

ECOLOGICAL MONITORING 

Controlling chemical compounds content in soil 
Controlling technogenic actions 
Controlling soil salinization and desertness 
Controlling soil erosion 

Controlling pollution of water surfaces with petroleum products 
Revealing underground unpermitted burials of industrial waste 

OCEANOLOGY, PILOTAGE, FISHERY 

• Evaluating ice conditions for pilotage 

• Controlling chlorophyll content, evaluating bio-productivity of the 
world's ocean regions 

• Studying the world's oceans and seas 

• Compiling ocean temperature maps, pollution recording, studying 
fronts, currents, water, wind, etc. 

INFORMATION SUPPLY DURING EMERGENCY 

• Controlling situation in geologically dangerous regions 
(earthquakes, mudslides, avalanches, volcanic eruptions) 

• Detecting forest and peat fires 

• Controlling situation during accidents at industrial sites (fires, 
explosions, construction crashing down) 


57 


• Controlling situation in large transport accident regions including 
underground product tubes 

• Controlling situation in hydrodynamic accident regions (dams, locks, 
sewage works) 
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SPACEBORNE IMAGING RADAR-C INSTRUMENT 

BRYAN L. HUNEYCUTT 
Jet Propulsion Laboratory 
California Institute of Technology 
4800 Oak Grove Drive, Pasadena, California 91109 


Abstract - The Spaceborne Imaging Radar-C is the next radar 
in the series of spaceborne radar experiments, which began 
with Seasat and continued with SIR-A and SIR-S. [i-4] The 
SIR-C instrument has been designed to obtain simultaneous 
multi frequency and simultaneous multi polarization radar 
images from a low earth orbit. £51 It is a multi parameter 
imaging radar that will be flown during at least two dif* 
ferent seasons. The instrument operates in the squint 
aligraent mode, the extended aperture mode, the scansar 
mode, and ‘the interferometry mode. The instrument uses 
engineering techniques such as beam nulling for echo track- 
ing, pulse repetition frequency hopping for Doppler centroid 
tracking, generating the frequency step chirp for radar 
parameter flexibility, block floating-point quantizing for 
data rate conpression, and elevation beamwidth broadening 
for increasing the swath i l lunination. 


I. INTRODUCTION 


The Spaceborne Imaging Radar-C (SIR-C) is a dual- 
frequency, quad-polarization imaging radar, designed to fly on the 
Space Shuttle in a low Earth orbit, and to operate from a stable 
platform located in the Orbiter payload bay. [6] Several flights are 
currently planned for SIR-C beginning in 1994. SIR-C has been 
designed to operate simultaneously at both L-band and C-band 
frequencies, and to utilize quad-polarization returns at each 
frequency. [6] SIR-C effectively has four separate radars: L-band 
horizontal Lh, L-band vertical L v » C-band horizontal C H , and C-band 
vertical C v . In the nominal mode, SIR-C operates each of its four 
radars by radiating at broadside, utilizing a fully focused antenna 
aperture, much like the previous experiments, except that now there 
are really four radars operating simultaneously. The important SER-C 
instrument parameters related to the operating modes and techniques 
are given in Table 1. 


Fully instrumented testing of these various modes and 
techniques will occur at the experiment developer’s facility, and not 
at the Shuttle launch site. Transmission from the full aperture and 
reception at far range is not planned even at the experiment 
developer’s facility, due to the complicated and costly logistics of 
erecting such a large structure in an environmentally protected range. 
The SIR-C antenna is presently planned to be delivered into the 
integration flow at the launch site in its ready-to-launch configuration. 


TABLE 1 

SIR-C Instrument Ciiaracttristtcs 


Parameter 

L-Band 

C*Band 

Azimuth Coalignment (arc min) 

±8 

±2 

Azimuth beam steering (degrees) 

±7 

; 1 

Extended aperture beams 

127 

127 

Extended aperture dwell time ( ms) 

20 

20 

Scansar beams 

4 

4 

Scansar dwell time ( ms ) 

30 

30 

Interferometry phase center offset (m) 

N/A 

4 

Beam null ratio 

l /PRF 

t /PRF 

PRF hopping dwell lime (s) 

1 

1 

Number of selectable PRF a 

16 

16 

Elevation beamwidtha 

5-18* 

5-18* 

Number of selectable beamwidths ' 

8 

8 

Digital chirp pulsewidth (jts) 

8.4, 16.9, 33.8 

8.4, 16.9. 33.8 

Digital chirp bandwidth ( MHz) 

10,20 

10,20 

Quantization (bit) 

8 

8 

BFPQ subset (bit) 

4 

4 

BFPQ block length (samples) 

128 

128 

PRF (Hz) 

1240-2160 

1240-2160 

Stalo frequency (MHz) 

89.994240 

89.994240 

Data steering 

offset video 

offset video 

Number of high rate record channels 

4 

4 


Functional tests of the SIR-C instrument will occur during integration 
with the Shuttle. 

II. ANTENNA OPERATION AND TECHNIQUES 

SIR-C will operate in modes, which were previously not 
possible, without electronic beam steering. Figure 1 illustrates the 
target illumination geometry corresponding to these operating modes. 

A. Squint Alignment 

Electronic beam steering in azimuth may be needed to align 
the H and V polarization swath illuminations, even though the 
antenna design calls for an azimuth electrical-to-mechamcal bo resight 
error of less than 4 arc minutes. Azimuth steering allows SIR-C to 
take data at a given squint angle, whereby the azimuth angle is fixed 
to one side or the other with respect to broadside. The H- 
polarization and V-polarization illuminations can be unintentionally 
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offset due to mechanical tolerances in manufacturing and 
environmental effects in space such as uneven thermal distribution 
and zero gravity unloading. The beams can then be commanded to 
different azimuth squint angles until they are coaligned to within 4 
arc minutes. The shape of the azimuth spectrum of the cross- 
polarization return can be used to verify proper co-alignment. 
Azimuth steering also allows experimenters to collect synthetic- 
apefture radar (SAR) data at selected fixed squint angles. 

Because of the limited area for locating the electronics on 
each panel, a single transmit/receive (T/R) module and phase shifter 
pair feeds a subarray, and not individual elements. The stick-like 
subarray contains 6 radiating elements for L-band and 18 radiating 
elements for C-band. This restricts the azimuth steering to only a 
few degrees. Electronic beam steering can be accomplished for only 
±1 deg in azimuth for C-band, and .+2 deg in azimuth for L-band; 
beyond these limits, the side lobes increase, the main lobe broadens, 
and grating lobes occur at unacceptable levels. [6] 

B. Extended Aperture 

Azimuth steering allows SIR-C to take data in the extended 
aperture mode, or "spotlight" mode, whereby the same scene on the 


ground is illuminated as the Shuttle passes the scene, and the 
effective aperture is thereby extended. This technique increases the 
effective azimuth Doppler bandwidth to improve the image quality. 
This increased bandwidth permits: (1) finer spatial resolution by 
utilizing the full bandwidth and (2) reduced speckle by incoherently 
adding more azimuth looks. To implement this extended aperture 
mode, as an example, the antenna beam is first squinted ahead a: the 
full range positive squint angle, and then decremented in small 
uniform steps of squint angle at uniform time increments, until the 
antenna beam is looking behind at the full range negative angle. 

C. Scansar 

Electronic beam steering in elevation allows SIR-C to 
image in a scansar mode, whereby the antenna beam is steered to as 
many as four previously selected elevation angles during each 
synthetic aperture interval. This dramatically increases the width of 
the swath which can be imaged during a data take. The beam is 
electronically steered in elevation as in azim uth by varying the 
phase-front across the aperture via the microprocessor-controlled 
phase shifters. The SAR processor must treat the data from each 
elevation angle as bum data. The azimuth resolution of the 
processed bursts of data is degraded by a factor of 4; however, the 
swath width illumination is increased by about the same factor. 
Since the synthetic aperture interval increases with the look angle off 
nadir, the dwell time at each of the four elevation angles increases 
with the look angle. 

D. Antenna Beam width Broadening 

For a uniformly illuminated aperture, the aperture 
dimension and RF wavelength determine the antenna beam width. 
SIR-C uses a technique of tapering the power in elevation, which 
provides lower side lobes (approximately -18 dB side lobes) than for 
the uniform illumination case. This configuration uses less T/R 
modules, and thus less power. Lower side lobes are achieved in 
■ elevation. The resultant antenna pattern is lower in amplitude and 
has a wider beamwidth. 

Within the constraints of the Shuttle payload bay volume, 
the aperture sizes were optimized to have sufficient aperture gain, yet 
sufficiently low azimuth and range ambiguities over the pulse 
repetition frequencies (PRFs) selected. The sizes of the full apertures 
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arc 12.1 m x 2.95 m for L-band and 12.1 m x 0.75 m for C-band, 
which yield a fully focused elevation beamwidth of 4.7 deg for each 
frequency. Variable beamwidths in elevation can be accommodated 
by selecting a pre-programmed phase function across the array. One 
of eight selectable beamwidths from 4.7 deg to 18.0 deg is selected 
by an uplink command. The amplitude decreases and the side-lobe 
levels increase slightly with increasing beamwidths. The beamwidth 
is optimized for swath illumination at the various look angles. 

E. Beam Nulling 

The phase shifters in one half the array will be periodically 
shifted 180 deg in elevation, thereby nulling the elevation beam 
pattern. This is an echo-tracking technique to determine the roll 
angle, and thereby accounts for the elevation antenna pattern in the 
raw SAR data during amplitude calibration. Results from simulations 
indicate that beam nulling during one receive interval in each one- 
second interval is sufficient to determine the roll angle and yet not 
degrade the image significantly. The roll angle can be determined 
with an accuracy of a few tenths of a degree. 

III. RF ELECTRONICS 

The RF electronics subsystem is located on the pallet in the 
cargo bay. The RF electronics subsystem accommodates several new 
techniques not previously used in the SIR series. 

A. Digital Chirp 

To attain fine resolution in range, SIR-C encodes each 
transmitted pulse, such that each pulse of duration T can resolve 
targets as if a much shorter pulse of duration T/TBW, where TBW 
is the time-bandwidth product, had been transmitted. SEAS AT , SIR- 
A, and SIR-B used the passive dispersive delay line (DDL) approach; 
however, for SIR-C a different DDL would be required for every 
TBW required. SIR-C attains this large time-bandwidth product by 
distributing the energy of each pulse over the frequency bandwidth, 
stepwise linearly with time. The SIR-C chirp signal is approximated 
digitally, in that a tone is successively stepped across the bandwidth 
within the pulse duration, approximating a linear FM (LFM), or 
chirp, signal. The number of frequency steps varies with the 
bandwidth, pulse-width, and frequency-step size combination. The 
10-MHz, full pulse- width case requires 76 steps; whereas, the 20- 


MHz, full pulse-width case requires 152 steps. The number of steps 
for SIR-C trades off the integrated side-lobe ratio with realistic 
switching rates. The integrated side-lobe ratio has been measured to 
be approximately -9 .5 dB and the measured impulse response width 
agrees closely with that of an ideal chirp. 

The frequency-step chirp technique provides flexibility in 
the selection of certain radar parameters: (1) pulse width, (2) 

calibration tone frequency, and (3) bandwidth. The transmitted pulse 
width can be decreased by a factor of 2 or 4 to lower the average dc 
power usage. Fixed frequency calibration tones can be generated 
from the same digital chirp device, and the frequency of the 
calibration tone can be selected. The RF bandwidth can be increased 
from 10 MHz to 20 MHz to improve the range resolution. 

B. Interferometry Mode 

The L-band and C-band apertures are subdivided into 3 
equal area sections called leaves. The two outer leaves are unfolded 
after attaining orbit in order to transmit. The received echo signals 
from the two extreme C-Band leaves can be coupled into two 
separate receiver channels in order to accommodate a special 
interferometry mode. This requires couplers in the C-Band RF feed 
system at the i nitial 3-way power combiner, such that the inputs from 
the two extreme leaves can be routed to two separate receiver 
channels. The technique provides returns from the two existing 
extreme leaf apertures, whose phase centers are separated by the 
distance of one leaf in the azimuthal direction. Processing of the 
interferometry data from ocean returns, for example, can provide 
ocean current direction. [7] 

IV. DIGITAL ELECTRONICS TECHNIQUES 

The Digital Electronics Subsystem uses engineering 
techniques required by SIR-C in order to obtain calibrated data and 
to be much more flexible in radar parameter selection than in 
previous imaging radar experiments in this series. These techniques 
are utilized during the normal operation of the radar. 

A. Alternating Pulse Mode 

The SIR-C exciter operates in an alternating pulse scheme 
such that the oppositely polarized pulse is delayed by half an inter- 
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pulse period. The radar effectively operates at twice the nominal 
PRR For quad- polarization operation, both like- and cross- 
polarization echoes are received in the same channel. 

B. Block Floating-Point Quantizing 

The offset video output from each of the four receivers is 
digitized to 8 bits per sample with uniform quantization at a constant 
rate of 45 MHz. A SIR-C feature allows formatting the raw samples 
as 4 bits, 8 bits, or an (8,4) floating-point block. The block floating- 
point quantizer (BFPQ) derives its name from "blocks" of data being 
uniformly quantized and subsets of the available bits being selected 
by a predetermined algorithm, equivalent to moving the "floating 
point" marker in binary data. If the data transmission rate to the 
onboard recorders is held constant, then doubling the number of bits 
per sample approximately halves the swath width, therefore trading 
off dynamic range for swath width. 

The (8 ,4) BFPQ is a data rate compression technique. The 
(8,4) BFPQ method provides an output rate similar to a 4-bit uniform 
scheme but with the dynamic range of an 8-bit system. For each 
block of data, a series of bits are transmitted, including the sign bit, 
followed by the optimally selected 3 contiguous bits per sample, 
followed by a common exponent for the block. The block size is 128 
samples of raw data. The formatted data are transmitted as 8-bit 
parallel words. The distortion noise for the (8,4) BFPQ is about 24 
dB higher than for the 8-bit uniform quantizer (UQ); however, the 
range of return signal level for which the distortion noise is almost 
constant is about 24 dB wider than that for the 8-bit UQ. [6] 
Therefore, the performance of the (8,4) BFPQ is less sensitive to the 
typically widely varying level of the return echo. 

C. Digital Data Steering 

In normal data takes , each receiver output is directed to a 
previously selected, separate digital data handling assembly (DDHA) 
for digitizing, buffering, and formatting. This steering occurs at 
offset video. Once the individual offset video signals have been 
digitized, formatted, and serialized, there is flexibility in steering to 
a redundant on-board Payload High Rate Recorder (PHRR). This 
data steering allows for recovery in case one or more of the DDHAs 
or the PHRRs should fail. 


In extended swath data takes, one receiver output is 
directed to two DDHA’s, such that the resulting data represent an 
image swath twice that of a single DDHA. The SAR processor will 
have all four channels combined on one tape cassette. 

D. Integer PRFs 

SIR-C has incorporated integer PRFs, whereby there are 
an exact integer number of PRF pulses in one second. SEASAT, 
SIR-A, and SIR-B used a countdown scheme to derive PRFs from the 
STALO, a stable local oscillator, without regard for whether the 
PRFs were integral. The SIR-C integer PRF facilitates timing for the 
radar since the PRF and timing circuitry provides basic control of the 
radar for PRF changes, receiver gain, exciter timing, digital window 
position, and data channel switching, all of which change state 
synchronouslyon one-second time ticks (OSTTs). Also, SER-C chose 
the PRFs such that there are an integer number of 8-bit bytes in a 
range line. This feature facilitates ground processing. Since the 
PRFs are derived from a very stable STALO, the requirement of 
integer PRFs affects the choice of the STALO frequency. By 
analysis, it was determined that by using a STALO with frequency 
89.994240 MHz, there were at least 16 integer PRFs in the range 
from 1240 Hz to 2160 Hz, which were compatible with the altitude 
and look-angle range of SIR-C. 

E. PRF Hopping 

SIR-C uses a PRF hopping technique, or multiple PRF 
technique, whereby the azimuth Doppler spectrum can be located 
unambiguously. For the PRF hopping technique, if the azimuth 
spectrum is offset in frequency because of pointing errors, for 
example, the centroid of the spectrum falls into different portions of 
the azimuth processing bandwidth when sampled at various PRF 
rates. By knowing the ambiguous location of the centroid for three 
different PRFs, the location of the unambiguous location, and thereby 
the yaw angle, can be calculated. The PRF hopping will occur at the 
beginning and end of each data take. There will be a one-second 
dwell at each of the three PRFs. 

V. Summary 

Several operating modes and engineering techniques used 
in the SIR-C instrument have been described: squint alignment, 
extended aperture, scansar, interferometry, beam nulling, PRF 
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hopping, antenna beamwidth broadening, generating frequency-step 
chirp signals, block floating-point quantizing, integer PRF’s, and 
digital data steering. 
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Int r oduc t ion 

X-SAR is the German/ Italian contribution to the NASA/JPL Shuttle 
Radar Lab missions as part of the preparation for the Earth 
Observation System (EOS) program. The Shuttle Radar Lab is a 
combination of several radars: an L-band (1.2 GHz) and a C-band (5.3 
GHz) multipolarisation SAR known as SIR-C (Shuttle Imaging Radar) 
and an X-band (9.6 GHz) vertical polarised SAR which will be 
operated synchronously over the same target areas to deliver for 
the first time calibrated multifrequency and multipolarisation SAR 
data at multiple incidence angles from space. 

A joint German/ Italian project office at DARA (German Space Agency) 
is responsible for the management of the X-SAR project. The space 
hardware has been developed and manufactured under . industrial 
contract by Dornier and Alenia Spazio. Besides supporting all the 
technical and scientific tasks, DLR, in cooperation with ASI 
(Agenzia Spaziale Italiano) is responsible for mission operation, 
calibration and high precission SAR processing. In addition, DLR 
developed an airborne X-band SAR to support the experimenters with 
campaigns to prepare for the missions. 


The main advantage of adding a shorter wavelength ( 3 cm ) radar to 
the SIR-C- radars is the X-band radar's weaker penetration into 
vegetation and soil and its high sensitivity to surface roughness 
and associated phenomena. The performance of each of the three 
radars is comparable with respect to radiometric and geometric 
resolution. 
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Instrument development 


As of about one year before the first mission, the flight hardware 
has been manufactured, and the integrated X-SAR instrument is 
undergoing a phase of thorough testing and characterization of its 
performance. Following to these tests, the X-SAR has been 
delivered to the U.S.A. and will be integrated in the SIR-C radars 
and subsystems at JPL (Jet Propulsion Laboratories in Pasadena) 
After the three X-SAR slotted waveguide array antenna panels are 
mounted on to the common mechanical antenna structure the antenna 
become for the first time a full 12 m in length. The X-band 
antenna will be mechanically pointed to its targets while the SIR-C 
antennas have the capability to point electronically to the same 
target area on ground. The passive X-SAR antenna is connected to 
its high power transmitter and its low noise receiver via a quite 
long waveguide running from the electronic boxes down in the pallet 
bay up to the top of the antenna structure. The three panels are 
fed via a rotary joint and a three-way power devider. The predicted 
mechanical loads and stresses, up to 15 g during launch and 
landing, are a challenge to the design of this waveguide, which has 
to be built out of special rigid and flexible pieces to take care 
of the dynamic and static deflections between the pallet-mounted 
equipment and the antenna as well as between the movable X-SAR 
structure and the core structure. 


X-SAR/SIR-C Blockdiagram 



PHRR's 


68 






Although it is planned to operate X-SAR and SIR C synchronously, 
each radar is also capable of autonomous operation. . To insure 
proper synchronism SIR-C acts as the master for all timing signals 
as well as for the local oscillator, which is the basis for the 
common frequency generation. Also part of the telemetry and 
command subsystems and the high rate data handling, recording and 
switching is under SIR-C control. Therefore, in addition to the 
testing of the subsystems itself, careful testing of this 
interfaces and the compatibility of both systems running, has to be 
verified before delivery to the Kennedy SpaceFlight Center for 
integration with the orbiter. The antenna will cover about 60 % of 
the whole cargo bay and is mounted with a fixed tilt of 14 degrees. 
jQl electronic boxes are installed on coldplates under the antenna 
structure on a spacelab pallet. Inside the shuttle cabin one of the 
three high rate data cassette recorders will be dedicated to record 

X-SAR data. 

Mission operations 

The three radars will be operated during the ten day missions 
synchronously over selected target areas based on the mission plan, 
worked out by the SIR-C /X-SAR science team in collaboration with 
the principal investigators .Highest priority is given to the 
so-called super test sites in which ground calibration and insitu 
measurements are made. In total, about 50 hours of X-SAR raw data 
are expected which will be recorded on board. A part of this data 
can in addition be transmitted to the ground in .real. time. The 
orbit altitude is planned to be 220 3cm and the inclination will be 
57* . The orbiter will fly in a 26’ roll and a 180 pitch angle 
which makes the radars looking 40* off nadar when no beam pointing 
is commanded. 
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To carry out the ambitious mission plans, a powerful and flexible 
mission operations and planning tool is being developed which is 
capable of performing joint, coordinated operations with SIR-C. 

The main elements of the MPOS (mission planning and operations 
system) are the mission planning and analysis system, the telemetry 
system, the command system and the SAR real time quicklook 
processing system. While the mission planning software tools such 
as the mission timeline planning, the data- take selection, the 
performance estimator and the command plan generator are also 
necessary for pre-mission preparation, these and the hardware 
elements mentioned above will all be integrated into the JSC POCC 
(payload operation and control center at the Johnson Space Flight 
Center) in Houston. 

The flow of the mission planning is normally the following: Input 
data provide latest information about the sites, the attitude 
timeline and the orbit state vectors. The mission planning software 
then calculates and displays the data take opportunities, which are 
then checked against the requirements and mission rules by the 
data-take selection software. To determine the optimum and jointly 
coordinated radar parameter settings, a performance estimator 
calculates radar parameters and the predicted image performance 
using known instrument characteristics for each target or 
data-take. These radar settings are then transformed by the command 
pl an generator into commands which are tagged with the appropriate 
times for execution in the instrument. Up to 256 command blocks may 
be sent to the onboard instrument sequencer memory, enough for a 
whole day of automatic operation if no updates would be necessary. 

Real time commanding is possible and might be necessary during the 
mission when the performance analyser system shows the operator 
improper radar parameters, or if the parameters are no longer 
optimum. This might happen because the Space Shuttle is not a very 
stable platform for a SAR and the attitude and orbit data avaiable 
during mission are not very accurate. The monitoring of the 
on-board instrument is performed with the telemetry system which 
shows the status and health of the X-SAR instrument as well as with 
the SAR processing system when realtime data transmission is 
scheduled for X-SAR. The real time data will be recorded on ground 
and in parallel routed to the deformatter and quicklook SAR 
processor, which show the processed image of the scene in a 
waterfall display. This will permit the monitoring of the 
performance of the instrument, the proper timing of the receive 
window and the pointing of the antenna to support the control of 
the test-site coverage. 
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After the mission an as-run timeline and a mission operations 
report evaluating the instrument functional performance will be 
supplied, together with a duplicate of all recorded data for the 
off-line processing systems. 

Performance critical areas 

A key instrument parameter for the performance is the in flight 
antenna pattern, but its prediction has a high uncertainty . Each of 
the three 4 m long antenna panels had to be measured seperately on 
a far-field test range. An ideal reconstruction of the full 12 m 
antenna performed using software tools showed excellent 
performance. 

The inflight antenna pattern will be more or less degraded by the 
mechanical performance of the antenna truss with its static 
deflections and temperature effects but also due to the presence of 
the SIR-C antenna and the orbiter elements. A further problem area 
is the knowledge of the pointing of the antenna to the target which 
is depending on the attitude measurement accuracy and the 
alignement accuracy of the antennas with respect to the orbiter 
coordinate system. Main regions of concern are the 1. sidelobe 
levels (up to -8.3 dB in a worst case) of the azimuth pattern which 
impacts on the ambiguity ratio and the peak gain variation (up to 
0.9 dB in worst case) which impacts on the gain stability and 
calibration accuracy. A verification of the inflight antenna 
pattern and its pointing will be measured with ground receivers 
distributed over the swath at least over the calibration test site. 

Gain variation of the transmitter and receiver chain over the 
temperature range is an other critical area which directly impacts 
on the radiometric accuracy or at least complicates the 
calibration. X-SAR was not successfull in implementing an automatic 
compensation of the gain variation over the temperature range of 
about 40 degrees Celsius. But for the postmission data processing 
a careful characterisation of the gain versus temperature under 
thermal vaccum conditions has been done resulting in a set of 
curves which have to be used together witch the temperature 
housekeeping recorded during the mission. The measured gain 
variation is about 0.2 dB per degree Celsius which should not be 
too bad taking into account the excellent temperature stability of 
the coldplates in earlier missions. 
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Performance prediction 


Based on jointly agreed upon definitions, X-SAR and SIR-C calculate 
their predicted performance with validated system perfomance 
software tools. The performance models use algorithms which 
transform the instrument parameters and assumptions about the 
platform, the- processor and the target into image performance 
parameters. The instrument parameters used were taken from the 
latest measurements of the flight hardware. The performance is 
calculated for point targets, providing geometric resolutions. For 
distributed targets radiometric resolution, ambiguity levels, 
signal-to-noise ratios and the performance swath width are 
calculated. The following table shows the performance for four 
typical incidence angles. 


Requirement 

for 

point 

targets 

distributed 

targets 

Incidence angle 

20' 

30* 

40 * 

50* 

20* 

30’ 

40’ 

50’ 

Range resolution 









fine < 10 m 

9.6 

9.8 







coarse < 20 m 



17.7 

18.0 





Range ISLR >11 dB 

12.1 

12.1 

12.2 

12.2 





Range PSLR >20 dB 

21.9 

21.9 

21.9 

21.9 





Azimuth resol. 30 

27.5 

27.5 

27.4 

27.4 





Azimuth ISLR >lldB 

15.47 

15.45 

15.51 

15.51 





Azimuth PSLR >20dB 

21.9 

21.9 

21.9 

21.9 





Ambiguity Azimuth 





19.4 

19.4 

20 . 0 

19.4 

Ambiguity Range 





>40 

>40 

>40 

>40 

Total > 18 dB 





19.4 

19.4 

20.0 

20.0 

Radiom. Res. <2. 7 5 





<2.1 

<2.8 

<2.0 

<2.1 

Image SNR 





>8.7 

>8.4 

>10.2 

>7.7 

Dynamic > 20 dB 





32.9 

34.7 

36.7 

36.8 

Swath width>15 km 





30.4 

20.7 

36 . 6 

32.0 


SAR Data Processing 

The processing of the X-SAR raw radar data will be performed in two 
processing and archiving centers one at ASI in Matera Italy (I- 
PAF) and one at DLR in Oberpfaf fenhofen Germany (D-PAF) . ASI will 
be responsible for processing the survey products of all X-SAR data 
within three months after the mission for the selection of the 
test sites for which high precision calibrated products will be 
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processed. The processing of the high precision products will be 
shared between ASI and DLR. The calibration tests ites will be 
processed first, so that the calibration team can determin the 
calibration factors to be used from both PAF's for the precision 
processing. Several data products like single look, multilook, 
geocoded and several product media can be selected by the data 
users which place their orders via the science team to the PAF's. 
All three processing centers JPL, ASI and DLR have jointly agreed 
on the data products, quality and the formats. The main challenge 
for the X-SAR data processing is the relative unstable Shuttle 
plattform and the low precision of the orbit and attitude data. To 
find the correct doppler and dopplerband a new algorithm to solve 
the doppler ambiguity problem in the X-SAR data has been developed 
and will be applied during the screening process of the data. The 
processor hardware and software development is in good progress and 
converted ERS-1 and simulated test data have been processed 
successfully. 

As an example for the SAR processing the D-PAF block diagram is 
shown . 


PAF Environment 



AA AnftragsabwickJang 

PI Principal Investigator 

JPO X-SaR Joint Project Office 

JPL Jet Propels ion Laboratory 

ASI Ag curia Spaiiak luliaoa 

HDDC High Density Digital Cassette 
SIR-C Shuttle Imaging Radar -C 


user copies 
delivery notes 
invokes 

l f accompanying 
sheets 
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polarisation 


Modes of operation 
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Global Wide swath, low spatial resolution, low data rate, HH and W 

Monitoring polarisations 
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MISSION LIFETIME 4 YEARS 

RELIABILITY .89 
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ASAR WIDE SWATH COVERAGE 
WITHIN 1 DAY 
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ASAR WIDE SWATH COVERAGE 
WITHIN 2 DAYS 
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ASAR WIDE SWATH COVERAGE 
WITHIN 3 DAYS 



96 


III 111 II 'll 1 1 li h III 



97 




RADARSAT Program 


$ 7-/8 _ 

/ 32% l /7 

^ i Q 

ihs } I 

N9 4 J 15 893 


J. McNally and S. Parashar 
Canadian Space Agency/Agence spatiale canadienne 
110 Rue O’Connor St, suite 200/piece 200 
Ottawa, Ontario K1A 1A1 
Canada 


ABSTRACT 


Work on the RADARSAT system is progressing well to meet the currently scheduled 
launch date of early 1995. The spacecraft bus and the Synthetic Aperture Radar (SAR) 
payload are at various stages of development. Requirements for the ground segment have 
been mostly established. The design of the ground elements such as mission control 
facility and SAR data processor is underway. The SAR applications development work 
is continuing and the chosen distributor, RADARSAT International Inc. (RSI) is making 
preparations to market RADARSAT data internationally. A plan for the follow-on to 
RADARSAT I is being finalized to ensure continuity of SAR data under the Radarsat 

program. 
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INTRODUCTION 


RADARSAT is a major Earth Observation (EO) program of the Canadian Space Agency 
(CSA). The program as presently approved includes construction, launch and operation 
of Canada’s first remote sensing satellite - RADARSAT I, and establishment and 
operation of associated ground facilities. RADARSAT I will cany one instrument or 
sensor, a Synthetic Aperture Radar (SAR) operating at C-band (5.3 GHj frequency) with 
horizontal transmit and horizontal receive (HH) polarizations. RADARSAT I is being 
designed to provide global SAR data operationally for five years after launch currently 
scheduled in early 1995. 

RADARSAT I will be launched into a dawn-dusk, sun-synchronous polar orbit by NASA 
using a Delta II rocket from the Western Test Range. In return for this contribution of 
the launch and related services, the USA government will receive for its use a pro-rata 
share of the SAR on-time available during the RADARSAT I mission. CSA will also be 
splitting the available SAR on-time with the private sector. A Canadian company, 
Radarsat International Inc. (RSI) has been assigned the worldwide distribution rights for 
the RADARSAT I SAR data. These rights are in return for financial contributions to the 
program such as for the ground SAR processing facility in Canada. RSI will develop the 
commercial data sales worldwide and pay revenues to the Canadian government to offset 
the cost of the RADARSAT I mission. 

One of the purposes of the RADARSAT program is to contribute to the creation of a 
viable international market for remote sensing data. This is to be accomplished by 
providing SAR data operationally through RSI to the worldwide remote sensing user 
community. Development of an international market by RSI, especially for operational 
applications, requires a continuing supply of SAR data. This need for SAR data 
continuity is now recognized and CSA accordingly is currently developing plans for 
follow-on RADARSAT missions. These plans should be finalized during 1993. 

RADARSAT I MISSION 

The RADARSAT mission is to supply and distribute global SAR data to the worldwide 
user community for such applications as ice mapping, ocean surveillance, agriculture and 
forestry monitoring, geological resource mapping and environmental studies. The 
development of applications and end users to utilize these applications and expand the 
economic benefits is an integral part of the mission. The Canada Centre for Remote 
Sensing (CCRS) has been managing a Radar Data Development Program (RDDP) for the 
last five years. This development work based on SAR data from the Canadian CV-580 
aircraft and European ERS- 1 satellite will be continuing so as to ensure users are capable 
and ready to utilize RADARSAT data. 
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To meet users needs the RADARSAT I SAR is to normally look right of the satellite 
track. This north looking configuration will provide essentially complete coverage of the 
Arctic but will leave a gap over the Antarctica. The spacecraft will be capable of 
undertaking manoeuvres so that the SAR can look to the left (south) for essentially a 
complete mapping of the Antarctic. This capability is being provided to meet the NASA 
requirement of mapping the Antarctic, once during a winter season and once during a 
summer season corresponding normally to the maximum and minimum formation of ice, 

respectively. 

The RADARSAT I SAR is designed to provide options to users in selecting swath-width, 
spatial resolution, and angle of incidence. As illustrated in Figure 1, these options will 
be provided through various SAR modes of operation selectable by the ground control. 
Within the accessability ground swath of 500km it will be possible to select individual 
beams or SCANSAR with the associated swath width and resolutions as given in 
Figure 1. Experimental coverage will be possible outside this range. The 500km 
SCANSAR swath will provide daily coverage of the Arctic, almost complete coverage of 
Canada and the USA each 3 days period, and global coverage over approximately 5 days. 
A 24 day exact repeat orbit has been selected for the RADARSAT I mission. 

In addition to the SCANSAR and multi-incidence beam observations, RADARSAT I is 
pioneering the dawn dusk orbit, i.e. the equator orbit crossing at approximately 6:00am 
or 6:00pm. This orbit offers advantages in solar power generation and associated design. 
It will avoid conflict in data reception on the ground with the other remote sensing 

satellites. 

SAR data will normally be acquired when RADARSAT is within view of one of the 
ground data receiving stations (Canadian Prince Albert and Gatineau stations, the USA 
Fairbank station, and stations as licensed by RSI). Two on-board tape recorders are being 
provided, however, to acquire SAR data for any part of the globe and dump the recorded 
data when the satellite is within the Canadian station mask. RADARSAT I will be 
controlled and programmed by Mission Control System (MCS) in Canada. The MCS will 
accept user requests from order desks in Canada, USA, and other regions for filling either 
from the archived data or newly acquired data. The SAR data will be processed using 
the processing facility being set up in Canada by RSI. The USA Fairbanks station is 
expected to have its own processing facility. The stations licensed by RSI to receive 
RADARSAT data will likely have their own processing facilities as well. The quality of 
the SAR data delivered from this RADARSAT system will be monitored and controlled 
by the MCS. To facilitate maintenance of SAR data quality a ground calibration site 
consisting of transporters and passive targets is being established in Canada. 
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RADARSAT I SYSTEM DEVELOPMENT STATUS 


The design and manufacture of flight hardware is well advanced at the prime contractor 
Spar Aerospace in Quebec and its major sub-contractors, COMDEV, CAL Corporation, 
MacDonald Dettwiler and Associates (MDA), SED Ltd., Fleet Aerospace and others 
across Canada. Ball Aerospace of Boulder, USA, the contractor for the spacecraft bus 
is in the final stages of assembling the flight hardware for the bus module. The other 
major sub-contractor is Domier who is producing the High Powered Microwave Circuit 
which is similar in design to the one used on the ERS-1. The on-board tape recorders 
being procured from Odetics Inc. are presently under test. 

Integration of the payload module (S AR sensor, tape recorders, and data handling system) 
is scheduled to begin mid 1993 at Spar. The final integration of the payload and bus 
modules plus die SAR antenna and the solar array will start early 1994 at CSA’s David 
Florida Laboratory in Ottawa. 

The existing Canadian receiving stations at Prince Albert and Gatineau, already 
operational for ERS-1, are being equipped to receive RADARSAT data. The existing 
Canadian facility for processing ERS-1 SAR data operationally is being upgraded for 
RADARSAT. The preliminary design review for this upgrade was held in 1992. The 
preliminary design reviews of mission control facility and telemetry, tracking and 
command stations of the MCS will be held soon. The requirement review for the mission 
planning and scheduling component of the MCS is planned for early 1993. The ERS-1 
transponder design will be used for the development of the ground calibration transponder 
for RADARSAT. The concept for order desks has been finalized. The work on 
operations planning, training, and utilizing development has been started. 


BEYOND RADARSAT I 


RADARSAT I should supply SAR data till the year 2000. The development of 
applications and associated utilization technology and the development of market for 
future sales of RADARSAT data are underway. To strengthen this development and 
commitment of users to invest in SAR data utilization, the users need to be assured of the 
continuity of RADARSAT SAR data. Accordingly, it has been realized that continuation 
of the RADARSAT program beyond RADARSAT I is vital. The Canadian Space Agency 
is presently developing the Long Term Space Plan (LTSP) for Canada. A plan for the 
continuation of the RADARSAT program is being developed as part of the LTSP. This 
plan is proposing that a second RADARSAT satellite be built as early as possible after 
the launch of RADARSAT I. This will provide a backup should there be a premature 
failure of RADARSAT I. If RADARSAT I operates for five years as anticipated, then 
the backup satellite will be launched as RADARSAT n around the year 2000. 
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In parallel with RADARSAT n, it is proposed that work on enabling technologies for 
RADARS AT n begins soon. This would entail advanced system and technology 
development for the SAR and key ground segment elements such as the data processor. 
It is also planned to support development of new SAR applications and value-added 
products. RADARSAT ffl is proposed to be ready for launch around the year 2004 so 
as to ensure continuity of the supply of SAR data under the RADARSAT program. As 
these SAR missions are costly and serve global user communities, mutually beneficial 
cooperation with interested countries or organizations will be explored as part of this plan 
for RADARSAT E and IE. 


Figure 1. RADARSAT SAR Modes 


Mode 

Width 

(km) 

Resolution 
R x A z (m) 

Looks 

Incidence 
Angie (deg) 

Standard 

100 

25x28 

4 

20-49 

Wide Swath 

(i) 

165 

48 - 30 x 28 

4 

20-31 


(2) 

150 

32-25 x 28 

4 

31*39 

Fine Resolution 

45 

11-9x9 

1 

37-48 

SC ANSAR 

(Narrow) 

305 

50x50 

2-4 

20-40 


(Wide) 

510 

100x100 

4-8 

20-49 

Extended 

(High) 

75 

22-19x28 

4 

50-60 


(Low) 

170 

63-28 x 28 

4 

10-23 
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SAR Operating Modes 
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RADARSAT Mission Concept 
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RADARSAT 1 Day Coverage 
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RADARSAT 1 Day Antarctica Coverage 
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RADARSAT 1 Day Antarctica Coverage 
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RADARSAT - General Layout 
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Abstract 

Results of conceptual design study and performances of 
developed key devices of BBM of the TRMM radar are presented. 
The radar , operating at 13.8 GHz and designed to meet TRMM 
mission objectives, has the minimum measurable rain rate of 0.5 
mm/h with a range resolution of 250 m, a horizontal resolution of 
about 4 km, and a swath width of 220 km. A 128-element active 
phased array system is adopted to achieve a contiguous scanning 
within the swath. The basic characteristics of BBM were confirmed 
by experiments. The development of EM has started with the 
cooperation of NASDA and CRL. 


1. Introduction 

TRMM(Tropical Rainfall Measuring Mission) is a joint space program of the 
USA and Japan(Ref. 1) to measure rainfall of tropics where about 60% of the global 
rainfall is concentrated. The feasibility study(Phase A) of TRMM was performed 
successfully under the joint efforts of the USA and Japan from January 1987 to 'March 
1988 and TRMM project was evaluated to be feasible. Phase B activities of TRMM were 
completed both in USA and Japan, and both countries are doing Phase C/D activities 
now. TRMM will adopt a circular orbit with the altitude of 350 km and the inclination 
of 35° . TRMM will be launched in August 1997 by the Japanese H-II rocket. 

The goals of TRMM are, (l)to advance the understanding of the global energy and 
water cycle by means of providing distributions of tropical rainfall, (2)to understan 
the mechanism through which the tropical rainfall affects the global circulation, and 
to improve the global circulation model, and (3)to evaluate satellite system for rainfall 
measurement. The primary mission product of TRMM is the monthly averaged 

rainfall over 5° x 5° grid boxes between the latitude of 37° N and 37° S with an 
accuracy of about 10% over three years. TRMM is required to understand tropical rain 
processes that play a key role in climate changes, particularly El Nino and Southern 
Oscillation. TRMM will contribute to international global change programs such as 
GEWEX. 

2. Conceptual Design Study of TRMM Precipitation Radar 

TRMM is the first space mission dedicated to measurements of tropical rainfall 
with the first precipitation radar in space. Communications Research Laboratory 
(CRL)performed a conceptual design study of TRMM precipitation radar in the 
feasibility study of TRMM in cooperation with NASA/GSFC (Ref. 2,3). The mission 
requirements for the TRMM precipitation radar basically specified by NASA are 
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summarized in Table 1. These requirements were used as the guideline for the design 
of TRMM precipitation radar. 

In the precipitation model, the rain is assumed to be uniform and extends to a 
height of 5 km. The thickness of the bright band or the melting layer is 0.5 km and the 
attenuation coefficient in the melting layer is assumed to be twice as large as that of 
rain below the melting layer. The following relations between the effective Z-factor 
Ze(mm 6 /m 3 ), attenuation coefficient A(dB/km) and rain rate R(mm/h) for 13.8 GHz 
are assumed: Ze = 372.4 R 1 - 54 (Ref. 4) and A = 0.032R 1124 (Ref. 5). 

In the TRMM precipitation radar, radar antenna beam is required to scan a swath 
width of 220 km in the cross track direction every 0.6 seconds in order to observe a 
raining area without any gaps between scanning lines which are perpendicular to the 
moving direction of the satellite. Therefore, high speed electric scanning becomes 
essential. 

Figure 1 shows the block diagram of the TRMM radar(Ref. 6). Active array radar 
is selected as a reliable candidate for the TRMM precipitation radar after trade off 
studies in the feasibility study. The system is reliable because it is still operational even 
if some parts of transmitters or receivers are damaged. The total system loss becomes 
small because transmission lines between the antenna and transmitter /receiver are 
short and loss of phase shifter can be compensated by the SSPA(Solid State Power 
Amplifiers) and LNA(Low Noise Amplifier). 

Main system parameters estimated in the feasibility study are shown in Table 2. 
The attained horizontal resolution is 4.3 km at nadir and the range resolution is 250 
m. The swath width is 220 km. The required number of independent samples of 64 is 
achieved by the frequency agility technique which uses dual frequency separated 6 
MHz. Dual frequency pulse will be transmitted mutually in the same pulse repetition 
period keeping a pulse repetition period long enough to insure the sufficient time for 
the data acquisition., The minimum detectable rain rate at the rain top is 0.5 mm/h 
and the measurable rain rate at the rain bottom is between 0.5 mm/h and 52 mm/h, 
with a signal to noise ratio of 0 dB for a single echo. If we use the attenuation data of 
the sea surface echo by rain, we can observe the average rain rate as large as 80 mm/h 
with a signal to noise ratio of 0 dB for a single echo. 

3. Development of Key Devices of the Bread Board Model of the TRMM 
Precipitation Radar 

CRL has developed key devices of the Bread Board Model of the TRMM 
precipitation radar since 1988(Ref. 7). These key devices are (1) slotted waveguide 
antenna elements, (2) 5-bit PIN diode phase shifters, (3) SSPAs and (4) LNAs. 
Integration of these components to form an 8-element Bread Board Model of the 
TRMM precipitation radar was also made. We developed 8-element slotted waveguide 
array antenna. Major required characteristics of the antenna are (1) the low sidelobe 
level to suppress the surface clutter, (2) the wide band width for the frequency agility, 
and (3)the high speed electric scanning with scan angle of ±17 degrees. We adopted the 
non-resonance type(wide band width)slotted wave guide antenna. The separation 
between slots cut on the narrow wall of the rectangular waveguide is 13.65 mm. Slot 
conductance is designed to attain the Taylor amplitude distribution(sidelobe = -35 dB, 
parameter n = 6) along the waveguide axis. Figure 2 shows outside appearance of the 
developed 8-element slotted waveguide array antenna part of the Bread Board Model 
of TRMM precipitation radar. The antenna length is 2.2 m. Measured antenna 
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Table 1. Mission Requirements 


Frequency 
Satellite altitude 
Scan angle 
Swath width 
Range resolution 
Horizontal resolution 
Minimum measurable 
rainfall rate 
Observation range 

Nadir 

Scan angle: 0 
Number of independent 
samples 

Antenna sidelobe level 


13.8 GHz 
350 km 
±17.0 degrees 
220 km 
250 m 

around 4 km at nadir 
0.5 mm/h at storm top 

15 km + 5 km 
15 km/cos0 
64 

around -30 dB 



* EM is developed by CRL 


Fig. 1 Block diagram of TRMM radar 



Fig. 2 Outside appearance of the 8-elenient slotted wave guide antenna part of 
the Bread Board Model of TRMM precipitation radar 
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radiation pattern in the plane of the wave guide axis of an 8-element slotted wave 
guide array antenna at 13.8 GHz is shown in Fig. 3 The half power beam width of 0.7° 
is attained. The peak sidelobe level about -30 dB is attained except for the shoulder 
level of -28 dB. The wide angle sidelobe level less than -35 dB is also attained. The 
antenna beam tilts about 4° from the direction normal to the slotted waveguide wall. 

We developed PIN diode 5 bit phase shifter at 13.8 GHz. In the PIN diode 5 bit 
phase shifter, 11. 25°, 22.5°, 45° and 90° bits are of loaded line type and 180° bit is of 
hybrid coupler type. We attained at the worst case between -10°C and +40°C, insertion 
loss values of 4-5 dB for the variation of the phase values. Phase errors about - 5°~ + 2° 
and amplitude errors less than 0.9 dB p . p are attained. As the phase and amplitude 
errors and variation of the insertion loss are directly related to the characteristics of 
phased array antenna, more detailed examination and evaluation of the phase shifter 
will be performed for the Engineering Model. Although values of the insertion loss 
are larger than those of the ferrite phase shifter, these loss values can be compensated 
by the SSPA and LNA. The mass of the developed phase shifter is about 150 g. 

In order to satisfy the peak transmitted power shown in Table 2, the output 
power of 10 W around the central part of the array is required. By using multi-stage 
GaAs FET amplifiers, output power of 10 W and an efficiency of more than 15 % is 
developed. In order to reduce the power consumption, the efficiency of SSPA itself, of 
the power supply and of the guard time margin for the pulsed operation must be 
considered. The output power variation is 0,2 dBp-p and the phase variation of the 
each SSPA is as large as 30° for the temperature variation between -10° C and +40° C. 
However, variation of the phase difference of 2 SSPAs whose phase difference is set to 
be equal at room temperature is less than 6° p . p for the temperature variation between 
-10° C and +40° C, because the temperature dependence of the phase variation of 2 
SSPAs is almost equal. The efficiency of the SSPA itself will be improved as large as 
25% in the development phase of the Engineering Model. The mass of the developed 
SSPA is about 400 g. 

LNA is composed of the multi-stage MIC amplifiers using the HEMT(High 
Electron Mobility Transistor). Good characteristics with gain of 30 dB and noise figure 
of 1.8 dB are attained. Variation of the phase difference of 2 LNAs whose phase 
difference is set to be equal at room temperature is less than 3° p . p for the temperature 
variation between -10° C and +40° C, because the temperature dependence of the phase 
variation of 2 LNAs is almost equal. It is important to attain the temperature stability 
of the phase of each SSPA and LNA, and also to keep uniform the temperature of the 
whole radar system in order to attain the stable antenna radiation pattern. 

Integration of PIN diode phase shifters, LNAs, phase adjusters and a 
divider/combiner is made to compose the 8-element transmitter and receiver part of 
the Bread Board Model of TRMM precipitation radar and the appearance is shown in 
Fig.4. The basic electric antenna beam scanning function by the switching of the 5-bit 
digital phase shifter was confirmed. The preliminary measurements of both the 
transmitting and receiving pattern of the 8-element active array antenna were also 
made to show good electric performances. 

4. CONCLUDING REMARKS 

Development of the Engineering Model of the radar has started with the 
cooperation of National Space Development Agency of Japan(NASDA) and the CRL. 
The Proto-Flight Model of the radar will be developed by NASDA, aiming at the 
launch of TRMM in August 1997 through joint efforts of USA and Japan. TRMM is 
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Table 2. Main System Parameters of TRMM 
Precipitation Radar 
(estimated in the feasibility study) 


Frequency 

13.8 GHz 

Antenna 


Type 

Slotted waveguide 
array(128 elements) 

Gain 

47.7 dB 

Beam width 

0.71 x 0.71 degrees 

Aperture 

2.2 m x 2.2 m 

Sidelobe level 

around -30 dB 

Scan angle 

±17.0 degrees 

Transmitter 

Type 

SSPA(x 128) 

Peak power 

577.8 W 

Pulse width 

1.67gs x 2ch 

PRF(#1) 

2778 Hz 

Receiver 

Noise figure 

2.3 dB 

IF frequency 

156 MHz, 162MHz 

Band width 

0.78 MHz x 2ch 

Characteristics 

logarithmic 

Dynamic range 

more than 70 dB 

Linearity 

less than ± 0.5 dB 

Smin 

-112.8 dBm 

Others 

Required S/N(#2) 

0 dB for single echo 

Total system loss 

2.0 dB 

Filter weighting loss 

1.5 dB 

Number of independent 

64(32 x 2) 

samples(#3) 

Data rate 

85 kbps 

Power consumption 

224 W 

Mass(#4) 

347 kg 


(#1) Fixed PRF; Frequency agility technique 

(e.g. 13.796 GHz and 13.802 GHz) is applied. 
(#2) 0.5 mm/h at storm top. 

(#3) Frequency agility technique is applied. 

(#4) Without margin 
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expected to improve the understanding of global atmospheric mechanism and water 

cycle, which is essential for human life. 
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ABSTRACT 

The data needs in scientific hydrology involve measurements 
of system states and fluxes. The microwave region is particularly 
well suited for measuring the system states of soil moisture and 
snow and the major flux into the earth as rainfall. This paper 
discusses the unique data needs for hydrology and presents some 
recent examples from AIRSAR experiments. 

INTRODUCTION 

Historically hydrology has developed as an engineering 
discipline to solve water resources problems such as flood 
protection and water supply. Evidence of the success of engineering 
hydrology can be found throughout well developed societies by their 
relatively high standards of living. Although there are still many 
water related problems to be solved worldwide, the realizations 
that water problems are no longer constrained to local drainage 
basins and the recent concern about climate change have asked 
completely new questions about hydrology; questions that 
traditional engineering hydrology is not equipped to answer. 

SCIENTIFIC HYDROLOGY 

According to Chahine (1992), the hydrologic cycle within the 
framework of climate change encompasses much more than the 
classical surface hydrologic framework (precipitation, evaporation, 
runoff, etc.). In addition to these land surface processes, scientific 
hydrology must focus on the interactive processes of clouds and 
radiation, precipitation, oceans, and atmospheric moisture. This 
conception of the hydrologic cycle not only addresses the transport 
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and storage of water in the global system, but also the energy 
needed and released through the phase changes. 

If we attempted to phrase these new concerns in the form of a 
single statement it would most likely be: "The major problem facing 
scientific hydrology is the tremendous spatial and temporal 
variability of hydrologic processes across the globe" (National 
Research Council, 1991). Thus it is primarily a scale effect that 
separates engineering hydrology from todays needs in scientific 
hydrology. However it is more than just a scale effect. We are now 
being asked more detailed questions about the intermediate stages 
of the hydrologic cycle in contrast to the engineering questions that 
have been frequently answered empirically. Thus in scientific 
hydrology we literally are asked where each drop of water resides 
and for how long and how this drop moves through the Earth's 
system. That is, we need to quantify fluxes and storages of both 
water and energy. 

ROLLOFSARS 

One of the more revolutionary aspects of SAR (Synthetic 
Aperture Radar) for hydrologists is the potential for measuring and 
monitoring various states of the hydrologic system. The major state 
variables that appear to be useful are the soil moisture, snow water 
content, snowpack condition, frozen soils. For the most part, 
hydrologists have modeled the hydrologic system pretty much as a 
"black box," using only input data (usually rainfall and maybe 
potential evaporation) and the output hydrograph. The unit 
hydrograph is a good example of a hydrologic "black box." Even 
though the development of the comprehensive hydrologic models 
exposed the interior of the black box and subdivided the rainfall- 
runoff process into a number of so called physical processes, this 
type of model was still pretty much a black box because there were 
no provisions for monitoring or measuring any system states or 
internal processes. 

Soil moisture is a system state that can be measured with 
microwave remote sensing. Techniques for measuring soil moisture 
include both the passive and active microwave approaches with each 
having distinct advantages. The theoretical basis for measuring soil 
moisture by microwave techniques is based on the large contrast 
between the dielectric properties of liquid water and dry soil. 

Future applications of soil moisture to hydrologic questions and 
applications are bound to become more common. Temporally 
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frequent spatial measurement of soil moisture will some day be 
available on a routine basis. Hydrologists are going to have to learn 
how to use these new data. In general, existing models have 
represented soil moisture in a way to make the model work but have 
not considered the possibility of independent determination of soil 
moisture or of soil parameters. For the most part, this approach has 
been justified because soil moisture data would not be available and 
hydrologists have not been able to deal with the spatial variability 
of soil moisture and soil properties. This may change with new 
microwave measurements. 

Snow, the amount and its condition, are important inputs to 
models the amount of snow in storage and its phase change to 
snowmelt runoff. Like soil moisture, microwave data appear very 
promising to the snow hydrologist. Not only can a microwave sensor 
be an all-weather instrument because it penetrates cloud cover, it 
can also penetrate the snow pack, which presents one with the 
opportunity of inferring many of the properties of the snow pack and 
the under-lying soil. These include depth and water content as well 
as the degree of ripeness, crystal size, and the presence of liquid 
water in a melting snow pack. As with soil moisture, the microwave 
measurement reflects several characteristics at once. 

A great deal of progress has been made in demonstrating our 
capability to measure these storage terms. All measurement 
approaches, from the laboratory, to truck, aircraft and satellite and 
shuttle platforms have been successful to some degree In addition 
there has been a great deal of increased understanding in the areas 
of theory, radar target interaction and algorithm development 

FUTURE SARS FOR HYDROLOGY 

The previous sections have discussed the basis of microwave 
remote sensing for soil moisture and snow and discussed their 
implications with respect to scientific hydrology. As promising as 
these new data seem, the future for using microwave data for 
scientific use is somewhat uncertain. For the next few years, 
researchers will be limited by the lack of suitable data. Only 
intensive and science-driven aircraft experiments will be available 
for collecting soil moisture and snow that hydrologists will find 
useful. Fortunately, there are a few experiments and satellites with 
SARs (ERS-1 and JERS-1) that should be invaluable for providing 
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sample data for developing and testing application models as well as 
answering some of the target-sensor questions. 

Looking ahead to when there may be microwave sensors on 
orbiting platforms, one confronts the basic differences between 
passive and active instruments and the intended use of the data. 
Comparing the instruments simplistically, the active sensors have 
the capability to provide high spatial resolution data (on the order of 
tens of meters) but their sensitivity to soil moisture may be 
confused more by roughness, topographic features, and vegetation 
than the passive systems. On the other hand, the passive systems, 
although less sensitive to target features, can provide spatial 
resolutions only on the order of tens of kilometers from a space 
platform. However, before the various potential microwave 
instruments can provide a stream of valuable data for scientific 
hydrology a number of questions must be addressed. 

RESEARCH NEEDED 

Although it appears that there will be more and more 
microwave measurements available in the future, there are a number 
of research questions that must be addressed before these data are 
available on a routine basis. The research questions can be split up 
into two categories for the purpose of discussion, one focusing on 
the microwave response and the other on hydrologic modeling. 

Microwave response: There are a number of unanswered questions 

regarding the microwave-target interactions that need to be 
answered before soil moisture and snow can be routinely determined 
with microwave instrumentation. There is a need to develop 
algorithms to abstract volumetric soil moisture directly from the 
microwave measurement (backscatter coefficient) . To do this, the 
other target characteristics of vegetation and surface roughness 
will have to be parameterized although recent work by Oh et al, 

(1992) indicate that polarimetric radar may provide a means for 
handling the roughness question. Connected directly to this need is a 
need to better understand the effects of surface roughness on the 
measured microwave response with respect to incidence angle, 
azimuth angle, wavelength, and polarization. Also, there is a need to 
understand the effect of the vegetation canopy on the microwave 
response. Vegetation variables include the geometry for the 
individual plant as well as the canopy as a whole, the water content 
(and perhaps the biochemical makeup) of the plant, and its stage of 
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growth. Microwave variables would include the incidence angle, the 
azimuth angle, wavelength, and polarization. 

There is a need to investigate the use of change detection 
algorithms for determining the relative soil moisture of an area and 
whether or not this information can be useful for hydrologists. 
Change detection should minimize the influence of target variables 
such as roughness and vegetation, at least over short time intervals. 
There is also a reasonable basis for expecting change detection 
methods to provide adequate data for hydrologic applications if the 
data are collected from a long term orbiting platform. Long term 
(multi-season or year) data will establish the upper (wet) and lower 
(dry) limits for the change algorithm. 

There is a need to develop software procedures for correcting 
the effects of terrain on the microwave response. Active microwave 
(SAR) is especially sensitive to this. This includes foreshortening, 
layover, and local incidence angle effects. Also, a potential issue is 
the relative accuracy of the DEM data with respect to the spatial 
resolution of the microwave data and the potential effect of 
subpixel variability on the measured signal. 

There is also a need to further investigate the potential for 
polarimetric SAR and its potential for abstracting target 
information such as the surface roughness and vegetation 
characteristics. Studies of this technique need to be carried out 
with carefully conceived ground data collection programs. 

HYDROLOGIC MODELING 

Because soil moisture has not been used as measured data in 
hydrologic modeling, there are a number of issues that hydrologists 
must address and solve in anticipation of actually working with 
spatial and temporal soil moisture. Existing hydrologic models and 
engineering techniques have been developed to solve specific water 
resource problems (i.e., water supply, flood protection, etc.). As 
such, many of these are lumped models that have been developed 
from point measurements. Models with this structure are not 
capable of using the spatial nature of remote sensing data. Before 
soil moisture data of this nature can be used in hydrology, there is a 
need to modify existing models or develop new models that reflect 
soil moisture in a way that is more physically realistic and so that 
soil moisture data can be used as input or to verify output. Also, 
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there is a need to develop procedures for modeling or estimating 
profile soil moisture from time series near surface measurements 
and/or using multifrequency data. 

Equally important to the research associated with modeling is 
a necessity to make the remote sensing data easily accessible and 
usable to the hydrologist. There is a need to develop software, 
perhaps within a work station concept, that can truly integrate GIS, 
remote sensing data from several sensors, and the hydrologic models 
so that all the calibration, scale, format, and so on, problems 
become transparent to the scientist and he can concentrate on the 
hydrologic analysis. 
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Introd uc tio n 

Regional air pollution, land-use conversion and projected climate change 
all affect ecosystem processes at large scales. Changes in vegetation cover 
and growth dynamics can impact the functioning of ecosystems, carbon 
fluxes, and climate (NRC, 1988; IGBP, 1990). As a result, there is a need to 
assess and monitor vegetation structure and function comprehensively at 
regional to global scales. 

To provide a test of our present understanding of how ecosystems operate at 
large scales we can compare model predictions of CO 2 * O 2 and methane 
exchange with the atmosphere against regional measurements of 
interannual variation in the atmospheric concentration of these gases 
(Mattthews and Fung, 1987; Tans et al., 1990). Recent advances in remote 
sensing of the Earth’s surface are beginning to provide methods for 
estimating important ecosystem variables at large scales. Ecologists 
attempting to generalize across landscapes have made extensive use of 
models and remote sensing technology (Tucker et al., 1986; Running et al., 
1989). The success of such ventures is dependent on merging insights and 
expertise from two distinct fields. Ecologists must provide the 
understanding of how well models emulate important biological variables 
and their interactions; experts in remote sensing must provide the 
biophysical interpretation of complex optical reflectance and radar 
backscatter data. 

Ecosystem models quantify the movement of carbon, water, and nutrients 
within and through terrestrial ecosystems and keep balance sheets on a 
daily, monthly, and annual basis. These models consider the source and fate 

of matter and energy that enter and leave ecosystems as a result of specified 
processes such as photosynthesis, respiration, decomposition, and the 
uptake and recycling of nutrients. 

Through the process of photosynthesis, leaves of plants capture CO 2 to make 
simple sugars while producing O 2 . Respiration returns CO 2 to the 
atmosphere using sugars in the growth process to construct complex 
organic compounds and to maintain living cells. Transpiration moves liquid 
water carrying dissolved nutrients from the soil through the vascular 

system of plants to the leaves where water evaporates into the atmosphere. 
Decomposition is performed by microbes that utilize organic substrates, 
recycling nutrients and releasing gaseous products such as methane and 
CO 2 to the atmosphere. To drive ecosystem models from space requires 

recognizing those structural features and environmental controls that 
affect the rates of key processes (the role of ecologists) and then finding 
ways of interpreting these variables remotely from space (the role of 

remote sensing experts). 



Because this is a SAR workshop, I will mention some special contributions 
that synthetic aperture radar may make in improving the estimates of 
fluxes from global-scale ecosystem models. My main theme, however, is to 
emphasize the value gained by combining information from a number of 
sensors to obtain confidence through independent estimates of critical 
variables and to predict relationships impossible to obtain from single 
sensors alone. 

Estimating Photosynthesis from Space 

The upper limit on photosynthesis, given a specified atmospheric CO2 level, 
is set by how much chlorophyll pigment can be packed into a landscape and 
by the amount of visible light intercepted by the chlorophyll. At any given 
latitude, on any given date, the solar radiation reaching the atmosphere is 
known. The challenge is to estimate the effects of clouds on attenuating the 
radiation through the atmosphere. Daily satellite-derived estimates of cloud 
cover can be made with a number of sensors. The TOMS (Total Ozone 
Mapping Spectrometer) satellite, which measures reflectance in ultraviolet 
waves lengths, provides monthly integrated estimate of visible radiation 
reaching the vegetation that correspondence directly with measurement 

recorded at ground stations (Eck and Dye, 1991). 

The amount of chlorophyll available to intercept light can also be estimated 
at monthly intervals from observing seasonal changes in the greenness of 
vegetation using the NOAA AVHRR series of weather satellites. The 
Normalized Difference Vegetation Index (NDVI), estimated from near- 
infrared and red reflectance measurements, is closely related to the 

chlorophyll content and the photosynthetic capacity of the vegetation 

(Yoder, 1992 ; Goward and Huemmrich, 1992). Adverse conditions for plant 
growth reduce this "greenness index" whereas improved conditions favor 
increases in the index in response to added synthesis of chlorophyll. 

Frequently, environmental conditions constrain photosynthesis below the 
potential established by the amount of visible light absorbed by 

chlorophyll. For example, stomatal pores close and halt all photosynthesis 
when leaves freeze at 0 to -2°C. Thermal sensors, coupled with 

measurements of NDVI, provide good estimates of ambient air temperatures 
on clear days when most extreme conditions occur (Fig. 1). Fundamental to 

predicting air temperature from satellites is the realization that (a) dense 

canopies of tall vegetation associated with maximum values of NDVI are 

within 2°C of ambient air and (b) that air temperatures can be estimated for 
less dense vegetation by linear extrapolation of the surface 

temperature/NDVI response to the maximum NDVI reference (Fig. 2). 
Because of a large change in dielectric constant when water freezes in the 
stems and branches of perennial plants, synthetic aperture radar provides 
an alternative means of evaluating freezing conditions (Way et al., 1990). 

In the unfrozen state, plant tissue varies in turgor and in water content. 
This variation in water status affects growth, photosynthesis, and 

transpiration in different ways. When water moves rapidly through the 

vascular system of plants some of the water columns in capillaries may 
break. As an adaptation, plants have evolved to close their stomata when 

flow rates approach critical levels during the day (Tyree and Sperry, 1988 ). 
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Because CO2 enters through stomata, closure of these pores reduces 
photosynthesis as well as transpiration. 

The driving force on transpiration and the movement of water through 
plants is the humidity deficit established between the water-saturated cells 
within leaves and the air immediately surrounding. The absolute humidity 
of the air can be estimated using two different thermal bands, one sensitive 
to water vapor in the atmosphere and one insensitive (Justice et al. , 1991). 
NOAA weather satellites (Advanced Very High Resolution Radiometer) 
provide the required combination of red, near-infrared, and two thermal 
bands to derive simultaneous estimates of the absolute humidity of the 
atmosphere along with estimates of ambient air temperature near the 
ground (Goward et al., 1993). With these two values the humidity deficit of 
the air around the vegetation can be calculated on clear days and a 
judgement made as to the extent that transpiration and photosynthesis are 
constrained below potential rates. 

Ecosystem models are particularly sensitive to the water status of vegetation 
at night when transpiration ceases and no (or very little) water uptake 
from the roots occurs. Under sustained drought, as progressively more 
water is extracted from the rooting zone, tension is placed on capillary 
water columns. Predawn measures of tensions on these water column allow 
prediction of how wide stomata can open the following day to permit 
photosynthesis and transpiration to proceed (Fig. 3). 

It is extremely difficult to estimate the rooting depth of vegetation and even 
more difficult to define the amount of water available, yet both of these 
variables plus daily precipitation are required by most ecosystem models to 
calculate a water balance and derive estimates of predawn plant water 
tensions. As an alternative, Goward et al. (1993) has demonstrated a 
correlation between the steepness of the Temperature/NDVI relationship 
(Fig. 2) and measured predawn water tensions (see also Nemani and 
Running, 1989; Goward and Hope, 1989). The relationship proved useful in 
estimating predawn water stress in a Douglas-fir/oak woodland in Oregon 
during summer months (Fig. 4). The relationship is empirical, however, as 
illustrated by lack of correlation at times when small amounts of 
precipitation affected surface wetness but did not reach the roots. 

Synthetic aperture radar, in principle, offers a more straight forward 
biophysical basis for estimating predawn tensions or a combined index that 
reflects major reductions in the hydraulic conductivity of branches and 
stems (Running and Waring, 1978). One would expect a reduction in the 
dielectric constant of the free water in conducting tissue, and accordingly, a 
reduced backscatter from L or P-Band radar (Narayanan and Vu, 1992). 
Diurnal variation in water tension is not nearly as useful in interpreting 

constraints of ecosystem processes (McDonald et al., 1992). On the other 

hand, if SAR can define progressive reduction in stem and branch water 
content during periods of high evaporative demand and restrictive soil 
water supply, it would provide a biophysically sound and more integrative 
basis for recognizing water stress on vegetation. 

Ecosystem models that include just those variables that constrain 
photosynthesis predict measured forest growth across a steep 
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environmental gradient in Oregon rather well (Fig. 6). In the analysis, the 
effective or utilized visible radiation ranged from over 90% of all that 
intercepted by maritime coastal forests (Sites la & Ih) to less than 25% of the 
radiation intercepted by juniper woodland (Site 6). 

Estimating Plant and Soil Respiration from Space 

Respiration associated with maintaining living cells increases with 
temperature and with the amount of protein in living biomass (Ryan, 1991). 
The leaves of plants contain the most active cells, the branches and stems 
contain progressively fewer live cells per unit weight. X-band (but not C- 
band) SAR shows potential for estimating the fresh weight of foliage, 
although seasonal and diurnal variation in leaf-water content may 
confound interpretations. Stem and branch biomass can be estimated with 

addition of L-band, but only up to levels of biomass in the range of 100-150 
metric tons of dry weight per hectare (Ranson and Sun, 1992; Dobson et al., 
1992). Other approaches using optical sensors offer a wider range of 

predictive abilities before saturating ( See Fig. 5; Li and Strahler, 1986). 

Soil respiration of C02 is also a function of temperature. Respiration from 
the surface litter decreases with drought; the drying of litter is related 
closely to the increased steepness of the Temperature/NDVI relationship 
described earlier. Energy balance equations can be applied to estimate soil 
temperatures and thermal sensors have been used to infer heat storage 
capacity (Luvall and Holbo, 1989). 

Methane production from soils is a special kind of microbial respiration 
particularly important in the tundra region of the Northern Hemisphere 
where a combination of standing water and hollow sedges offer a uniquely 
active environment for the generation of methane (L.A. Morrissey, NASA 
Ames Research Center; S.L. Durden, J.P.L., personal communication). The 
ability of SAR (C-band, VV polarization) to distinguish standing water and 
the upright stems of sedges provides an accurate measure of where methane 
production in the tundra is concentrated. SAR should offer a similar ability 
to identify subtropical and temperate rice paddies and other vegetation 
types that are periodically inundated and under such conditions produces 
methane in some abundance. 

Summary 

Ecosystem models predict many variables that can not be easily sampled 
over broad areas in a timely manner that reflects inter- and intra-annual 
variation in the rates of important processes. The net exchange of CO2 and 
other important trace gases, however, can be generated from ecosystem 
models. What is lacking is the means of initializing and driving these 
models at regional and global scale under continually changing conditions. 
To this end, satellites carrying a variety of remote sensing instruments 
appear the best means of testing the general reliability of regional-scale 
models. Some ecosystem variables require almost continuous monitoring 
such as cloud cover, as this affects the daily solar radiation reaching the 
Earth's surface. Other variables such as the "greenness" of the landscape 
may be recorded less frequently, dependent upon the degree of change 
seasonally. 
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There are advantages in synchronous coverage as shown in estimating 
vegetation greenness, ambient air temperature, surface wetness, and the 
humidity deficit of the atmosphere . By combining these series of remotely 
sensed observations limits on rates of photosynthesis, transpiration, 
respiration, nutrient uptake, and other ecosystem processes may be possible 
with present satellite technology. An evaluation of standing biomass and 
other structural characteristics of vegetation is required less frequently 
and therefore can be made at finer spatial resolution. Future satellites 
equipped with synthetic aperture radar (L and P bands) should be able to 
provide estimates of standing biomass below levels of about 150 metric 
tons/hectare. SAR appears ideal for recognizing some unique combinations 
of standing water and low stature vegetation associated with high rates of 
methane production in tundra and other landscapes. SAR also appears well 
suited to identify when vegetation is frozen. One of the most intriguing 
possibilities is that SAR technology may be able to provide a more direct 
biophysical measure of vegetation experiencing drought-stress than any 
other remote sensing system. This latter possibility warrants careful testing 
and documentation across a range of vegetation. 
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Figure Captions 

Fig. 1. Ambient air temperatures for twelve months were estimated from 
remotely sensed temperature measurements and vegetation 
"greenness obtained from NOAA AVHRR weather satellite for a series 
of sites ranging from temperate rainforest to juniper woodland 
(Goward et al. (1993). 

Fig. 2. On a 9 X 9 km grid, AVHRR thermal sensors were used to estimate 
surface temperature (Ts) on a drought-prone ponderosa pine forest 
with variable stocking density. When the 81 separate estimates of 
surface temperature were analyzed in conjunction with paired 
estimates of vegetation greenness (NDVI), the ambient air 

temperature could be estimated by extrapolating the relationship to 
the maximum NDVI, representing a canopy intercepting more than 
99% of all visible radiation (Goward et al., 1993). 

Fig. 3. Under static conditions of no transpiration at night, tension measured 
on the water columns of twigs (predawn water potential) predicts the 
extent that stomatal will open during the day (Waring and 
Schlesinger, 1985). 

Fig. 4. Seasonal comparison of predawn leaf water potential measurements 
and the NDVI/Ts slope were derived from near-infrared, red, and 
thermal channel measurements obtained from an ultralight aircraft 
flying at 300 m altitude. A steeping slope of the NDVI/Ts relationship 
correlates well with increasing drought stress measured on Douglas- 
fir trees, but some deviation occurs following small storms that only 
wet the surface litter (McCreight et al., 1993). 

Fig.5. Classification of shadow in video images of forests correlate well with 
forestry estimates of above-ground foliage and woody biomass up to 
levels near 700 Mg/ha (McCreight et al., 1993). 

Fig. 6.An ecosystem model that evaluates the light utilized by photosynthesis 
on an annual basis correlates well with observed above-ground 
growth of forests across a transect in Oregon. Only where forests 
have reached maximum height (Site 1, old-growth) was it necessary 
to take into account the increased respiration and reduced hydraulic 
conductivity of older trees to match observed growth (Runyon et al., 
1993). 
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I -INTRODUCTION 

On a global scale, forest represents most of the terrestrial standing biomass 
(80 to 90%). Thus, natural and anthropogenic change in forest covers can 
have major impacts not only on local ecosystems but also on global 
hydrologic, climatic and biogeochemical cycles that involve exchange of 
energy, water, carbon and other elements between the earth and 
atmosphere. Quantitative information on the state and dynamics of forest 
ecosystems and their interactions with the global cycles appear necessary to 
understand how the earth works as a natural system. 

The information required include the lateral and vertical distribution of 
forest cover, the estimates of standing biomass (woody and foliar volume), 
the phenological and environmental variations and disturbances 
(clearcutting, fires, flood) and the longer term variations following to 
deforestation (regeneration, successional stages). To this end, seasonal, 
annual and decadal information is necessary in order to separate the long 
term effects in the global ecosystem from short term seasonal and 
interannual variations. 

Optical remote sensing has been used until now to study the forest cover at 
local, regional to global scale. Radar remote sensing, which provides 
recently SAR data from space on a regular basis, represents an unique 
means of consistent monitoring with different time scales, at all latitudes 
and in any atmospheric conditions. Also, SAR data have shown potential to 
detect several forest parameters that cannot be inferred from optical data. 
The differences -and complementarity- lie in the penetration capabilities of 
SAR data and their sensitivity to dielectric and geometric properties of the 
canopy volume ; while optical data are sensitive to the chemical 
composition of the external foliar layer of the vegetation canopy. 
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Basically, the radar backscatter is sensitive to the quantity and state of water 
distributed in the complex vertical and horizontal structure of the canopy. 
Thus, it is expected that the radar measurement contains information on 
the forest type, tree architecture, tree biomass by parts (trunk, branches, 
leaves), canopy structure and also in certain conditions, the ground 
topography, soil moisture, soil surface roughness and understory 
vegetation. The physical relationships between radar backscatter and forest 
parameters depend on the radar frequency, polarization, incidence. 


Past results using SEASAT, SIR-A, SIR-B and airborne data demonstrated 
the use of SAR data for mapping of forest extent, mapping of wetland extent 
in forested areas, detection of disturbances due to deforestation, fires or 
flood. Some research attempts have been made to relate the radar 
backscatter to forest height, age class and biomass but the approach used was 
almost empirical. Given the complexity of the interaction mechanisms, to 
develop and validate algorithms which will provide the biophysical product 
from SAR data require a quantitative approach evolving both experiments 
and theory. 


II - RECENT ADVANCES 

During the past three years, significant progress has been observed in radar 
remote sensing of forest. There has been an emerging international 
concensus on the need for worldwide studies on processes associated with 
climate change, in which forest covers play an important role. Also, 
research in radar remote sensing has been driven by the numerous 
spacebome SAR missions scheduled in the 90's, starting with ERS-1 
launched in 1991 and J-ERS-1 in 1992. 

Technical advances in the development of SAR systems have provided in 
the past years multiparameter airborne SARs. Such systems, when 
calibrated, are optimum tools for algorithm development. Among those the 
current NASA/JPL system is an operational, three frequency polarimetric 
SAR which has taken part in many data acquisition campaigns in the world 
since 1988. Observation of forest covers is currently one of the most 
important research topics of such airborne missions. Experiments have been 
performed on a diversity of forest ecosystems in the US, Canada, Europe and 
in the tropics. Experimental procedures have been refined and more 
relevant ground data have been collected. Meanwhile, significant 
improvement in theoretical modeling of forest backscatter has been 
observed. This resulted from the availability and the relevance of 
experimental data and a better integration of the experimental and 
modeling approaches in the research projects. 
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Since 1991, ERS-1 has open a new era of SAR data that can be used on a 
regular basis. The radiometric and geometric quality of the data currently 
provided by ERS-1 makes possible the use of multitemporal observations 
either for terrain identification or for changes detection. 

Recent results in remote sensing of forests concern several aspects including 
forest type classification, determination of phenological and environmental 
states of the forests (Way et al. 1992) or the extent of flooded forests. 
However, the emphasis has been put on two main subjects : deforestation 
monitoring and biomass determination. 

n-1 DEFORESTATION 

The conversion of forests to agriculture and pasture and the harvest of 
forest for timber and fuelwood cause a net release of C02 and also other trace 
gases such as methane and nitrous oxide. Uncertainties in the estimates of 
the amount of carbon released globally from vegetation and soils result 
from different assumptions concerning 1/ the rate of land use change, in 
particular the rate of deforestation and 2/ the biotic response to disturbances. 


Even today although much progress has been made in recent years, 
estimates of the current rate of conversion of tropical forests to agricultural 
land vary from 70.000 to 100.000 km2 per year. 

The SAR is expected to contribute significantly to determine accurately the 
perturbation characteristics including the rate and location of 1 /ecosystem 
conversion to agriculture and pasture, 2/forests degraded for timber harvest 
of fuelwood and 3/ afforestation and natural regrowth. 

To measure the extent of deforestation has been well demonstrated with 
SEASAT, SIR- A, SIR-B and with airborne SARs operating at low frequencies 
(e.g. L and P bands). The contrast ratio between radar backscatter intensity of 
forests and other land uses including clear cut, agriculture, pasture is of the 
order of several dB in most cases. 

At higher frequencies (e.g. C and X bands), the radar backscatter of deforested 
areas varies as a function of the terrain characteristics including soil surface 
roughness, soil moisture, crop, grass or regrowth vegetation parameters. 
Thus, deforested areas can have a large range of backscatter values. With 
ERS-1 (5.3 GHz, 23° of incidence, VV) the conversion of forest to wide 
spreading pasture can logically be detected as shown on ERS-1 images in the 
Amazone basin (ESA, 1993). On the other hand, small patches deforested for 
shifting cultivation, e.g. in Africa, might be difficult to identify on the basis 
of their radar backscatter intensity at a single date. 
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The availability of reliable multitemporal data by ERS-1 offers possibilities 
to develop classification methods based on their temporal changes. In 
general, forest covers exhibit minor temporal change compared to dear cut, 
agriculture and other land uses. Through several examples of overlaid 
multitemporal images of ERS-1, the use of temporal changes for forest and 
land use classification and monitoring appear obvious. 

In the past few years, a new approach of dassification methods based on 
scattering behaviour of different surface types instead of their backscattering 
intensity have been developed ( Kong et al 1988, Van Zyl 1989, Evans et al 
1990, Freeman et al 1992). The development and validation of such methods 
have been possible since reliable airborne, polarimetric SAR data have been 
collected). In general, the polarization state of the received wave is 
compared to that of the transmitted wave in order to deduce the properties 
of the scatterer. The dassification can be unsupervised, based solely on 
comparing general properties of the Stokes parameters of the scattered wave 
to that of simple scattering models. As an example, results at L band ( Van 
Zyl 1989) indicated that scattering from dear cut areas and agricultural fields 
is mostly similar to that predicted by the odd number of reflection dass 
while the scattering from forest covers generally is classified as being a 
mixture of characteristics odd and even reflections and diffuse scattering. 

For future spacebome SARs, the use of polarimetric data may improve 
significantly the robustness of algorithms of forest/non forest delimitation. 

U- 2 FOREST BIOMASS 

Besides the rate of deforestation, one of the greatest uncertainties 
concerning the global carbon budget arise from a lack of information on 
forest biomass. First the assessment of the biomass of existing forests is 
necessary to determine the amount of CO 2 released to the atmosphere by 
forest burning. Also, to determine whether forests are a source or a sink of 
carbon, accurate estimate of biomass and carbon storage is needed, espedally 
that of regenerating forests which accumulate carbon quickly during the first 
20-30 years after disturbance. 

Optical remote sensing of forests is based on spectral reflectance data 
governed mainly by tree foliage properties induding leaf internal structure, 
leaf chemical composition and leaf angle orientation. With longer 
wavelengths, radar waves can penetrate into the forest canopy and the radar 
backscatter can give information on the trunk, branches, which represents 
90% of the total above ground biomass. 

In the past years, it was observed that the intensity in a SAR image at low 
frequency (L-band) is proportional to several forest parameters including 
tree height, stand age and above ground biomass of the forest stands. 
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However, to develop a robust algorithm for inversion of SAR data into 
forest biomass requires a thorough understanding of the interaction 
mechanisms between microwave and forest canopy. This is necessary to 
define the validity domains of inversion algorithm regarding the SAR 
parameters, the forest characteristics and the environmental conditions. 

The last few years have seen several experiments conducted on forest covers 
using multifrequency, polarimetric SARs together with an appropriate 
ground data collection plan ( Kasischke et al, 1991, Dobson et al. 1992 Ranson 
et al. 1992, Le Toan et al. 1992). To reduce uncertainties in measurements of 
forest parameters most experiments have been performed on coniferous 
forests, which present a relatively simple structure at the tree and stand 
level. 

The observations confirmed past results on the strong correlation between 
radar backscatter at low frequencies (P and L bands) and forest parameters. 
With well calibrated radars, quantitative comparison between frequencies 
and polarizations is now possible. P-band backscatter intensity yields the best 
correlation and the best sensitivity to forest biomass. Also cross polarization 
gives better results compared to like polarizations (Le Toan et al. 1992, 
Dobson et al. 1992 ). Polarimetric features, such as the polarization phase 
differences can be also related to forest age or biomass. However, their 
contribution compared to the backscatter intensity is not significant in terms 
of the retrieval of forest biomass. Nevertheless, polarimetric measurements 
have been very useful for the modeling of forest backscatter at different 
frequencies. 

Theoretical modeling will constitute the foundation of inversion 
techniques. The requirement is that the models must include all the 
physical properties of the target involved in the interaction process so that 
they can cope with simultaneous observations of a large variety of active 
(and passive) systems using multiple frequencies and multiple 
polarizations. 

In general, two main approaches have been used to solve the scattering 
problem : the field approach and the intensity approach. Also, 
approximations have been made to characterize the vegetation canopy 
either as a continuous medium or a discrete medium. 

Most models developed to describe the backscatter of forest canopies use the 
radiative transfer theory with a discrete medium approach. The models 
consider vegetation as a collection of lossy scatterers representing vegetation 
components (leaves, stems, branches, trunks). The forest canopy is in 
general modeled as a multilayer medium, each with a specific distribution 
of dimension, orientation and dielectric constants of the scatterers of 
various shapes. Although the radiative transfer theory provides only the 
incoherent terms, with the vector formulation of the theory multiple 
scattering problems and polarimetric backscattering can be accounted for. 
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Recent improvements in modeling of forest canopies included : 

a) the consideration of more complete structure of the trees instead of 
simplifications in the past assimilating the trees as leaves alone or branches 
and soil alone etc. 

b) more accurate calculation of the scattering of vegetation elements 
assimilated to finite dielectric cylinders or dielectric discs, 

c) the consideration of coherent interaction between elements of structured 
vegetation, 

d) the addition of second order scattering and, 

e) the consideration of the interaction between vegetation and rough soil 
surface. 

However, the most important progress concerned the model validation 
which benefited from more concertation between experimenters and 
modelers. More appropriate vegetation data have been measured in view of 
their transformation into model input parameters. On the other hand, 
experimental observations on the polarization and frequency behaviour of 
the radar basckscatter have led to substantial improvements of the 
theoretical models. 

Several results of model validation have been published recently ( Durden 
et al. 1989, Ulaby et al. 1990, Sun et al. 1991, Karam et al. 1992, Hsu et al 
1993). Once validated, the models can be used to identify different scattering 
mechanisms at different frequencies and polarizations. It was found for 
example that at P-band, the most significant mechanism between radar 
wave and a pine forest is HH polarization from the interaction between 
trunk and ground and for HV and VV, the direct scattering from the 
branches ( Hsu et al. 1993). Since there exists a strong correlation between 
different parts of the tree biomass, the correlation between radar 
backscattering coefficients and trunk biomass, branch biomass and total 
biomass as found experimentally are now understood ( Beaudoin et al. 1993). 

IH - CURRENT AND FUTURE WORKS 

More works are still needed to conduct experimental / theoretical studies on 
different forest ecosystems in order to generalize the results obtained on the 
retrieval of forest parameters. Also, the use of multifrequency, 
multipolarization data to retrieve different forest parameters including 
biomass by parts, canopy structure and ground information must be 
undertaken, at present with airborne SARs and in near future with SIR- 
C/X-SAR. 

On the other hand, the use of data from current satellites (ERS-1, J-ERS-1) is 
presently assessed for forest observations. 
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One possible research direction concerns a reexamination of the types of 
forest parameters of interest retrievable from SAR data in relation with 
ecosystem modeling. The sensitivity of radar backscatter to tree water 
content, tree water status and tree structure can be taken into account to 
derive some condensed descriptors of forest ecosystems. 
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RADAR RESPONSE OF VEGETATION: 
AN OVERVIEW N 9 4 



f 


Fawwaz T. Ulaby and M. Craig Dobson 
The University of Michigan 


• Vegetation Classes 

• Soil Scattering: (1) Backscatter 

(2) Forward Scattering 


• Radar Response 

• Vegetation Biomass 

• Vegetation Structure 

• Temporal Variations: (1) Short Term (hours to days) 

(2) Long Term (Seasonal) 

• Effect of Rain 

• Emergence of a User Community 

• Concluding Remarks 



The University of Michigan 
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STRUCUTRAL CLASSIFICATION OF VEGETATION FROM SAR 
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• Direct Ground Backscatter 

• Vegetation-Ground Bistatic Scattering 

• Trunks 

• Leaves (needles) 

• Branches 

• Direct Crown Backscatter 


University of Michigan 




MOD-3 

SOIL BACKSCATTERING 


A. Theoretical Models 

• Small Perturbation Model • Phase Perturbation Method 

• Physical Optics Model • Full Wave Model 

• Geometric Optics Model • Integral Equation Model 

Models agree with experimental observations only under 
certain conditions. Overall, models not useful. 


B. Michigan Empirical Model 


• Frequency Range : 1-10 GHz 

• Angular Range : 20° - 70° 

• Roughness range : s = 0.32 cm to s = 4.0 cm 
(expected validity for any s >0.3 cm) 

• Moisture range : 0.05 g/cm 3 to 0.31 g/cm 3 

Moisture Sampling Depth 

L-Band (1.25 GHz) : Average Moisture of Top 10 cm layer 
C-Band (5.3 GHz) : Average Moisture of Top 3 cm layer 
X-Band (9.5 GHz) : Average Moisture of Top 1 cm layer 





154 


Backscaitering Coefficient, o° (dB) Backscattering Coefficient, o° (dUB) 


MOD-3 


Mode l Verification For A smooth Surface (s =0.4 cm 


1.-5 GHz 



1 ' 1 » \ 

I ■*“ 1 

1.5 GHz 

5 = 0.4 cm 


o N 


O 

c w , Measured 

\ . 


s 

°°wi* Measured 

. 1 - 1 - • . « - t - 1 . ^ v . \ 


A 

o°^, Measured 

0. 10. 20. 30. 40. 50. 60. 70. 8 

Incidence Angle (Degrees) 

4.75 GHz 

0. 

▼ 

°°hr* Measured 

9.5 GHz 


4.75 GHz 


9.5 GHz 


10. 20. 30. 40. 50. 60. 70. 80. 
Incidence Angle (Degrees) 


0. 10. 20. 30. 40. 50. 60. 70. 80 

Incidence Angle (Degrees) 














/o°„ (dB) 
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Inversion Algorithm 

If radar measures cr^, and at a given frequency 
and angle, both s and m v can be determined from the ratios: 

P=<4./C'w 

Q ^hv l ^vv 


p-q Diagram (1.25 GHz, 40°) 



Note: p(dB) = lOlog^ / o£,) = cr^(dB)- <r° v (dB) 
q(dB) = <7hv(dB)- <7® v (dB). 
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ERS-1 SAR Response 
9 = 23 ° 

vv Polarization 





RADIATION LABORATORY 
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RADIATION LABORATORY 



Laser Profiler 



RADIATION LABORATORY 
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Reduction of Reflectivity by Surface Roughness at L-band 
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RADAR RESPONSE TO VEGETATION 
• OBJECTIVES 

• To Discriminate/Classify Vegetation Classes 

• To Estimate Biomass 

• To Estimate LAI 

• To Estimate Soil Moisture 

• To Monitor Changes (deforestation, growth, stress, etc.) 

• Other 




• VEGETATION CANOPY 

• Structure: (1) Macro (tree or plant scale): Tree height, 

density, ground cover 
(2) Micro (wavelength scale): Leaves, 
branches 

• Dielectric Properties 

• Ground Cover (soil, debris, undergrowth, etc.) 


• TOOLS 

• Wavelength 

• Polarizations 

• Phase Statistics 

• Incidence Angle 

• Time 


. APPROACH 

• Theory 

• Observations 

• Lab 

• Field 

• Air S AR 

• Satellite 


i 64 . The University of Michigan 
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ERS-1 RESULTS 


• Class Statistics 

• Observation versus Theory (MIMICS) 

• Biomass Response (Deciduous and Coniferous) 

• Seasonal Variation (LAI) 

• Deciduous 

• Coniferous 



The University of Michigan 
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Normalized Frequency 
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Observed Backseat ter Coefficient (dB) 


Comparison of SAIL Observations with MIMICS Simulations 
C-band, W-polarization 
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Dependence Upon Above Ground Forest Biomass 
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Backscatter vs. LAI for Deciduous Forests 
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Leaf Area Index (m A 2/m A 2) 


Backscatter vs. LAI for Closed-Canopy Conifer Forests 
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CONCLUSIONS & RECOMMENDATIONS 


1. Retrieval of Soil Moisture and 
Surface Roughness 

• |L-Band Quad-Pol for bare soil 


Status 

Field tested, 
Some AirSAR 
Verification 


• P-Band Quad-Pol: extends to agricultural crops 

2. Effects of Organic Debris 

Demonstrated 

• Extinction depends on size IX in lab 

At P and L-Bands, only trunks and large 
branches are significant 


1. In general o° = f (biomass, structure) 

2. Extinction by crown layer increases 

with frequency 

3. Scattering by foliage and small branches: 

• negligible at P and L Bands 

• dominates at C and X Bands 

4. Scattering by trunks and large branches: 

• dominates at P and L Bands 


MIMICS, 

AirSar, 

ERS-1 


The University of Michigan 
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CONCLUSIONS & RECOMMENDATIONS 

Even P-Band is insensitive to high biomass 
forests (Pacific NW = 500 tons/ha) 


6. Innundation under Forest Cover 


L-Band HH 


7. Effects of Intercepted Precipitation 

• negligible at P Band 

• = 1 dB increase or decrease at L-Band 


= 2 dB increase at C-Band 


SIR-B 

AirSAR 

Verified 


AirSAR, 

Scattering 

Verified 


8. Freezing of Vegetation Leads to 
Significant changes in c° at all Bands 


AirSAR, 

MIMICS 

Verified 


9. Deforestation Readily Detectable at 
P and |L-Band| 


SIR-B, 

AirSAR 


10. LAI Foliar Biomass Estimation 


C-Band Quad or X-Band 

11. Multi-Date Observations: Very Powerful Tool 


MIMICS, 

Field, 

AirSAR 


Requires good Relative Calibration (Stability) = ± 1 dB 


Requires good Absolute Calibration 


= ± 1 dB 


The University of Michigan 
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IMAGING RADAR STUDIES OF POLAR ICE 


Frank Carsey 
Jet Propulsion Laboratory 
California Institute of Technology 
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Topics 


• Scientific Overview 

• Radar Data Opportunities 

• Sea Ice Investigations 

• Ice Sheet Investigations 

• Conclusions 


Scientific Overview 


Sea Ice Scientific Objectives: 

• To estimate globally the surface brine generation, heat flux, 
and fresh water advection (as ice). 

• To monitor phasing of seasonal melt and freeze events and 
accurately estimate melt and growth rates. 

• To develop improved treatment of momentum transfer and 
ice mechanics in coupled air-sea-ice models. 


Kev Radar Observations: Ice Type and Velocity 
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STUDY PHILOSOPHY, UNDERSTANDING THE ROLE OF SEA ICE IN CLIMATE 
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SATELLITE DATA SETS FOR SEA ICE OBSERVATIONS 
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SEASONAL CHANGE IN ERS-1 BACKSCATTER 
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Radar Opportunities 


i 

/ 

• Systems Now Deployed in Space: 

ERS-1, JERS-1 

• Ground Station SAR Processors, Geophysical Processors 
International Collaboration--Could Improve 

• Systems Approved for the Near Term, Far Term 
RADARSAT (Coverage) 

• Proposed Systems for the Far Term 

• Airborne Systems, Polarimetric Data 

• In-Situ and Laboratory Scatterometers 

• Historical Data 

In Summary: Data Opportunities are Excellent. 

Will Funded Research Opportunities Keen Up? 


Scientific Overview 

Ice Sheet Scientific Objectives 

• To Map and Classify the Ice Sheets According to Dominant 
Processes in the Mass budget 

• To Monitor the Calving Ice Flux From Greenland and 
Antarctica 


Kev Radar Observations: Surface Conditions, Ice Velocity 
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GEOLOGIC REMOTE SENSING WITH RADAR 
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SURFACE AND SUBSURFACE MAPPING : landforms, 

PALEODRAINAGES 
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RADAR INTERFEROMETRY/ THERMAL IR 

LASER ALTIMETERS/STEREO IMAGING SPECTROMETERS 
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Develop and test the correlation via field experiments by analyzing wind data (zo), radar 
data (o°), and roughness data 


TOPSAR: HEKLA LAVA FLOW 
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ERS-1 image of eastern Libya showing ancient drainages, some of 
which may be buried, as evidenced by the bright linear sand dunes 
trending NW-SE over the dark drainages. North is approximately 
up and the image is about 100 km wide. 
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THE CLIMATE RECORD: HOW GLOBAL 
CHANGES IN CLIMATE ARE 
MANIFESTED LOCALLY 
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REGIONAL CHANGES IN THE REGIONAL ATMOSPHERIC 

CORRELATIONS LAND SURFACE CIRCULATION THROUGH 

THROUGHTIME TIME 
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MULTIFREQUENCY DUAL-POL. TOPOGRAPHY COVERAGE 
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Application of Cloude's target decomposition theorem to polarimetric imaging 

radar data 


Jakob J. van Zyl 

Jet Propulsion Laboratory 
California Institute of Technology 
4800 Oak Grove Drive, Pasadena, CA 91109 
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ABSTRACT 

In this paper we applied Cloude's decomposition to imaging radar polarimetry. We 
show in detail how the decomposition results can guide the interpretation of scattering 
from vegetated areas. For multifrequency polarimetric radar measurements of a clear- 
cut area, the decomposition leads us to conclude that the vegetation is probably thin 
compared to even the C-band radar wavelength of 6 cm. For a forested area, we notice 
an increased amount of even number of reflection scattering at P-band and L-band, 
probably the result of penetration through the coniferous canopy resulting in trunk- 
ground double reflection scattering. The scattering for the forested area is still 
dominated by scattering from randomly oriented cylinders, however. It is found that 
these cylinders are thicker than in the case of clear-cut areas, leading us to conclude 
that scattering from the branches probably dominates in this case. 


1. INTRODUCTION 

In the analysis of polarimetric imaging radar data, one is often faced with the situation 
where some geophysical parameters must be inferred from an area that exhibits 
significant natural variability in the scattering properties. In such cases, the resulting 
average Stokes or covariance matrix differs significantly from that of a single scatterer. 
It has long been speculated that if a unique way could be found to decompose such an 
average Stokes or covariance matrix into a sum of matrices representing single 
scatterers, one would not only be able to more accurately interpret the scattering 
processes, but the problem of inferring geophysical parameters from the measured 
radar data would also be significantly simplified. 

One of the first examples of such a target decomposition technique was provided by 
Chandrasekhar [1] in his book on radiative transfer. Considering the case of lateral 
scattering of light by small anisotropic particles, he decomposed the total average phase 
matrix into the sum of a phase matrix that represents dipole scattering and a phase 
matrix that represents pure random scattering, or noise. As such, his target 
decomposition followed the same principle that breaks the Stokes vector of a partially 
polarized wave into the sum of a Stokes vector representing a fully polarized wave, and 
a Stokes vector representing a completely unpolarized wave (see Papas [2] for example). 
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Huynen 131, in his thesis published in 1970, also introduced a target decomposition 
theorem in which he decomposed an average Mueller matrix into the sum of a Mueller 
matrix for single a scatterer and a “noise" or N-target Mueller matrix. However, unlike 
Chandrasekhar’s decomposition, Huynen’s N-taiget is not polarization independent, 
and has been shown to be just one of an infinite set of such residue matrices, meaning 
that Huynen's decomposition is not unique. 

In 1988 Cloude [4] introduced a target decomposition based on an eigenvector 

decomposition of the target covariance matrix. This decomposition was shown to be 

unique and, in the monostatic case, breaks the average covariance matrix up into the 

weighted sum of three covariance matrices representing three different single scatterers, 

analogous to the decomposition of a Stokes vector of a partially polarized wave into the f 

sum of the Stokes vectors of two full polarized waves that are orthogonally polarized [2]. 

In this paper we examine the application of Cloude's decomposition to the analysis of 

polarimetric synthetic aperture radar (SAR) data. P 


2. CLOUDE’S DECOMPOSITION 

Cloude showed that a general covariance matrix [T] can be decomposed as follows: 

4 

[TJ = X*< k £kJ (1) 

i=i 

In (1), K t :i = 1,2, 3,4 are the eigenvalues of the covariance matrix, k t ; i = 1,2, 3,4 are 
the eigenvectors of IT] , and kj means the adjoint (complex conjugate transposed) of 
k ( . In the monostatic case, for reciprocal media, the covariance matrix has one zero 
eigenvalue and the decomposition results in at most three covariance matrices on the 
right-hand side of (1). 

In (1) the elements of the eigenvectors are the elements of an equivalent scattering 
matrix in the basis In which the decomposition is performed. For example, if the 
covariance matrix elements are expressed in the normal HV basis, then the four 
elements of the eigenvectors are simply the four scattering matrix elements in the linear 
basis. Also, since the covariance matrix is, by definition, a hennitian matrix, it follows 
that the eigenvalues are all real, and the eigenvectors are orthogonal. Therefore, one 
can interpret (1) to mean that this decomposition breaks the covariance matrix into the 
weighted sum of covariance matrices from “orthogonal" scattering matrices. 

Cloude did his decomposition in the basis formed by the Pauli spin matrices, but 
pointed out that the decomposition can be done using any four complex matrices that 
satisfy the constraint of completeness and normalization. Here, we shall interpret our 
results in the natural linear basis in which the radar measurements are performed. 
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Also useful In our discussions later is Cloude's definition of target entropy, 

4 

Ht= -X p i l0 g4^ (2) 

t=l 

where 


Aj 4 A>2 + A3 4 * A4 

As pointed out by Cloude, the target entropy is a measure of target disorder, with 
H t = 1 for random targets and H T = 0 for simple (single) targets. 


3. DECOMPOSITION OF SCATTERING FROM AZIMUTHALLY SYMMETRIC NATURAL 

TERRAIN 


Before getting into the radar data analysis, let us illustrate the decomposition using a 
general description of the covariance matrix for azimuthally symmetrical natural terrain 
in the monostatic case. Borgeaud et aL [5] showed using a random medium model, 
that the average covariance matrix for azimuthally symmetrical terrain in the 
monostatic case has the general form 


where 


[TJ 


1 

0 

KP* 


0 

V 

0 


p) 

0 

cj 


c 


p 


( S hh S hh) 

(Shhfk) 

(Shh S hh) 


(4) 


(5) 

( 6 ) 


V = 



(7) 


c = 



( 8 ) 


The superscript * means complex conjugate, and all quantities are ensemble averages. 
The parameters C, 77 , £ and p all depend on the size, shape and electrical properties of 
the scatterers, as well as their statistical angular distribution. It is easily shown that 
the eigenvalues of [T] are 


209 



*1 = jc+i + ^{f -i) 2 + 4 |p| 2 | 

*2 = §{ c + i -^- 1 ) 2 + 4 > l 2 } 


A 3 = C TJ 

The corresponding three eigenvectors are 


4 = 


s 2 = 


I [C-1 + V£] 2 

r 2p/[C-l + VAp 
n 

|[C-i + Va] 2 + 4|p| 2 

u 

1 J 

[c-i-vsf 

to 
d 1 

H- * 
1 

|[C-i-Va] 2 +4|p| 2 

i J 


L 3 


ro^i 

1 

\°j 


In these expressions, we used the shorthand notation 
4 - [C - 1 ] 2 + 4|p | 2 


(9) 

( 10 ) 
(ID 


( 12 ) 


(13) 


(14) 


(15) 


We note that A is always positive. Also note that we can write the ratio of the first 
elements of the first two eigenvectors as 


fell 

1 [C - 1 + Va] 2 (jf • - 1 ■ - Va f + 4[p| 2 ] _ 1 _ Va ] 2 

ff\ 

1 

1 

[f-i-V4f([c-i + VJf + 4|p| 2 ] 


(16) 


which Is always negative. This means that the first two eigenvectors represent 
scattering matrices that can be interpreted In terms of odd and even numbers of 
reflections. In the rest of this paper, we shall examine decompositions of the form 
shown above for different types of terrain using measured polarimetric radar data. 
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4. EXAMPLES OF CLOUDE DECOMPOSITION} 


4 . 1 Randomly oriented dielectric cylinders. 


Before analyzing measured data, let us examine a special theoretical case of unifor ml y 
randomly oriented dielectric cylinders. In general, the scattering matrix of a single 
dielectric cylinder oriented horizontally can be written as 


IS] = 


a O' 

0 bj 


(17) 


where a and b are complex numbers whose magnitudes and phases are functions of 
cylinder dielectric constant, diameter and length (6). Generally, when the cylinder Is 
thin compared to the wavelength, |a| > |b|. To calculate the average covariance matrix 
for randomly oriented cylinders, one has to apply a rotation operator to [S] and average 
the results over ah angles of orientation. Assuming a uniform distribution In angles 
about the line of sight, one can easily show that the resulting average covariance matrix 
for the monostatic case is of the form given In (4) with 


c 

= [3|a| 2 +3|b| 2 + 29l(a*b)]/8 

(18) 

p 

[|a| 2 + |b| 2 + 69l(a*b)] 
[3|a| 2 +3|b| 2 + 29t(a*b)] 

(19) 

V 

2| a-b| 2 

[3|a| 2 + 3|b| 2 +29t(a*b)j 

(20) 

C 

= 1 

(21) 

It Is easily shown using (9) — (14) that 


h 

= C (i+W) 

(22) 

^2 

= c(l-|p|) 

(23) 


= C T] 

(24) 


The corresponding three eigenvectors are 


ki 


1 

V2 


W |pp 

0 


V 


) 


(25) 


211 


( 26 ) 


t 2 = 


1 

12 

1 


-P/\P\) 


0 

1 


(27) 


V 


0) 


From (19) and (20), it Is easily shown that 
l-|p| = 77 


(28) 


This means that for the case of randomly oriented cylinders, the second and third 
eigenvalues are the same. 


In the thin cylinder limit, |b| approaches zero, and we find that 

Pthin = 1 / 3 

Vthin = 2/3 


[T] = C 


r 

(\ 

0 

1 > 


1 
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- 1 > 


r 0 
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0^ 
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+ — 
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0 

0 
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+ — 

0 

1 
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3 



0, 



u 

0 

1 , 


l-l 

0 

1 , 


<0 

0 



(29) 

(30) 

(31) 


In this case, equal amounts of scattering are contributed by the matrix that resembles 
scattering by a sphere (or a trihedral comer reflector) and by the cross-polarized return, 
although a significant fraction of the total energy Is also contained in the third matrix 
resembling scattering by a dihedral comer reflector. The entropy In this case Is 0.95 
Indicating a high degree of target disorder or randomness. 


On the other hand. In the thick cylinder limit, |b) approaches |a| and we find that 


P thick = 1 


(32) 


thick ~ 0 


(33) 


In this case, only one eigenvalue is non-zero, and the average covariance matrix Is 
identical to that of a sphere (pr a trihedral comer reflector). Tlie entropy Is 0. Indicating 
no target randomness, even though we have calculated the average covariance matrix 
for ra ndo mly oriented thick cylinders! The explanation for this result lies In the fact 
that when the cylinders are thick, the single cylinder scattering matrix becomes the 
identity matrix, which is Insensitive to rotations. Therefore, even after rotations, we still 
added up only identity matrices, and no apparent target randomness Is introduced. 


For cylinders that are neither thick compared to the wavelength, nor thin compared to 
the wavelength, 0 < |b / a\ < 1. Figure 1 shows the behavior of the normalized 

the eigenvalues as predicted by (22)-(24). We notice, as mentioned before, that the 
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second and third eigenvalues are always the same magnitude, but they get smaller 
relative to the first eigenvalue as the ratio of the two scattering matrix elements 
approaches 1. 



Figure 1. Ratio of eigenvalues as a Junction of the ratio of b to a. The ratio of 
b to a depends on the cylinder dielectric constant and cylinder size compared 
to the wavelength. 


4.2 Measured scattering from clear-cut area. 

Using the NASA/JPL AIRSAR system [71, the complete Stokes matrix [8] can be 
measured for each pixel in a scene. The actual measurements in general compare well 
with the general form of the covariance matrix assumed In (8), except that the four 
terms assumed to be zero are not exactly equal to zero. However, these four terms are 
typically much smaller than the other terms, so that to a good approximation, the 
assumption In (8) Is valid. Here, we shall set those four terms to zero. Table 1 below 
shows the measured parameters at P-band, L-band and C-band for a clear-cut area in 
the Shasta Trinity National Forest in California. 


The clear-cut area is covered with short shrub-like vegetation. We note that at all three 
frequencies the scattering is dominated by an odd number of reflections. At P-band the 
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even number of reflections and the cross-polarized returns are almost the same 
strength, and about half that of the odd numbers of reflections. This is very similar to 
the thin randomly oriented cylinder case discussed earlier. As the frequency increeses, 
the even number of reflections and the cross-polarized returns become more different, 
and also become a smellier fraction of the total scattering. This is consistent with the 
randomly oriented cylinder case where the radius of the cylinder increases. The same 
conclusion is reached when considering the entropy. The highest value (most 
randomness) is observed at P-band, and the randomness decreases with Increasing 
frequency. As pointed out before, as the cylinder radius increases, the entropy 
decreases. Therefore, from the decomposition results we conclude that the vegetation in 
the clear-cut Is randomly oriented, and that most of the scattering comes from 
vegetation that has branches that are thin compared to the three radar wavelengths. 


Table 1. Measured ixilues for a clear-cut In the Shasta Trinity National 
Forest Also shown are the results of the Claude decomposition. 


Parameter 

P-band 

L-band 

C-band 

n 

0.5261 

0.5308 

0.4083 

C 

0.5642 

0.7580 

0.7159 

p 

0.0928+10.0582 

0.2324+10.1057 

0.3558+10.0440 


1.0260 

1.1615 

1.2437 

a 2 

0.5382 

0.5964 

0.4722 

a 3 

0.5261 

0.5308 

0.4083 

Entropy 

0.95 

0.94 

0.88 


4.3 Measured scattering from a forested area. 

Table 2 below shows the measured parameters at P-band. L-band and C-band for a 
forested area in the Shasta Trinity National Forest in California. We note that at all 
three frequencies the scattering from this coniferous forest is dominated by an odd 
number of reflections. At P-band the even number of reflections contribution is 
stronger relative to the cross-polarized returns when compared to the case of scattering 
by randomly oriented cylinders. This is likely a result of the P-band signals penetrating 
through the canopy and then suffering a double reflection when scattering from the 
ground and then from the tree trunks before returning to the radar. We notice that the 
imbalance of the even numbers of reflections and the cross-polarized returns decrease 
as the frequency increases, and also become a smaller fraction of the total scattering. A 
possible interpretation of this is that as the frequency increases, the penetration 
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through the canopy decreases and the scattering from the randomly oriented branches 
becomes more important. The highest value of the entropy (most randomness) is 


Table 2. Measured values for a forested area in the Shasta Trinity National 
Forest Also shown are the results of the Cloude decomposition. 


Parameter 

P-band 

L-band 

C-band 

V 

0.3301 

0.3485 

0.2416 

£ 

0.6529 

0.7122 

0.4685 

p 

0.2803+i0.0167 

0.3950+i0.0582 

0.3669+i0.0016 

Ai 

1.1566 

1.2805 

1.1873 

A 2 

0.4963 

0.4316 

0.2812 


0.3301 

0.3485 

0.2416 

Entropy 

0.87 

0.84 

0.75 


observed at P-band, and the randomness decreases with increasing frequency. From 
the decomposition results we conclude that the penetration through this coniferous 
canopy cannot be neglected. Most of the scattering can still be contributed to the 
randomly oriented branches, however. The relative strengths of the different 
mechanisms, as well as the lower values of the entropy, suggest that the branches are 
thicker than those of the vegetation in the clear-cut area. 


5. DISCUSSION 

In this paper we applied Cloude’s decomposition to imaging radar polarimetry. The 
decompositions illustrated here have been implemented on a pixel-by-pixel basis to 
analyze polarimetric radar images. Space does not permit the inclusion of more 
examples, however. 

We have shown how the decomposition results can guide the interpretation of scattering 
from vegetated areas. For clear-cut areas, we concluded that the vegetation is probably 
thin compared to even the C-band radar wavelength of 6 cm. For forested areas, we 
noticed an increased amount of even number of reflection scattering at P-band and L- 
band, probably the result of penetration through the coniferous canopy resulting in 
trunk-ground double reflection scattering. The scattering for the forested area is still 
dominated by scattering from randomly oriented cylinders. It is found that these 
cylinders are thicker than in the case of clear-cut areas, leading us to conclude that 
scattering from the branches probably dominates in this case. 
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The quantitative analysis of scattering requires detailed analysis of the decompositions 
that will result from model predictions. Future work will concentrate on analyzing the 
decompositions of model predicted covariance matrices. The results can then be used 
to analyze polarimetric radar images. 
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MULTIFREQUENCY OBSERVATIONS 

R. Glitz, F. Heel, H. Kietzmann, H. Ottl, and C. Schmullius 


1. Introduction 


N 9 4 


) 

ldV(^2 


Analyses of SAR image data were performed at DLR to classify various kinds of 
vegetation and different terrain types. The data were collected both with the 
DLR experimental synthetic aperture radar (E-SAR) in X-band, C-band and L- 
band and with the NASA/JPL DC-8 SAR in C-band, L-band and P-band. 

E-SAR is a single frequency and single polarization system (both parameters 
can be selected) but several flights were used to collect multifrequency/mul- 
tipolarization data which were geometrically matched after processing. Clas- 
sification of different crop types was based on comparison of the backscatter 
coefficients of calibrated SAR data in different frequency bands and polari- 
zations . 


The DC-8 SAR collects polarimetric data in different bands simultaneously. 
Data acquired with the NASA/JPL DC-8 SAR are qualified for scientific inve- 
stigations by reducing the cross-talk and channel imbalance to a tolerable 
extent. The data are absolutely calibrated by using reference targets with 
known backscattering cross-sections. The signatures and polarimetric features 
of terrain types, such as grassland, concrete, sea, forest (coniferous; deci- 
duous) and urban areas, are extracted and discussed with respect to frequency 
and incidence angle dependence. A multifrequency polarimetric feature vector 
was applied for classification. The results of this new approach for separa- 
ting and classifying different object classes are presented here. 


2. Crop Classification Based on E-SAR Data 

During the Multi-Sensor Airborne Campaign (MAC) Europe ' 91 on July 12th, the 
test site Oberpfaffenhofen was covered by E-SAR flights in order to collect 
X-HH, X-W and C-W SAR data. 

The latest capabilities of the DLR E-SAR are described in [1], and the re- 
sults of the MAC '91 agricultural classification in [2]. The high spatial 
resolution of E-SAR (2 to 3 m) allowed even the analysis of small fields with 
very good results. An unsupervised migrating means classification (k-distri- 
bution) resulted in a rejection of only about 5 % of the observed area. 

On May 20th, 1992, another airborne campaign was performed with the E-SAR, 
after it had been modified to provide an additional L-band capability. L-HH, 
C-W, X-W, and X-HH data were collected. A selected agricultural test site 
is shown in Fig. 1 corresponding to these four data sets. 

It should be mentioned that the spatial resolution of the L-HH image is lower 
than in C-band and X-band images. The approximately equally spaced white 
spots in the images represent electric pylons. Test site statistics of the 
six main crop types are illustrated as scatter grains in Fig. 2 for C-W and 
X-W and in Fig. 3 for L-W and C-W. The ellipses include +1 standard devia- 
tion around the mean. 

The plots in Fig. 2 permit a proper discrimination of the various crop types, 
while Fig. 3 only shows small differences in L-band for a number of crops. 
This observation is valid for the vegetation in late May. 
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Fig. 1: 


Top: The SAR images based on X-HH, X-VV, C-VV and L-HH data. 

Bottom. Color composite image consisting of the following co-registered channels: 
L-HH (red), X-VV (green), and C-W (blue). 
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Scattergram, E-SAR 20. May '92 

Fig. 2 : Scattergram of C-Wversus X-W bands 

showing proper crop type discrimination. 
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Fig. 3: Scattergram of C-Wversus L-HH bands. 

showing overlapping crop type statistics 
in the L-band. 


220 



3. Airborne Campaign with the NASA/JPL DC-8 SAR 


In August 1989, an airborne SAR measurement campaign took place at the Ober- 
pfaffenhofen (OP) test site area. The NASA/JPL DC-8 SAR system is equipped 
with three simultaneously operating SAR sensors, comprising C-, L- and P-band 
with wavelengths of 5.6 cm, 24 cm and 68 cm [3]. This system permits fully 
polarimetric measurements [4] which means acquiring the scattering matrix 
elements HH, HV, VH and W. Three parallel tracks were flown to cover scenes 
of the same area under different incidence angles of 35°, 45°, and 55°. 

In order to qualify these multi frequency data products for comprehensive sig- 
nature analysis, the appropriate preprocessing had to be carried out [5] . 

The data prepared in this way, with their diversity in frequency, incidence 
angle and polarization states, form a good basis for signature and classifi- 
cation studies. 

4. Data Preparation and Terrain Type Selection 

JPL delivered single look complex SAR image data. These fully polarimetric 
data covered the P-, L- and C-frequency bands and had been radiometrically 
corrected. 

As a first step to qualify all the data products for signature analysis, po- 
larimetric calibration (cross-talk, channel imbalance) was performed at the 
Institute for Radiofrequency Technology followed by absolute calibration [6], 
[7]. Finally, the data were transformed to a 4-look image format. 

After adequate image data preparation, a total power C-band image ( Fig. 4 ) re- 
presenting the central part of the Oberpfaffenhofen test site area (12.5 km x 
5 km) was selected for terrain type investigations. This image includes the 
contours of the specific terrain types to be analysed, such as grassland (Wl, 
W2), concrete (Bl, B2), lake water surface (SI, S2), coniferous forest (Nl, 
N2), deciduous forest (LI, L2) and urban areas (Al, A2) . The terrain type 
sizes vary from about 10 000 m 2 to 30 000 m 2 depending on the availability of 
the corresponding terrain types. 


5. Multifrequency Terrain Type Signatures and Classification 

The simultaneous P-, L- and C-band measurements of each selected terrain type 
enabled the study of their backscattering behaviour. Fig. 5 shows the compa- 
rison of their copolarized backscattering coefficients Ohh- It is generally 
seen that urban areas and forests give rise to considerably higher returns 
than grass, concrete or lake targets. This can be explained by an increased 
roughness and reflection-like interactions. As expected for surface targets, 
especially in the case of grassland, a a^h decrease could be found in the 
C-L-P sequence. An inversion of this sequence can be recognized for all fo- 
rest areas. C-band signals mainly interact with needles or leaves in the up- 
per forest region. Medium sized twigs and branches predominantly scatter the 
L-band signals, and the P-band signals couple essentially with the larger 
branches and trunks. This demonstrates that the size of the scattering ele- 
ments considerably influences the target return. The urban area data show no 
significant wavelength dependence because of similar scattering mechanisms 
for the aforementioned frequencies. 

The corresponding copolar phase differences are shown in Fig. 6 . We recog- 
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Figure 4: Total power C-band image with selected terrain types 




nize a phase difference for urban araas of approximately jfl80 , indicating' a 
relatively clear double bounce scattering in contrast to all analysed surface 
type targets. In the case of concrete and lake areas, there are some phase 
anomalies due to interferences from stronger targets. Regarding forests, the 
phase increases with increasing wavelength due to scattering mechanisms, as 
shown in Fig. 6 . 

Further investigations deal with the influence of the incidence angle on the 
backscattering coefficient Ch h of coniferous and deciduous (Fig. .7) forests. 
For coniferous forest, a clear decrease in Ohh is observed in the incidence 
angle range from 40° to 50° for all concerned wavelengths, whereas deciduous 
forest shows this dependence only for C-band. This latter target type reveals 
no obvious incidence angle variation of Ohh for.P- and L-band which could be 
explained by the orientation of branches and twigs in this angle range. The 
Ohh behaves like that of a rough surface consisting of a layer of needles or 
leaves. The bars in Fig. 7 indicate the standard deviation around the mean 

values . 

A statistical procedure [5] has been developed, based on the maximum-likeli- 
hood method, to take into account the ramdomness of distributed targets. 

Analyses of a variety of polarimetric quantities allow the extraction of the 
most significant components to establish a multifrequency polarimetric fea 
ture vector. This feature vector consists of the scattering amplitudes HH, 

W, HV and the HH/W phase difference for the applied P-, L- and C-band fre- 
quencies resulting in a total of twelve dimensions. In order to meet the sta- 
tistical nature of each terrain class of interest, the mean and the covarian- 
ce have to be determined for the single features and the overall feature vec- 
tor. Adequate processing of this feature representation permits a separation 
of different terrain classes depending on their specific scattering proper- 
ties. Thus, a procedure is established for classifying different terrain ty 
pes. Fiq. 8 illustrates a classification result extracted from the scene in 
Fig. 4 concerning coniferous forests. The black areas in the scene are clas 
sified as coniferous forest. The corresponding test areas Nland N2 are cor- 
rectly assigned whereas negligible misclassification occurs in all other test 
areas. Results with similar quality have also been obtained for other terrain 

classes. 


6. Conclusions 

The comparison with ground truth observations shows that terrain classifica- 
tion from calibrated radar image data in different frequency bands exhibits a 
high degree of reliability. The amplitude relationship of backscattered sig- 
nals in X-band and C-band contains sufficient information for crop classifi- 
cation. 

In general, for the classification of different terrain types, the use of a 
multifrequency polarimetric feature vector showed promising results. The 
large number of degrees of freedom associated with this vector allows an in- 
creased variety of signature investigations. 

Therefore, to enhance the potential of remote sensing, such methods should be 
further investigated and applied in the future. 






Figure 8: Coniferous forest classification image 
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THE NEED FOR RADAR SIGNATURE MEASUREMENTS 


Alois J. Sieber, Carlo Lavalle 
Institute for Remote Sensing Applications 
Joint Research Centre 
1-21020 Ispra - Italy 


Parameters which are characteristic for selected objects are known as signatures for these objects. As 
an example one may think about the handwriting as being such a characteristic signature for the indi- 
vidual writer. One important task in working with radar data is to search for information about targets 
of interest which are unique for these objects. However, previous work with air- and spacebome 
radar images demonstrated that the resulting object information does not only depend on the prime 
measurement parameters like the frequency, the polarization and the illumination geometry adopted 
in the measurement. Further important effects are due to the illumination technique, the length of the 
synthetic aperture or the changing illumination angle throughout the synthetic aperture, the band- 
width and,as a consequence, it seems that each system, air or space-borne, shows different object 
features. 

The upcoming two years will provide for the first time unique opportunities to compare target infor- 
mation collected almost simultaneously over the identical sites by SAR systems on the successfully 
working ERS-1 of ESA, the Japanese ERS-1, the unique multi-temporal Shuttle Imaging Radar 
(SIR) programme, the NASA/JPL multiparameter AIRSAR, the Canadian multi parameter SAR 
sensor of the CCRS and INTERA, the European multi-parameters SAR sensor EARSEC as well as 
many ground based radar systems in different parts of the covered part of the earth. All these data 
collection efforts will be complemented by a device dedicated to radar signature measurements un- 
der fully controlled environmental conditions, the European Microwave Signature Laboratory 
(EMSL). 

The paper will present a systematic approach to investigate the signatures of selected objects using 
controlled experiments and translating them into the real world by physically sound models. 

While in-field experiments are locally fixed, while air-borne measurements are regionally bound 
and while space-borne systems are often discussed in the frame of global issues, the multi-parame- 
ters Shuttle Imaging Radar programme is a perfect tool for the data collection of geophysically cali- 
brated radar signature on a large scale - if the interpretation of its data is physically based and sup- 
ported by controlled signature measurements. The paper will demonstrate that the Shuttle Imaging 
Radar (SIR) programme and the work in the European Microwave Signature Laboratory (EMSL) 
are complementing each other and are linked by the North American and the European airborne SAR 
systems. 
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SPACEBORNE IMAGING RADAR 



ANTENNA TRUNNION STRUCTURE (ATS) 

X-SAR SUPPORT STRUCTURE (XSS) & TRI-DRIVE 




SIR-C/X-SAR FREE FLYER CONCEPT - OVERVIEW 
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- MUST BE UNDERSTOOD AS IT AFFECTS ANY OPERATIONAL 
ALGORITHM 

MEANS OF ASSESSING TRADEOFFS BETWEEN SEASONAL 
MULTITEMPORAL OBSERVATIONS AND FREQUENCY/POLARIZATION 
DIVERSITY 


GEOPHYSICAL PARAMETERS FROM SAR 


DC 

< 

CO 

CO 

o 

U! 


DC 

UJ 

> 

-J 

u. 

• 

UJ 

UJ 

DC 

Li. 


9 < 

gco> 

<n>k 


CO 

C 

2 < 

DC H DC 

<<£ 

|Si 

OH< 

DCZ 0 - 

— U1LJJ 

< cc -J 
DC O 
D 2 
Oco 



c 

© 

4_ X 

C UJ 

>. © c 
(5 x £ _ © 

O “ Si 2 

^ ! §>! I 

§i?-M 

> 3 O O 

x c <0 C 
' - w CO 


CP) 

(0 

o. E <o 

> O s 

> H 5 c 

05 r S 

Xo £ o 
O — > Q_ 
-j2 o t 

X o) o t; 
U ® o 5 
Lii > a 5 



c 

o 


T 5 

O 


o 

3 

55 

© ® 
o o 

t 1 

3 © 


si 


O 


a 

< _ 

DC © 

a ® 

O © 

VL 

LU © 

OS 

oo 


o 

3 

"D 

O 


« 

o 

3 

■D 

O 


© 

V © 
co " 

>N 

x: 

CL 
O 

© CL 

a z 
_ © 
© o 

8 £ 

S | 

to § 

« o 
8" 3 

2 o 

CL 

o E 

j= © 

cn H 

Q. -L 
© ~ 

¥> I 


© 

-Q 

O 


© 

a 

o 


O O 


1 1 


c 

o 


© 

*D 


238 


PROGRAMMATIC GOALS 
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CONCEPT DESCRIPTION 
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SIR-C/X-SAR FREE FLYER CONSTRAINTS AND ASSUMPTIONS 
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COST IS THE PRINCIPAL DETERMINANT OF FEASIBILITY 



A CONFIGURATION CONCEPT 
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Deployed Solar Array 




Deployed Solar Array — 
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Electronics Boxes 



Thruster Cluster 



Deployed Radiator 









A CANDIDATE MISSION DESCRIPTION 
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AN IMPLEMENTATION SCENARIO 
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SIGNIFICANT SCIENCE VALUE. BEGINS EOS SAR GLOBAL DATA SET 
EARLY 
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Earth Observing System Goals: 
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initially with Germany 

currently exploring broader international partnership 



Geophysical Properties 
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Surface roughness 

Subsurface structure and drainage 
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Operational Geophysical Products from SAR 
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Snow moisture and extent L-quad, C-quad, X weekly, day-time 

Vegetation water potential/moisture L-quad, C-quad, X diurnal 

Thin ice type L-quad ? 3 days 

Inundation extent (forests) L-HH weekly 



EOS SAR Mission Strategy: MultiSAR 
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change in "classification" properties 


Split Concept for Global Data Products 
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Potential Synergism of International MultiSAR Mission 
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diurnal coverage 
increased mapping 
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JPL POLAR REGION APPLICATIONS 
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• ICE FLOW IS RELATED TO SURFACE HEIGHT AND THICKNESS 

• REPEAT, HIGH RESOLUTION ELEVATION DATA WOULD ALLOW MONITORING OF 
ICE TRANSPORT AND ABLATION 
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application Planimetric Planimetric Altimetric Extent Site Repetitevity Vegetation Sensor Mission 

Resolution Error A/P error A/P Q/R/L E/A/SS V/G I(m)/N Lifetime 
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Obtain complete global coverage in less than 6 months and 
continue measurement capability for at least 12 months in 
order to monitor seasonal and shorter period changes of ice, 
vegetation, wetlands, and time-varying landforms. 
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SPACE-BASED STEREOSCOPIC 

MISSIONS 
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GLOBAL TOPOGRAPHY MISSION 

OPTICAL STEREO APPROACH 
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• "MISSION” COSTS > $500 M 

• VERTICAL ACCURACY DOES NOT MEET SCIENCE REQUIREMENTS 

• COVERAGE NOT GLOBAL 

• DIFFICULT TO ENSURE UNIFORM QUALITY, C0NTIGNU0US DATA SET 
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RADAR INTERFEROMETRY 

THEORY 
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5 = §\I2k 

sin(0-a) = ((p+5)**2 - p**2 -B**2)/(2*p*B) 
z(x) = h - p cos(0) 


TOPSAT SATELLITE IMPLEMENTATIONS 
REPEAT PASS / SINGLE / TWIN SATELLITES 
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ALONG TRACK INTERFEROMETRY 
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In the case of along-track interferometry, only the line of site velocity can be measured and 
therefore the along track velocity component is unknown. 



ON-ORBIT CONFIGURATION 




TOPSAT - VISTA 
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Number of looks: 8 azimuth / 1 range 
Height accuracy: ~3 m 



VISTA Advantages 
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Larger antennas are more 
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DUAL SPACECRAFT CONCEPT 

VIEW FROM BEHIND VELOCITY VECTOR 



RADAR SWATH 
3S KM 



GEOMETRY NEEDED TO OBTAIN NEAR POLAR COVERAGE 
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L-band Technology Assessment Summary 
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Calibration Subsystem 




TWO S/C / DELTA II 
LAUNCH CONFIGURATION 
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TOPOGRAPHIC MAPPING LASER ALTIMETER 
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GLOBAL TOPOGRAPHY MISSION 
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SLOPE MEASUREMENT ACCURACY at the 1 0 LEVEL 
ALONG-TRACK AND ACROSS-TRACK SLOPES 



ROLE OF MULTI-BEAM LASER ALTIMETER 
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HIGH-ACCURACY CONNECTION OF LAND AND OCEAN TOPOGRAPHY 

< 20 cm SURFACE HEIGHT ERRORS. IN COASTAL REGIONS 

< 1 0 cm WITH MULTI-PULSE AVERAGING 


S20- *7 


SOME THOUGHTS ON FUTURE RAIN MAPPING MISSIONS 
(TRMM FOLLOW-ON) _ . 

N 9 4 “ 


John S. Theon 


1 5 9 0 ti 


National Aeronautics and Space Administration 
Washington, D.C. 20546 


INTRODUCTION 

The release of latent heat of condensation is the largest internal energy source of the atmosphere. 
Latent heating is most significant during the precipitation process. Our knowledge of the 
distribution of precipitation is poor. It is only well observed within limited areas of the globe. 
Over the oceans, for example, it is known only to within a factor of two. Thus, there are strong 
scientific requirements for observations of precipitation from instruments on board a satellite. 
The Tropical Rainfall Measuring Mission (TRMM) will be the first satellite to measure rainfall 
with adequate accuracy and provide information about the vertical distribution of precipitation, 
not in only in a tropics and sub tropics, but to plus or minus 35 degrees of latitude. There is a 
need for the continuity of rain observations for climate modeling purposes and to expand 
observations to cover more of the globe than simply plus or minus 35 degrees latitude. 

TRMM FOLLOW-ON MISSION REQUIREMENTS 

One of the requirements that can be met by a TRMM follow-on mission is the need for five 
additional years of tropical rainfall measurements for short term climate studies. A three year 
data set simply may not include both normal and El Nino Southern Oscillation (ENSO) 
conditions. ENSO occurs about twice in a decade, maybe three times, but we need to extend the 
data set to include El Nino; La Nina which is a cold eastern pacific oscillation; and normal 
conditions. A total of eight years should achieve this goal. Shukla from the University of 
Maryland has shown that the results of large scale climate models are immensely sensitive to the 
amount and vertical profiles of latent heating in the tropics. The only way to obtain adequate 
quantitative values for latent heating is from TRMM and its hydrometeor profiling. Shukla also 
showed that the impacts of tropical circulations on mid latitude weather and climate requires a 
much better knowledge of tropical latent heat release than we have at present. 

We also have a need to extend measurements to middle and high latitudes both for analysis of the 
global hydrological cycle and to understand precipitation and its impact on the environment. 
Measurements are urgently needed for the Global Energy and Water Cycles Experiment 
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(GEWEX). We need to begin the measurement of frozen precipitation from space. Over land, 
only radar can make these measurements. Even radar will have problems detecting light snow 
fall. The GEWEX Continental Scale Experiment (GCIP) will greatly help testing and validating 
snow retrievals over land. By the year 2000 we need rainfall measurements in the tropics and 
mid latitudes on a near real-time basis to assimilate in large scale forecast models. This means 
rain and latent heat data simulations will have been developed during TRMM-1. This could lead 
to significant increases in accuracy for one to five day weather forecasts and greatly improved 
assessment and prediction of global change, particularly global warming, by improving the 
climate models. We need better models which couple the land surface hydrology with 
precipitation. Most hydrological studies are on mid latitude sites. (Personal communication 
from J. Simpson, GSFC, January 1993 


Shown in Figure 1 are the satellite data requirements from the World Climate Research Program. 
The third entry is liquid water. For measuring liquid water, microwave imaging radiometers are 
needed such as on SSM/I; Eos, which is the MIMR; and TRMM microwave imager TMI. For 
the liquid water profile, only a rain radar can make that measurement, and of course, the only 
precipitation radar available will be on TRMM. Therefore, since TRMM is only a three year 
mission, we need a follow-on in order to fully extend the coverage of GEWEX. 


Measurements 

Instrument type 

Candidate sensor 

EM339HHHMHN 




■mhhi 

Wind profile 

usEEMmmmm 

LAWS 

« i ? ii i ! j* immmm 

Temperature & humidity 
profile 

ER & microwave sounder 

KM irt ll " 1 


Liquid water (total) 

Microwave imaging 
radiometer 

SSM/I, MIMR, TRMM 
microwave imager (TMI) 

QSHBH 

Liquid water (profile) 

Rain radar 

Precipitation radar (PR) 


Cloud amount and cloud top 
temperature 

Imaging radiometer 
(vis & RI) 

AVHRR (3), MODIS 

ESHHI 

Cloud particle properties 

Polarimeter 

POLDER, EOSP 


Cloud particle profiles 

Millimeter wave radar 
altimeter 

(undetermined) 

(undermined) 

Atmospheric chemistry: minor 
constituents including ozone 

Spectrometer 

TOMS, GOME, GOMOS, 
SCIAMACHY, MOPITT, 
IMG, ILAS, HIRDLS, TES 

NOAA, ERS-2, ADEOS, 
POEM, EOS 

Aerosol: stratospheric 

Limb-scanning 

spectrometer 

SAGE, SCIAMACHY, 
ILAS 

EOS, POEM, ADEOS 

Aerosol: tropospheric 

Muld-directional imager 

MISR 

EOS 

Solar radiation 

Solar radiometer 

ACRIM 

EOS 

Earth radiation (top of the 
atmosphere) 

Wide -band scanning 
radiometer 

ScaRaB, CERES 

««■ 

Earth surface/ cloud multi- 
directional reflectance 

Multi-directional imager 

ATSR, MISR 

ERS, EOS, POEM 

Earth surface reflectance and 
colour 


avhrr <3) avnir, moths 

NOAA, ADEOS, EOS, 
POEM 


Figure 1. World Climate Research Program (WCRP) Satellite Data Requirements (From a 
presentation by P. Morel, Director, WCRP, Hamilton, Bermuda, January 1991). 
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Ocean surface colour 

Imaging radiometer (vis) 

SeaWIFS, OCXS, MODIS 

SeaWIFS, ADEOS, EOS, 
POEM 

Sea surface temperature 

Imaging radiometer or 
sounder (IR) 

avhrr (3> atsr, octs, 
MODIS, AIRS. IASI 

NOAA, ERS, ADEOS, EOJ 
POEM 

Ocean surface wind vector 

Microwave scatterometer 

AMI, NSCATT, 
STIKSCAT 

ERS, ADEOS, EOS, POEK 


Radar altimeter 

RA 

ERS, EOS, POEM 

Ocean surface topography 

Precision radar altimeter and 
orbitometry 

Single or 2-frequency 
altimeter, GPS + water 
vapour radiometer or 
DORIS 

TOPEX-POSEIDON, 

alt (4) 

Sea ice cover 

Microwave imaging 
radiometer 

SSM/I, MIMR 

DMSP, EOS, POEM 

Sea ice texture, edge and 
motion 

Imaging radar 

AMI, SAR 

ERS, RadarSat, JERS-1 

Snow cover 


AVHRR, MODIS, SSM/I, 
MIMR 

NOAA, EOS, POEM, 
DMSP 

Snow water equivalent and soil 
moisture 

(undetermined) 

(undetermined) 

(undetermined) 

Glacial and iceberg discharge 

High-resolution imager (vis) 
and radar altimeter 

TM, HRV, ASTER AND 
RA 

LANDSAT, SPOT, EOS 
and ERS, EOS, POEM 


Figure 1 (Continued) 


Footnotes: (1) The accuracy and vertical resolution of HIRS in the troposphere are insufficient 
for climate studies 


(2) GEWEX dedicated non-sunsynchronous earth observing mission on a 55*60 
degree inclined orbit (undetermined). 

(3) AVHRR calibration only marginally adequate for cloud classification, 
vegetation mapping and sea surface temperature measurements. 

(4) TOPEX follow-on high precision ocean altimetry mission (undetermined) 


GLOBAL HYDROLOGICAL CYCLE 


PROBLEM IS TO COMPLEMENT 
THE SURFACE-BASED NETWORK FOR 
DETERMINING GLOBAL PRECIPITATION 
AND SOIL WETNESS 


• GLOBAL PRECIPITATION ESTIMATES EXIST, BASED ON IR CLOUD 
IMAGES OR PASSIVE p-WAVE RADIOMETRY AND EMPIRICAL 
CLIMATIC CORRELATIONS 

• VERTICAL DISTRIBUTION OF RAIN DROPS NEEDED FOR DIRECT 
QUANTITATIVE ESTIMATION 

• RESOLVING THE DIURNAL CYCLE IS ESSENTIAL 

• TRMM IS ESSENTIAL FIRST STEP. TECHNIQUE TO BE FURTHER 
DEVELOPED IN FOLLOW-ON "EARTH-SCANNING" MISSION 


Figure 2. Global Hydrological Cycle Problem 
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Figure 2 presents the global hydrological cycle problem. We see that the problem is to 
compliment the surface based network for determining global precipitation and soil moisture. 
Global precipitation measurements estimates exist based on infrared cloud images or passive 
microwave radiometry and empirical climatic correlations. The vertical distribution of raindrops 
is needed for the direct quantitative estimation of precipitation. Resolving the diurnal cycle is 
essential which TRMM certainly can do. We need a non-sun-synchronous orbit in order to 
measure the diurnal cycle. TRMM is the essential first step, but the technique can be further 
developed in follow-on Earth scanning missions. 

TRMM FOLLOW-ON PLANNING 

Planning for tropical rainfall measurements to follow TRMM has been a part of the French Space 
Agency (CNES) planning since Bilan Enegetique du Systeme Tropical (BEST), which means the 
tropical energy budget, was first proposed in 1989 as a candidate to meet GEWEX space 
observing system requirements and for an operational mission to begin at the turn of the century. 
The proposers invited international cooperation, but the mission has never really been approved. 
In Hamilton, Bermuda where the GEWEX Joint Scientific Steering Group met in January of 
1991, the scientific need for continued precipitation measurements from space after the TRMM 
mission was presented by Professor Pierre Morel, Director of the World Climate Research 
Program. About that same time the Japanese discussed a plan to begin a TRMM follow-on 
mission. At the fourth session of Joint Scientific Steering Group for GEWEX, which was held 
January 1992 in Tokyo, Japan, a resolution was agreed upon that said, "TRMM follow-on 
missions, that could continue the cooperation between USA and Japan, include discussions of a 
satellite launched after the year 2000 in a 55 degree inclined orbit with an advanced multi-beam, 
multi-frequency precipitation radar, and an advanced passive microwave imager." (From 
presentation by Thomas La Vigna, GSFC, January 1993). 

A TRMM follow-on was also discussed at the International Workshop on the Processing and 
Utilization of the Rainfall Data Measured from Space which was held March 1992 in Tokyo, 
Japan. At that time it had been established that precipitation is an important component of the 
Earth system which influences global weather and the climate. Moreover, there now exists 
strong evidence that rain rates can be measured from space platforms to within the scientific 
requirements for improving weather forecasting and for improving the models used in the 
prediction of climate change. Continuation of the measurements is crucial for detecting and 
monitoring global change of precipitation. In view of these considerations, the workshop 
participants strongly endorsed the concept of a follow-on mission to TRMM. They urged that 
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the necessary steps be taken by NASA, NASDA, ESA, and CNES to begin pre-phase A studies 
immediately for such a cooperative mission. A candidate configuration for a 55 degree 
inclination orbit, which allows sampling over the diurnal cycle and full coverage of the mid- 
latitude weather, was put forward. The payload was suggested to be similar to that of TRMM 
except for upgrades in the instruments, such as the dual frequency radar, multi-channel 
microwave radiometer, and visible and infrared radiometers. 

At the International Symposium on Active Sensors and Non-Synchronous Missions Dedicated to 
GEWEX, held June 1992 in Jouy en Josas, France, strong science recommendations came for the 
continued measurement of precipitation from space after TRMM. The meeting proposed a non- 
sun-synchronous mission to extend the results of TRMM beyond the tropics in order to realize 
the global objectives of GEWEX. The various space agencies were encouraged to define the 
scope and specification of a TRMM follow-on mission and sensors for appropriate interactions 
with the scientific community, to coordinate and conduct the corresponding feasibility studies 
and development, with a view to creating the conditions for a program decision, and accomplish 
studies and trade-off assessments of several options for the rain radar, for example, dual beam, 
two-frequency, and dual polarization. 

At the Asia-Pacific International Space Year Conference, "The Earth in Space", held November 
1992 in Tokyo, Japan, the opening key lecture by Professor P. Morel discussed a follow-on 
TRMM for continued measurements from space of rainfall and importance to climate models. 
The opening day lecture by S. Tilford entitled "U. S. Earth Observation Satellite Program" 
dis cussed future potential space missions including a TRMM follow-on, TRMM-2. 

In the GEWEX Dedicated Space Mission Study Workshop, held November 1992 in Tokyo, 
Japan, the focus was on a TRMM follow-on mission. The workshop recommendations included 
unique and essential scientific requirements for GEWEX and global climate science as being the 
observations of radiation, clouds, and rain which show diurnal variations, and for the latter two 
parameters, vertical resolutions are indispensable for characterizing relationship among these 
quantities. In order to meet these requirements, a TRMM follow-on mission is required to 
provide data continuity over the tropics and to extend coverage to higher latitudes. In view of 
these considerations NASDA's proposal for a TRMM follow-on mission aiming at a launch in 
the year 2000 was endorsed as a basis for development In addition to the original TRMM 
instruments, inclusion of a cloud profiling radar was recommended to be considered for 
monitoring the global 3-dimensional distributions of clouds. 


295 



At the 6th Standing Senior Liaison Group (SSLG) Meeting, held December 1992 in Washington, 
D.C., the NASDA presentation entitled "Study on TRMM Follow-On" also recommended that 
the TRMM follow-on mission to meet GEWEX requirements be started. There are a number of 
improvements which could be made to the basic TRMM in the TRMM follow-on. However, in 
order to meet the schedule laid out for continuing the observations, which would mean a launch 
in the year 2000, it appears that changes for the TRMM follow-on will be minimal. They will 
include only the increase in the inclination angle to 50 to 60 degrees of latitude, in order to cover 
the globe more geographically, and an increase in the altitude of the orbit so that the space craft 
will be able to have a five year lifetime. This basically means that the radar might have to be 
increased in power and antenna size. The passive microwave imager might also have to be 
increased in antenna size in order to keep the resolution of the footprint small enough to observe 
rainfall accurately, especially in convective situations. 

Other improvements, such as the inclusion of a radar with multiple frequencies, or a push-broom 
type scan pattern, or polarization, or all of these, will have to wait until a successor mission to the 
TRMM follow-on. Likewise, a cloud radar whose feasibility has not been demonstrated yet will 
also have to await development at the laboratory and aircraft levels to prove the concept Then, 
possibly, it can be implemented on a mission sometime after the year 2000. 

SUMMARY REMARKS 

If we are to maintain data continuity, we need to initiate a TRMM follow-on mission very soon. 
NASA is considering this action and may fund a study on the subject during this fiscal year. The 
resources to fund the phase C/D portions of the mission remain to be identified. They will have 
to compete with all the other ideas that have been proposed for specialized Earth Probe Missions 
in the years following the turn of the century. The measurement of rainfall from space is and 
ought to be a high priority in studying the Earth system and global change. It remains to be seen 
whether we can improve upon the very fundamental measurements that TRMM will make. How 
we go about this should be the subject of continuing research in the years between now and the 
end of the decade. 

Acknowledgments: The author wishes to acknowledge the able assistance of Dr. William W. Vaughan in 

reviewing and editing the manuscript and Ms. Jean Taylor for typing the manuscript 
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CALIBRATION 

PRACTICAL EXPERIENCE WITH ERS-1 


• Introduction 

• Radiometric Calibration 

• Geometric Calibration 

• Phase Calibration 

• Polarimetric Calibration 
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Introduction 


Basic SAR measurement parameters are : 


• Radar backscattering 

• Target position 

• Target speed 

• Polarisation 


SAR calibration facilitates quantative 
measurements needed to derive geophysi- 
cal parameters of the area under observa- 
tion from basic SAR measurements (e.g. 
soil moisture, biomass, ocean wave en- 
ergy, ocean currents, ice type, ice flow, ) 
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Radiometric Calibration 


Radiometric calibration is relating the SAR 
output in terms of digital number to o. 

• Design for stability 

• Sensor calibration 

— pre-flight characterisation 

— internal calibration 

— external calibration 

• Data calibration 
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Characterisation parameters for ERS-1 

• Transmitter Parameters 
— Frequency 

— Pulse power 
— Pulse duration 
— Chirp bandwidth 

• Antenna Parameters 

— Relative azimuth pattern 
— Absolute elevation pattern 

• Receiver Parameters 
— Sensitivity 

— Bandwidth 
— Power transfer function 

• Calibration Parameters 

— Internal calibration stability 


303 


A Instantaneous Amplitude 


Corresponding to 
Peak RF Power 


Time 
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Internal Calibration 


The objective of internal calibration is moni- 
toring system parameter changes with 
time, such as 

• Output power 

• Receiver gain 

• Linear distortion (amplitude & phase) 

• Non-linear distortion ( amplitude & 
phase) 

• Noise level 

• Level of spurious signals 
Internal calibration methods include 

• Calibration using transmit pulse sample 
(ERS-1 method) 

• Reference signals : noise, cw 

• Transmit power measurements 
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External Calibration 


The objective of external calibration is to 
measure and monitor system components 
outside the internal calibration loop. 


External calibration uses made-made or 
natural targets. 


External calibration can be used to deter- 
mine the radar image calibration factor 
either in combination or without reference 
to internal calibration. 


AMI Transponder Specifications 



Wind Mode 

Image/ Wave 

Radar Cross Section 

87. bdBm 2 

65.0 dBm 2 

Adjustment Range 

-j-0, — bdB 

+0, —bdB 

Calibration Accuracy 

±.bdB 

±.5dB 

Cross- calibration 
Accuracy 

±.2dB 

i.2dB 

Stability (Over 3 Years) 

±.ldB 

±.ldB 
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Gain (dBi) 


ERS-1 Transponder Antennas 
Gain Values for 12 Antennas 



' Wgg mte rs // serial rotors 565 to 576 

Average Gain = 27.81 dB 
Standard Deviation = 0.06 dB 
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Degrees 


Degrees 


E and H plane patterns, co and cro$s-polarisatk>n. 
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2432A OSCILLOSCOPE 
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FIGURE 2. RADAR CROSS SECTION CACIPRATI 







WELATIV* HO* «CL. MULTIPATH AMPLITUO* 


miwm ESTIMATE 
RF#2, Tx571, Rx572 



RELATIVE RANGE (a* 


ERS1; SAR AR0T2 
MUUFA7H ESTIMATE 
RF«,Tx571,Ri572 
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Transponder Calibration Accuracy 




SAR 

SCATT 




full 

range 

full 

range 


A 

Antenna Gain Error 

0.028 


dB 


Pointing Error 

0.53 

« 

deg 

B 

Electronic Gain Error 

10131 

0.056 

dB 

C 

Antenna Coupling Error 


- 

dB 

■ 

Transponder Stability 
A+B+C 

.08 

0.088 

dB 

TT 


1 a 

1 a 


Transponder Stability 
1 <j — (A+B-f C)/3 

0.027 

0.029 

dB 

E 

Clutter & Noise 
against <7° = 0 background 

0.0495 

0.0217 

dB 

F 

Coherent Multipath 

0.0207 

0.0388 

dB 

G 

Incoherent Multipath 

B 

EZHBI 

dB 


Single Transponder Observation Error 
V~D 2 + E 2 + F 2 + G 2 


0.093 

dB 

« 

Atmospheric Loss Uncertainty 



dB 

■ 

Combined Error = 

VD 2 + E 2 + F 2 + G 2 + W 2 



dB 

111 

Transponder Calibration Accuracy 

0.2 

0.2 

dB 
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Representative SAR Calibration Budgets 


SAR Radiometric Accuracy 





1 a 


B 

Absolute Calibration Error 

a 

0.14 

dB 

B 

Cross Swath Calibration Error 

b 

0.08 

dB 

C 

Radiometric Stability Drift 

b 

0.25 

dB 


Nominal Calibration Accuracy 


0.3 

dB 

D 

Radiometric Stability Drift 
Nominal - Observation 

c 

0.25 

dB 

E 

Across Swath Characterisation Error 

c 

0.1 

dB 

F 

Across Dynamic Range Characterisation 

c 

0.1 

dB 

G 

Atmospheric Loss Uncertainty 

c 

0.07 

dB 


Total Radiometric Error = 

VA 2 + W + C 2 + \[D 2 + E 2 + F 2 + G 2 


0.6 

dB 


(a) Dominate^. by traaspoeder calibration 

(b) Single sample error ^reduced by in-fligbt monitoring 

(c) Random error per observation 


Absolute Calibration Error 


1 


Number of 
Samples 

Single Sample 
Error(l a) 


H 

Transponder Calibration Accuracy 

3 

0.2 

dB 

I 

Antenna Gain Characterisation 

3 

0.1 

dB 

J 

Single Transponder Observation Error 

3 

0.07 

dB 

K 

Atmospheric Loss Uncertainty 

3 

0.07 

dB 

1 

Combined Error = 
jH 2 /3 + P/3 + J 2 /3 + K 2 / 3 

1 

] 

0.14 

dB 
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SAR- SWATH 
ESA/ESTEC, ERS 


VENICE ORBIT 

-1, 3-day repeat orbit ( 43), 24.365 deg SAR 


Figure 7: Transponder Sites Flevoland, Commissioning Phase Coverage 
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SAR CALIBRATION HISTORY 

COMMISSIONING PHASE & ICE PHASE (from 7-SEP-1991 to 11 -MAR- 1992) 
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Calibration of other areas 


The application of the calibration factor to 
other areas requires knowledge of geome- 
try (range & antenna angle) and stability 
over time. 


High resolution information about the an- 
tenna diagram can be obtained from im- 
agery of extended uniform targets such as 
the tropical forest or ice shelves. 


No absolute information is required for 
these ’targets of opportunity’. Stability over 
time and a smooth dependence of g° on 
incidence angle is necessary. 









5JJEC-1992 14:42:18597 









Data calibration 


Data calibration include the following steps 

• Correct raw data for a number of system 
parameters 

— Linear distortion (amplitude & phase) 

— Non-linear distortion 

— Detector imperfections (DC-biases, 
amplitude & phase imbalance) 

— time variations as measured by inter- 
nal calibration 

— noise bias 

• Application of calibration factor to the 
entire image, taking into account ge- 
ometry and antenna angle. 
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Geometric calibration 


Geometric calibration can be achieved fol- 
lowing the same principles as for radiomet- 
ric calibration. 

• Design for stability 

• Characterise SAR system in terms of 
delay, phase error and frequency offset 

• Calibrate using point target with accu- 
rately known position 

Geometric image calibration is carried out 
in the data calibration stage by correction 
for knpwn system biases and applying ge- 
ometry information (platform position & atti- 
tude and terrain models). 
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Phase calibration 


During the formation of an intensity image 
the phase information in the radar data is 
lost. For applications in interferometry and 
polarimetry complex images are formed 
which provide two data points per pixel, 
amplitude & phase or real & imaginary 
parts of the complex signal. 


Careful and elaborate phase characterisa- 
tion of both the SAR sensor and the SAR 
processor are required. Time stability is 
critical for the above applications. 



Polarimetric calibration 


Polarimetric calibration can be treated es- 
sentially as a multi-channel extension of 
the methods discussed above. Polarimetric 
SAR differ from single channel SAR in the 
following ways. 

• There are two orthogonal transmit chan- - 
nels to consider and four receive chan- 
nels. 

• Radiometric calibration can be sepa- 
rated into absolute calibration and rela- 
tive between channel calibration 

• Geometric calibration can be separated 
into absolute calibration and relative be- 
tween channel calibration 

• Phase calibration between channels is 
essential 

• Crosstalk has to be taken into account 
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SYNTHETIC APERTURE RADAR SIGNAL 
PROCESSING: TRENDS AND TECHNOLOGIES 
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OVERVIEW OF SAR PROCESSING 

SAR DATA SYSTEM OVERVIEW 



DAT A UT I LI ZAT ION AND I NFORMATION 
EXTRACT ION SYSTEM 



SAR CORRELATOR: TRENDS AND TECHNOLOGIES 
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- CURRENT, FUTURE TRENDS 


SAR DATA SYSTEM OVERVIEW 
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SAR CORRELATOR: TRENDS AND TECHNOLOGY 

SAR PROCESSOR ALGORITHM OVERVIEW 
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APPROACH 2: APPLY TWO 1-D FILTERS FORMING A RANGE COMPRESSED IMAGE 




SAR CORRELATOR: TRENDS AND TECHNOLOGIES 

ALGORITHM OVERVIEW 
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APPLICATION IN INTERFEROMETRY AND POLARIMETRY 

- WAVE EQUATION: HIGH SQUINT (75 DEGREES), PHASE PRESERVING 

- CHIRP SCALING: REQUIRES LINEAR FM CHIRP, NO INTERPOLATION 


SPECAN ALGORITHM FLOWCHART 


CHIRP 

PARAMETERS 

I 

RANGE REF 
FUNC GEN 


DOPPLER 

PARAMETERS 


AZIMUTH REF 
FUNC GEN 


SIGNAL DATA 



MULTI-LOOK IMAGE 
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RANGE/DOPPLER ALGORITHM FLOWCHART 

WITH SECONDARY RANGE COMPRESSION (SRC) 



AZIMUTH REF 

■1 

AZIMUTH REF FUNC 

FUNC GEN 


CMPLX MULT 


INVERSE 
AZIMUTH FFT 


COMPLEX IMAGE 
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SAR CORRELATOR: TRENDS AND TECHNOLOGIES 

RANGE DOPPLER ALGORITHM PERFORMANCE 



( it ) sjojj9 esBL|d o;iejpenb 
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RESIDUAL PHASE ERRORS FOR THE SECONDARY RANGE COMPRESSION (SRC) 
AND THE SQUINT IMAGING MODE (SIM) ALGORITHMS (CHANG, ET AL, 1992) 




2-D TRANSFORM (WAVE EQUATION) 
ALGORITHM FLOWCHART 



COMPLEX IMAGE 


ROGCA, 1987 ; CAFFORIO, 1991 


346 




2-D TRANSFORM (CHIRP SCALING) 
ALGORITHM FLOWCHART 


SIGNAL DATA 



COMPLEX IMAGE 


RANEY, 1991 
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SAR CORRELATOR: TRENDS AND TECHNOLOGIES 

ARCHITECTURE OVERVIEW 
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- UNLIMITED COMPUTATIONAL CAPACITY, GRACEFUL DEGRADATION, LOW 
COST PER PROCESSOR, DIFFICULT I/O, COMPLEX SOFTWARE DEVELOPMENT, 
FUTURE FOR ALL LARGE COMPUTATION TASKS 

- MIMD (DUAL PARALLELISM): THINKING MACHINES CM-5, INTEL l-WARP 

- SIMD MACHINES: DEC MASSPAR 


SAR CORRELATOR ARCHITECTURES 
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DISADVANTAGES 

- HIGH COST PER MFLOP, PERFORMANCE LIMITED BY CENTRAL BUS DATA RATE 




ADVANCED SAR CORRELATOR ARCHITECTURE 

BRITISH AEROSPACE ERS-1 PROC 
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SAR CORRELATOR: TRENDS AND TECHNOLOGY 

SAR CORRELATOR ARCHITECTURES 


CO 

GC 

O 

CO 

CO 

LU 

o 

o 

GC 

Q. 

LU 


LU 

CL 

Q. 


O 

H 

CO 

3 

O 





CL 

o 

_l 

u_ 

s 

ft 

o 

o 

£ 

o 


CL 

X 

o 

3 

o 

CE 


UL 

O 

CO 

s 

x 

Ui 

I- 


oc 

< 

CO 


CO 

Ui 

o 

< 

H 

z 

< 

> 

D 

< 


oc 

O 


U_ 


LU 

CE 

3 




LU 

I- 


X 

o 

c 

< 


a 

o 


CO 

UJ 

X 


3 


a 


UI 

X 


> 


< 

o 

a 


z" 

o 


< 

Q 

< 

X 

o 


Ui 

o 



uT 

X 


H 

O 

UI 

I- 


CO 

UI 


a 

< 

H 


< 

> 

a 

< 

co 

5 


o 

x 

< 


ui 

_i 

CD 

X 

UJ 



z 


■ 


• • 


351 


ALGORITHM APPROXIMATIONS 




SAR CORRELATOR: TRENDS AND TECHNOLOGY 

SAR CORRELATOR ARCHITECTURES 
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DISADVANTAGES - PRIMITIVE OPERATING SYSTEMS, REQUIRES SPECIAL ALGORITHMS FOR 

EFFICIENCY, POOR I/O TO EACH PROCESSOR 


SAR CORRELATOR DESIGN 

EMMA MIMD CONCURRENT 
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SAR CORRELATOR: TRENDS AND TECHNOLOGIES 

OPERATIONAL SAR PROCESSING SYSTEMS 
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- ERS-1 PAF AND FAST DELIVERY PROCESSORS 

- NASA/JPL SEASAT, SIR-B AND SIR-C PRECISION PROC 

- DLR XSAR, ASF ERS-1 AND JERS-1, NASDA JERS-1 
2-D TRANSFORM 

- RESEARCH PROCESSORS: ESA, CCRS, JPL, DLR 
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ASSUME 10 SHIFT PER WEEK OPERATIONS WITH 50% REPROCESSING 
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SAR CORRELATOR: TRENDS AND TECHNOLOGIES 

SUMMARY 
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LEVEL PRODUCTS 

- USER INTERACTION IN MISSION OPERATIONS, TELESCIENCE 



Recent Advances and Plans in Processing and Geocoding 
of SAR Data at the DFD 

W. Noack 

German Aerospace Research Establishment (DLR) 
German Remote Sensing Data Center, DFD 
W— 8031 Oberpfaffenhofen, Germany 





1. Introduction 

Because of the needs of the future projects like ENVTSAT and the experiences made with the 
current operational ERS— 1 facilities a radical change in the SAR processing scenarios can be 
predicted for the next years. At the German PAF several new developments have been initial- 
ized which are driven mainly either by user needs or by system and operational constraints 
(’’lessons learned”). 

At the end there will be a major simplification and uniformation of all used computer systems. 
Especially the following changes are likely to be implemented at the German PAF: 

transcription before archiving, processing of all standard products with high throughput direct- 
ly at the receiving stations, processing of special ’’high -valued” products at the PAF, usage of 
a single type of processor hardware, implementation of a large and fast on-line data archive 
and improved and unified fast data network between the processing and archiving facilities. 

In the following a short description of the current operational SAR facilities as well as the fu- 
ture implementations will be given. 


2. Archiving and Catalogue System 


At present, archiving of SAR products from ERS-1 and later on also from X-SAR is done 
using Archival Optical Disks (WORMS) with a capacity of 2 GByte per disk. The disk drives 
are installed at each processing system and the disks are carried by hand to the facility where 
they are needed. The main disadvantages are the slow throughput, the limited storage capacity, 
the susceptibility to operator errors and the high maintenance requirements of the drives. 

An improvement of this situation is represented by the Epoch- 1 system, which consists of a 
mixture of non-erasable and erasable disks with a total on-line capacity of 400 GByte. The 
system is used for many projects, among them are the storage of AVHRR data and the back- 
ground storage for the image processing laboratories. But there are still the disadvantages of 
relatively slow throughput and the limited storage capacity. Aso, this system would not be able 
to guarantee the on-line access of one year ERS-1 PAF production of SAR and Geocoding 
which is actually more than 400 GByte. 
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Therefore, in view of ERS — 2 and ENVISAT the next generation archiving system is being stu- 
died. During a first analysis of the market two systems have been selected for a detailed inves- 
tigation, of which the one system allows a total on-line capacity of up to 10.000 TByte. The 
study is carried out in 1993 with the goal of having the first unit installed and operational end 
1994. 

Due to historic reasons there are a variety of catalogues for the different kind of sensors. The 
biggest catalogue is the one for ERS— 1. In near future these catalogues will be accessible by 
a recently developed common user interface. This interface is also considered to be available 
for external SIR— C /X— SAR users and shall be implemented before the second SIR— C / X— 
SAR mission, depending on availability of funding. 

3. Processing 

With respect to processing two major milestones have been achieved within the last year: first 
the Multisensor SAR Processor MSAR has become operational for ERS— 1 beginning ’92 and 
second the precision processor of the X— SAR PAF has been accepted in December ’92. The 
following two figures are showing the hardware and the software block diagram of the XPS sys- 
tem. A SUN UNIX computer enhanced by an array processor of STAR Technologies with di- 
rect access to the solid state Corner Turn Memory are the major hardware units of the proces- 
sor. The X— SAR data will be read in from a SONY recorder using a real time capable Frame 
Synchronizer and Formatter. 
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figure 1: XPS hardware block diagram 
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The software structure of the XPS is shown in figure 2 above. The Control System allows among 
many other options especially automatic and parallel processing of three different orders which 
are divided in the steps Transcription, Processing, and Archiving. The Control System has been 
accepted already and is since then in operational use. The Correlator has been accepted also 
and is in the integration phase. It uses the slightly modified and supplemented Range— Doppler 
algorithm. The other subsystems are either in the final acceptance phase (Formatter, Archiver) 
or in the design and development phase (Screener, Transcriptor). 

In opposition to ERS- 1, where the MSAR Screener is mainly used for archiving and quality 
estimation of raw data, the Screener of the XPS system is designed as a preprocessor of the Cor- 
relator and covers many different tasks. Among these are the Doppler Screening, the Doppler 
Ambiguity Resolving, the calculation of the internal calibration parameters, the extraction of 
header data, the inventarisation of the tape contents, the analysis of the mission data base and 
the geolocation of the raw data. 

Control System, Screener, Processor, Formatter, and Archiver have been written in Ada except 
the modules running on the array processor. This will allow for a reuse of components in a dif- 
ferent environment in future projects. 

Beginning this year a new project has been started with the goal of the development of a high 
precision SAR processor making use of a powerful workstation only without external devices 
like an array processor or a Corner Turn Memory. The system will be quite flexible with ’’open” 
interfaces between the modules allowing the easy exchange of components. It shall mainly sup- 
port the development of new, higher -valued products and shall be used as a prototype proces- 
sor for future sensors like ENVTSAT. The programming language will again be Ada and existing 
software shall be reused. Especially the Control System of the XPS shall be used for the para- 
metrisation and monitoring of processes as well as other components. Because of the foresee- 
able increase in compute power of workstations throughput requirements will be investigated 
during the first phase with low priority only. This project is in a preparatory phase. The develop- 
ment will be done inhouse with possible cooperation with universities or institutes. 

4. Geocoding 

The present Geocoding System GEOS which is in operational use for ERS— 1 and which will 
be adapted to X— SAR runs on SUN workstations (SPARC2) and will be implemented on the 
new generation (SPARC10) in early 1993. The system generates mainly ellipsoid or terrain cor- 
rected products with a throughput of several products per day. Currently no batch processing 
facility is implemented and therefore, each product has to be generated individually. 

In future it is foreseen to include the generation of the ellipsoid corrected product directly in 
the SAR processor and additionally direct at the receiving stations. At present the throughput 
of the terrain geocoded products is one item per day. The limiting factors are the time consum- 
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ing determination of tiepoints, the quality control and the compute power. To improve the situ- 
ation the tiepoint interface will be upgraded and the quality control will be automised. 

Additionally, the list of geocoded products shall be extended by a radar mosaic, a geocoded 
map quadrant and coregistered data sets. Each extension requires software upgrades and will 
be performed in part with the support of industry and universities. 
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Tend to be restricted to expert users 


JPL . Third Spaceborne Imaging Radar Symposium 

Multimission SAR Processing System Architecture 
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HDDR-HIgh Density Digital Recorder 
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• Facilitates Intersensor comparison, global model comparison 
Preferred format for the Earth Observing System Data and Information System 
Version 0 (and probably for versions in the late 1990's) 

HDF is supported by a large software tool base for science data visualization 
Translator developed by the Alaska SAR Facility for converting CEOS to HDF 
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DOMSAT with image data compression (1 Mbps data channel) 
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* EOSDIS Handbook, 1992 
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*EOSDIS Handbook, 1992 
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Complex = FACILITIES FOR LEVEL 0 AND 1 DATA HANDLING 
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EOSDIS - Earth Observing System Data and Information System 

DAAC - Distributed Active Archive Center 

IMS - Information Management System 

PGS - Product Generation System 

DADS - Data Archive and Distribution System 
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CONCLUDING REMARKS 
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SAR T/R Module Architecture 



T/R Module is a Unique Assembly of the SAR Functional 

RF Components 



Simplified Block Diagram of SAR System 
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T/R Module is a Key Element in the Overall SAR system 



Requirements Traceability 
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Requirements Traceability (cont’d) 
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- Antenna Size and Weight 

- Deployment and Stowage 



Typical SAR Antenna Configuration 
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T/R Module is an Integral Part of Antenna Design 

















Typical SAR Antenna Panel Cross Section 
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T/R Module Characteristics Are Key Drivers of Electrical, Mechanical 
and Thermal Designs of Antenna and SAR System 




Desired Module Characteristics 
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The Better the Module, the Better the Antenna and Overall 

SAR System 




Typical T/R Module Specifications/Requirements 
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T/R Module Specifications are a Combination of Derived, 
Self-Imposed, and Direct Flow Down Requirements 



L- Band T/R Module Development 
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L-Band T/R Module is Ready For Insertion into SAR Applications 

With Little or no Modifications 




SEASAT L-Band 500 W SSPA 
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SIR-C L-Band T/R Module 
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SBR L- Band T/R Module 
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GBR C-Band T/R Module 
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SIR-C C-Band T/R Module 




COBRA C-Band T/R Module 
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Airborne X-Band T/R Module 
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X-Band T/R Module 




Concluding Remarks 
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Advanced Antennas for SAR Spacecraft 


William B. Gail 

Ball Communication Systems Division 
P. O. Box 1235 
Broomfield, CO 80038-1235 
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Abstract 

Single and multi-frequency antenna concepts have been 
developed to evaluate the feasibility of building large 
aperture polarimetric SAR systems to be launched in 
low cost launch vehicles such as the Delta II. The 
antennas are 18.9 m long by 2.6 m wide (L-band) and 
achieve single polarization imaging to an incidence 
angle of 55° and dual/quad imaging to 42°. When 
combined with strawman spacecraft designs, both 
concepts meet the mass and volume constraints 
imposed by a Delta II launch 

Introduction 

Advanced spacebome SAR missions generally require 
large aperture antennas to meet the needs of the user 
community. SAR antenna concepts presented to date* 
however, have not fully realized the potential of 
available antenna technologies. The two EOS SAR 
industry briefing studies, for example, indicated that 
the largest L/C/X SAR that could be launched in a 
Delta II has an antenna aperture of only 2.6\10.9 m. 
The performance of this antenna, particularly with 
regard to maximum incidence angle ranges, is 
considerably less than that desired by the science 
community. 

To evaluate the limiting performance of a mass, 
volume, and cost constrained SAR. an antenna 
feasibility study was performed. Four constraints were 
placed on this study: 

1) Satisfy known science requirements so as to place 
the user community' ahead of the technology 

2) Base the design on the proven S1R-C heritage, 
minimizing design changes but incorporating 
new technologies where most useful 

3) Design to a Della 7920 launch to limit launch 
cost 

4) Generate a design which minimizes the number 
of parts and processes to reduce cost 

Alternate design approaches were then examined, 
based on these constraints, to identify feasible (not 
necessarily optimal) designs. 


Requirements 

The science requirements specified by the user 
community for EOS SAR were used to drive the design. 
Mission/system requirements were derived from the 
science requirements. An orbit at 700 km was selected 
to minimize atomic oxygen exposure and drag and thus 
extend the spacecraft lifetime. Requirements imposed 
by selection of the Delta 7920 as a launch vehicle are 
strong drivers for the spacecraft and antenna designs. 
While other launch options could be more effective, 
designing to the Delta option effectively identifies the 
key technology and cost drivers for the antenna. Both 
single and multiple frequency design options were 
considered. 

L-band Advanced SAR (LASAR) 

The single frequency' quad-polarization L-band 
Advanced SAR (LASAR) concept developed during 
this study is shown in Figure 1. Characteristics of this 
antenna are summarized in Table i. The LASAR 
antenna is 18.9 m long and consists of seven leaves 
that fold for launch. Bus and payload subsystems are 
contained on the backside of the square vertical 
structure in this figure. 



Figure 1 Conceptual illustration of the LASAR 
spacecraft showing the deployed 18.9 m antenna. 
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Table 1 Characteristics of the LASAR antenna. 


CHARACTERISTIC 

VALUE 

Array size (m) 

18.9x26 

Antenna mass (kg) 

511 

Peak radiated power (W) 

4900 

Number of T/R modules 

224 

Peak output power per T/R circuit (W) 

11 

Number of leaves 

7 

Leaf size (m) 

2.68x2.6 

Leaf thickness (cm) 

6.4 

Panels per leaf 

2 

Polarization 

quad 

Elements per stick 

6 

Sticks per panel 

16 


The aperture uses the same six-element-per-stick 
microstrip patch architecture used on SIR-C, as 
illustrated in Figure 2. This approach maximizes SIR- 
C heritage in manufacturability as well as performance. 
The peak radiated power of 4900 W requires a module 
peak power of only 28 W, well within the capability of 
current devices. Because the antenna support structure 
is integral to the panels, total panel thickness is only 
6.4 cm. 

The use of thin panels allows the 18.9 meter antenna to 
be folded and stowed within the tight launch shroud of 
the Delta, as shown in Figure 3 Panels are connected 
by simple hinges and deployed with reliable motorized 
actuators. Figure 4 illustrates the first step of this 


deployment sequence. 

1 

2 A 

2 

i 

m 


i — ui — u — nm i-m cHZDHrrra-cztn 


nrrrTf>n-n 





i = >o-cir^ 

t~H h rvi-n-c^ i HTTcrnoo 

rT^-T-r-rv-i-CT-ra n-r—w— ro-cn-o 

rrr>r-rH-H-T-ri mz>T^T>CZH=l 


r-w-^CTv-w— i-ct 





-< 2.68 m- >■ 


Figure 2 Layout of a leaf showing the microstrip patch 
elements. 


The pattern characteristics of the LASAR antenna are 
very similar to those of the SIR-C antenna, the most 
obvious difference being that the azimuth beamwidth is 
considerably smaller. Figure 5 shows the H- 
polarization transmit pattern when the beam is scanned 
15° off boresight 



Figure 3 Stowed configuration of the LASAR 
spacecraft within a Delta 7920 shroud. 



Figure 4 Initial deployment sequence for the LASAR 
antenna. 
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THETA ELEVATION CUT DEGREES 


Figure 5 Sample antenna pattern for the LASAR 
antenna showing the H-polarization transmit beam 
scanned at 15°. 

Antenna flatness is maintained through the structural 
integrity of the panels, which are designed with 
sufficient stiffness to withstand the launch load using 
materials carefully selected to minimize bending when 
the antenna is subject to thermal gradients. A thermo- 
mechanical analysis showed that the antenna maintains 
the required /V20 flatness with a 20% margin. 
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Figure 6 LASAR performance versus incidence angle: 
a) ambiguity and beamwidth limited swath for single 
and dual/quad modes, b) thermal noise equivalent a 0 . 


SAR performance of the LASAR design meets or 
exceeds all of the known science requirements for EOS 
SAR. Figure 6 demonstrates both incidence angle and 
sensitivity performance. The 18.9 m aperture supports 
imaging to incidence angles of 55° for single 
polarization and 42° for dual/quad polarization. This 
satisfies the EOS SAR science requirements of 50° and 
40° respectively and is well beyond the 40°/30° 
performance of the 10.9 m aperture presented in the 
EOS SAR industry briefing. Thermal noise equivalent 
a 0 of this antenna is -41 dB at 45°. The maximum 
SCANSAR swath is 670 km for single polarization and 
400 km for dual polarization 

Using mass estimates developed during this study and 
the EOS SAR industry studies, a mass budget was 
derived for the LASAR design. This budget is shown 
in Table 2. Mass savings are obtained primarily from 
the antenna mass reduction and use of a single- 
frequency instrument. Additional mass reductions 
from the industry study values, particularly in the 
instrument, are expected with further analysis. The 
LASAR approach meets the Della launch requirement 
with considerable margin including design 
contingency. 


Table 2 Mass budget for the LASAR spacecraft 


ITEM 

MASS (kg) 

Antenna 

511.39 

Flight Electronics 

263.00 

RF electronics 

50.00 

DDHA 

60.00 

CTTA 

75.00 

PC DA 

20.00 

DDRE 

20.00 

Harness 

30.00 

Adapter plate 

8.00 

Bus System (dry) 

1268.00 

Structure 

298.00 

AACS 

102.00 

Cabling 

15.00 

Propulsion 

62.00 

Thermal control 

30.00 

Power 

358.00 

Command and data 

152.00 

Mechanisms 

148.00 

Communications 

103.00 

Miscellaneous 

698.48 

Launch adapter/misc 

100.00 

Propellant 

90.00 

Project reserve 

100.00 

Contingency 

408.48 

TOTAL (kg) 

2740.87 

MARGIN (kg) 

559.13 
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Muitiband Advanced SAR (MASAR) 


L-BAND LEAF 


The performance advantages obtained through the 
LASAR design were extended to a multi -frequency 
system in the MASAR design. The design employs the 
same general seven-leaf configuration used shown in 
Figure 1 to obtain an 18.9 meter L/C/X aperture. The 
characteristics of the MASAR antenna are summarized 
in Table 3. Peak transmitted power levels were selected 
to satisfy 7 the science requirements while maintaining a 
feasible spacecraft DC power budget. Module peak 
powers are 1 1 W at L-band, 4 W at C-band. and 4 W at 
X-band. Panel thickness is the same as for LASAR, 
providing for the same stowage and deployment 
procedures. 

Table 3 Characteristics of the MASAR antenna 


CHARACTERISTIC 


VALUE 


Array size (m) 
Antenna mass (kg) 
Number of leaves 
Leaf size (m) 

Leaf thickness (cm) 
Panels per leaf 


18.9x2.6 

711 

7 

2.68x2.6 

6.4 

2 


L 

C 

X 

Polarization 

quad 

quad 

dual 

Peak radiated power (W) 

1900 

1400 

1200 

Number of T/R modules 

224 

224 

224 

Peak power per T/R circuit (W) 

11 

4 

4 

Sticks per panel 

32 

64 

64 
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Figure 7 Layout of a shared aperture leaf showing L- 
band elements and regions of C- and X-band elements. 

Table 4 Mass budget for the MASAR spacecraft. 


To obtain the largest possible aperture size for each 
frequency within the 18.9 x 2.6 m area constraint, a 
shared aperture approach was used. As shown in 
Figure 7, the L-band array uses the full area to achieve 
an aperture size that is the same as that for LASAR. 
The C-band and X-band apertures share the space 
occupied by the L-band aperture. This approach is 
known to be feasible at the element level: further work 
is in progress to understand the impact on performance 
at the array level. 

Antenna flatness is maintained through 1) the 
structural integrity' of the panels, which are designed to 
be individually flat to the required X-band 1/20 flatness, 
and 2) use of active hinges between panels. A thermo- 
mechanical analysis showed that, with the inclusion of 
active hinges, the antenna maintains the X-band a/20 
flatness across the enure array. 

Table 4 provides a mass budget for the MASAR 
concept, showing that a MASAR spacecraft meets the 
Delta launch mass requirement, although with only a 
small margin. 


ITEM 

MASS (kg) 

Antenna 

711.08 

Flight Electronics 

383.00 

RF electronics 

91 .00 

DDHA 

75.00 

CTTA 

114.00 

PC DA 

27.00 

DDRE 

23.00 

Harness 

45.00 

Adapter plate 

8.00 

Bus System (dry) 

1314.00 

Structure 

298 00 

AACS 

102.00 

Cabling 

15.00 

Propulsion 

62.00 

Thermal control 

30.00 

Power 

404.00 

Command and data 

152.00 

Mechanisms 

148.00 

Communications 

103.00 

Miscellaneous 

771.62 

Launch adapter/misc 

100.00 

Propellant 

90.00 

Project reserve 

100.00 

Contingency 

481 .62 

TOTAL (kg) 

3179.70 

MARGIN (kg) 

120.30 
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The importance of the shared aperture approach is 
shown in Figure 8, which illustrates the incidence 
angle performance of comparable shared and separate 
aperture antennas. Both antennas are assumed to have 
a fixed 18.6 x 2.6 m area available for the aperture. In 
the shared aperture approach, the L-band aperture uses 
the full 2.6 m width. In the separate aperture approach, 
the 2.6 m width is apportioned between the L. C, and 
X-band apertures so that the L-band aperture has a 
width of only 1.9 m. As a result, the maximum 
incidence angle decreases from 55° to 47° for single 
polarization and from 42° to 33° for dual/quad 
polarization. 

Conclusion 

The LASAR/MASAR feasibility study shows that SAR 
antennas large enough to satisfy the needs of the 
science user community can indeed be built using non- 
exotic SIR-C heritage technology and launched with 
the low-cost launch vehicles. 
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Figure 8 Comparison of ambiguity and beam width 
limited swath versus incidence angle for a) shared 
aperture, and b) non-shared aperture. 


Table 4 Science requirements and compliance of the MASAR spacecraft. 


REQUIREMENT 

VALUE 

DESIGN VALUE 

Frequency 

L. (C.X) 

L. (C. X) | 

Polarization 

quad, (quad, dual) 

quad, (quad, dual) 

Resolution 



- Local mode 

20-30 m 

30 m (3 looks) 

- Regional mode 

50-100 m 

50-100 m 

- Global mode 

250 m 

250 m 

Incidence angle range 



- Single polarization 

15-50 deg 

TBD-55 deg 

* Quad polarization 

15-40 deg 

TBD-42 deg 

Thermal noise equiv o° @ 45° 

-36 dB. (-27 dB, -22 dB) 

-37 dB. (-28 dB, -23 dB) 

Max single swath width 



- Single polarization 

30-50 km 

within ambiguity limits 

- Quad polarization 

30-50 km 

within ambiguity limits 

Max SCANSAR swath 

600 km 

>600 km 

Calibration 



* Relative amplitude 

+/- 1 .0 dB 

+/- 1 .0 dB 

- Absolute amplitude 

♦/- I.OdB 

+/- 1 .0 dB 

- Relatvie phase 

TBD - 

TBD 

- Geometric 

+/- 0.5 pixel 

+/- 0.5 pixel 

Ambiguity level 

<-20 dB 

<-23 dB 

Sidelobe level 



- Elevation 

<-14 dB 

<-15 dB 

- Azimuth 

<-12 dB 

<-13 dB 

Polarization isolation 

<-25 dB 

<-30 dB 

Orbit exact repeat interval 

<5 days 

2 days 

Orbit near repeat interval 

16 days 

16 days 
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ABSTRACT 

The Synthetic Aperture Radar (SAR) has viewed, in the last years, a huge increment of interest from 
many preset and potential users. 

The good spatial resolution associated to the all weather capability lead to consider SAR not only a 
scientific instrument but a tool for verifying and controlling the human daily relationships with the Earth 
Environment. New missions have been identified for SAR as spatial resolution became lower than three 
meters : disasters , pollution , ships traffic ,vo!canic eruptions, earthquake effect are only few of 
possible objects which can be effectively detected, controlled and monitored by SAR mounted on 
satellites. This paper deals with high resolution radar design constraints and dimensioning . 


INTRODUCTION 


The recent years have been characterized by a increasing of interest on remote sensing by SARs both 
for civil and not civil applications. The technology development can now assure new features to SAR 
sensors and spacecraft missions can be better tailored to the specific user requirement. 

The two key parameters for new missions are: revisit time and resolution. 

The first parameter deals with the mission definition, coverage, orbits selection and the number of 
satellites which can be allocated in orbit at the same time. The high resolution permits the analysis of 
manmade or natural objects with the capability of their detection and in same case their identification. 
Furthermore high resolution leads to better radiometric resolution for the conventional applications 
where pixel dimensions requirement are less stringent. A resolution better than three meters is generally 
requested for these Survelliance Missions but the target is to reach the one meter resolution within next 
decade. Civil applications could be identified in the controlling and monitoring of environment areas and 
phenomena like: 

- Pollution 

- Fires 

- Volcanic activity 

- Ships traffic control 

- Land traffic control 

- Environmental disasters 

- Crisis monitoring 

- Urbanization control 

- coast erosion 

The high resolution is generally attractive in countries like Europe where we find high population density 
and limited land extension so the number of objects to be discriminated is very large. 
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Alenia Spazio has been involved , in the last few years, in a number of studies for high resolution SAR 
sensor and their potential applications. In this note only few considerations will be provided identifying 
the most significative aspects of high resolution SAR design. 

REQUIREMENT OVERVIEW 

The main requirement of an high resolution SAR sensor are here below identified: 


parameter 

unit 

high resolution mode 

survey mode 

range resolution 

meters 

1 to 3 

5 to 100 

azimuth resolution 

meters 

1 to 3 

5 to 100 5 

swath extension 

Kmeters 

5 to 10 

50 to 400 

sensitivity 

dB 

-20 

-20 

number of looks 

- 

1 

1 to 16 i 

azimuth ambiguities ratio 

dB 

-22 

-22 

range ambiguity ratio 

dB 

-22 

-22 

radiometric accuracy 

dB 

< 2 

< 1 

access area 

degree 

30 tO 60 + 

20 to 60 + 

polarization 

- 

HH or VV 

HH or VV/ HH and 
VV 

altitude 

Kmeters 

550 to 800 

550 to 800 


The high resolution is obtained in the range direction by wide bandwidth (BW) chirp transmission. The 
figure 1 . shows the bandwidth needs as function of angle of incidence for a defined resolution. The 
minimum BW for one meter resolution is 300 MHz if 30 Deg incidence angle is requested. Now if we 
look to the International Regulation about frequency allocation for RADAR applications ( S: 200 MHz, 

C: 100 MHz, X: 300 MHz ,Ku: 600 MHz ), then we discover that only X and Ku bands can be used 
for the scope. 

The azimuth resolution can be obtained by a conventional stripmap mode with two meters antenna length 
or the most promising spotlight mode which consists of increasing the conventional synthetic aperture by 
tracking with the antenna beam the ground swath. This implies azimuth steering capability in the 
antenna. Figure 2. gives the frequency dependence of resolution in case only the quadratic term of the 
SAR echo is compensated into the processing. The present compression algorithms should be improved 
for high resolution processing looking in particular to the depth of focus and to the combination of the 
discrete synthetic sub apertures produced by the antenna steering. 

The Swath dimensions for high resolution SAR is limited mainly by the data rate but should be pointed 
out that in case of spotlight mode the antenna dimensions are heavily affected because the azimuth 
beam width should fit with swath dimension in the along track direction. Then the antenna height is 
dictated by ambiguities (minimum antenna area) and link budget considerations. The result of 
dimensioning appears in an antenna with 15/20 square meters area, narrow in the along track 
direction. Therefore in same cases, when it is requested to combine high resolution with a survey- 
medium resolution / medium swath mode, it becomes difficult to find swath coverage solution without 
using the complex Scansar technique . Alternatively antenna beam forming by TR modules 
weighting or switching-off a certain number of TR modules, so reducing the antenna height , are 
techniques to be attempted. The drawbacks of spotlight mode are the need of an active antenna with 
azimuth steering and the gaps on ground between an images and the next one ( fig 4 .). The gap for an X 
band sensor is in the order of one/ two times depending on satellite altitude and incidence angle, the 
Swath dimension in the azimuth direction. The access area ( fig 3. ) is important in determining the 
revisit time of the sensor and the incidence angle at which a "spotswath" could be viewed. The access 
area lower limit comes from resolution/ BW requirement while the upper limit from link budget/ antenna 
dimension constraints. 
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The multi polarization doesn’t seem essential to surveillance mission SAR adding complexity to the 
system without dramatically improving detection capability. The spacecraft altitude is defined from 
mission requirement and impacts specifically high resolution in antenna dimension and data rate. 

The data rate of an instrument like we are outlining in this notes has a data rate in the order of one 
Gigabit/sec. This high amount of data can be transmitted to ground directly or by a Data Relay Satellites 
which are under study for the next future. 

CONCLUSIONS 

ALS has studied the high resolution implementation on SAR sensor . The result shows that the design is 
feasible in X band where resolution can be achieved and technology is available. The selected solution is 
based on spotlight mode: active antenna is requested with azimuth steering capability (1/2 Deg) for 
azimuth resolution and elevation steering capability ( 20/ 30 Deg ) for imaging within the access area. 


Range resolution 


Max achievable resolution (azimuth) 

Frequency dependence 
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Microwave Remote Sensing Investigations 
Planned for BOREAS 

N94-15914 

Kyle C. McDonald 
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4800 Oak Grove Drive 
Pasadena, CA 91109, U.S.A. 




Summary 

The Boreal Ecosystem - Atmosphere Study (BOREAS) is a multidisciplinary field and 
remote sensing study that will be implemented jointly by the United States and 
Canada. The goal of BOREAS is to obtain an improved understanding of the 
interactions between the boreal forest biome and the atmosphere in order to clarify 
their roles in global change. Specific objectives are to improve the understanding of 
the processes that govern the exchanges of water, energy, heat, carbon, and trace 
gases between boreal ecosystems and the atmosphere, and to develop and validate 
remote sensing algorithms for transferring the understanding of these processes 
from local to regional scales. 

Two principal field sites, both within Canada, have been selected for this study. The 
northern site is located near Thompson, Manitoba, and the southern site encompasses 
Prince Albert National Park in Saskatchewan. The growing season in the northern 
site tends to be limited by growing-degree days while the southern site is limited by 
soil moisture and fire frequency. Most of the field work will occur at these two sites 
during 1993 and 1994 as part of six field campaigns. The first of these campaigns is 
scheduled for August 1993 and will involve instrument installation and an 
operational shakedown. Three large scale Intensive Field Campaigns (IFCs) are 
scheduled for 1994, along with two smaller scale Focused Field Campaigns (FFCs). The 
first 1994 campaign will be an FFC designed to capture the biome under completely 
frozen conditions during the winter. The second FFC and the first IFC are scheduled to 
capture the spring thaw period. Another IFC will take place in the summer during a 
period of maximum water stress. Finally, the third FFC will be scheduled to capture 
the collapse into senescence during the fall. 

The BOREAS science team has been organized into six disciplines, one of which 
centers on remote sensing science investigations. The five microwave investigations 
have been included in this group. These activities consist of four radar studies and 
one passive microwave investigation. Titles of these studies are listed below together 
with the principal investigators and a brief statement of their primary objectives: 

(1) Microwave scatterometer observations of boreal forest species. 

M. C. Dobson, The University of Michigan 

Objective: Obtain c- and x-band scatterometer measurements of 2-3 tree 

species within the two test sites. Backscatter from each species 
will be observed for 4-5 days continuously during each of the 
three IFCs. These measurements will then be used to verify 
linkages between biophysical controls and SAR observations. 
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(2) Distribution and structure of biomass in boreal forest ecosystems. 

K. J. Ranson, NASA Goddard Space Flight Center 

R. H. Lang, George Washington University 

Objective: Determine the amount and spatial distribution of above-ground 

biomass. The planned activity includes development of 
algorithms for determining the amount standing biomass and its 
apportionment from SAR backscatter, development of model 
inversion techniques, and determination of the optimum radar 
parameters for use in estimating canopy biomass. 

(3) Estimation of hydrological parameters in boreal forest using SAR. 

S. Saatchi, Jet Propulsion Laboratory, California Institute of Technology 

Objective: Estimate hydrologic parameters of the forest and underlying 

ground surface. The water content of the vegetation and ground 
surface and surface cover will be examined with SAR. Hydrologic 
parameters of soil, snow, moss and permafrost will be studied. 

(4) Monitoring environmental and phenologic state and duration of state with 

SAR as input to improved C02 flux models. 

J. B. Way, Jet Propulsion Laboratory, California Institute of Technology 

K. C. McDonald, Jet Propulsion Laboratory, California Institute of Technology 

Objective: Classify boreal forest by functional group and determine the 

environmental and phenologic state and water status within each 
functional group. In so doing, SAR will be used to determine the 
photosynthetically active period of these functional groups. This 
information will be provided as new input for estimating the CO2 
flux of the biome. 

(5) The diurnal and annual radiobrightness of boreal forest. 

A. W. England, The University of Michigan 

Objective: Monitor and model diurnal and seasonal radiobrightness 

signature of boreal forest. This information will be used to update 
existing cold region radiobrightness models to account for 
seasonal and diurnal signatures of the boreal forest. 

These studies have been designed to investigate both the static and temporally 
varying microwave signatures of the boreal forest biome. Multifrequency 
polarimetric SAR imagery will be collected by aircraft-borne systems during the 
field campaigns. This will allow detailed analysis of the static and temporally varying 
aspects of backscatter. During the intervals between these campaigns, imagery from 
the ERS-1 satellite will provide more of a continuous monitoring capability. 

This article was prepared at the Jet Propulsion Laboratory, California Institute of 
Technology, under contract to the National Aeronautics and Space Administration. 
The research described in this article is co-sponsored by the National Aeronautics 
and Space Administration and Energy, Mines and Resources, Canada. 
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Institute of Remote Sensing Applications, Academia Sinica 

Wu Guoxiang 

National Remote Sensing Center of China 
Wang Zhen-song 

Institute of Electronics, Academia Sinica 


The Shuttle Imaging Radar-A and -B (SIR-A and SIR-B) carried on the Space 
Shuttle Columbia in November 1981 and the Challenger in October 1984 
acquired images of test sites of North China. The Russian ALMAZ SAR also 
acquired imagery of part of this test site in September, 1992. In November of 

1990, the airborne SAR developed by the Chinese Academy of Sciences 

(CAS/SAR) covered this area for the purpose of Chinese spacebome radar 
development. By studying and analyzing these SAR data we have achieved 
positive results in geoscientific applications. 

Penetration. It is discovered that SIR-A has penetrated through the thin 
dry sand sheet at the Altengaobao area in Inner Mongolia. A triangular 

geologic body, which has a strong radar return, is clearly shown on the SIR-A 
image (Fig. 2a). The central part of the triangular geological body is totally 
covered by dry sand which is shown on the Landsat MSS image (Fig. 2b). 
Digging a shallow hole at the central part, the bedrock is exposed at about 1 
meter of sand in the field verification. The geological map of this region 
shows that the upper part of the triangular body is Precambrian metamorphic 
rock, the lower part is comprised of Mesozoic granites, and the central part, 

covered by dry sand, is gneisses. The theoretical study indicates that the 
penetrating ability of radar should at least meet three conditions: 1) fine 

grain size; 2) less thickness; 3) extremely dry. The surface conditions of the 
test site are content with the former two conditions. The meteorological 
records of the test site show that the annual mean rainfall is less 150 mm, 

some places even less 10 mm, dryness is 4 - 12 (the dryness of adjacent desert is 
7 - 12), and evapotranspiration is 31.5 times as rainfall, which means that the 
test site is content with the third condition. 

Rock Type and Fault Recognition. The distribution and range of 

Caledonian plagioclase-granite at Wang Yang area in Gansu Province has been 
detected by using SIR-B data. Due to adverse natural circumstances and poor 
traffic, the test site is geologically under-studied, with only a few plagioclase 
granite bodies being displayed on the geological map. The large 

Bayinnorgong strike slip fault, 80 km long, was recognized on CAS/SAR and 
SIR-A images at the Bayinnorgon area of inner Mongolia, and the field work 
proved its existence. Several bunches of faults were interpreted from the SIR- 
A image of the Keping area in Xinjiang Province. The longest one is 50 km. 
Many secondary faults are well developed along the two sides of the major 
fault which cuts through all of the strata in the test site. The strata on the two 
sides are strongly distorted, deformed, and displaced, and the displacement is 
over 5 km. 
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Mineral Exploration. SAR is useful for detecting ore-controlling 
structures. A gold-bearing structure belt was directly discovered by the 

interpretation of SIR-A images in Honguyulin area. The belt is about 650 km 
long and composed of Sinian sandstone, limestone and quartz-diorite. 
Sampling along the belt, the chemical analysis results show that the gold tenor 
of most samples reaches 1 ppm and the highest sample is 7.94 ppm. The 
discovery of the gold-bearing structure belt has two leading factors: 1) The 

linear strong return anomaly. On one hand the strong return comes from the 

rough surface of the belt. On the other hand, it is caused by the high dielectric 

constant of altered rocks. The dielectric constant of the country rock is less 
than 7, but the mineralized rock is over 20. 2) The relevant analysis result. 

The gold-bearing structure belt is a fractured structure belt, distributing 
along and contacting with the granite body, which is an ideal location for 
metal mineral deposit. 

Multi-parameter SAR Data Comparison. Jilantai area of Inner Mongolia 
is our spacebome radar science test site. SIR-A, SIR-B, ALMAZ SAR and 
CAS/S AR data of the test site have been acquired so far (Fig. 1). Their 
parameters are as follows: Operation Bands: L, L, S and X; Spatial Resolution: 

40mx40m, 32mx28m, 10-15m and lOmxlOm; Look Angle: 47°, 26°, 45°, and 19°; 

Flight Direction: NW 80°, NE 20°, NW 20°, and NE 90°. All four SAR systems 

operate in HH polarization. With the diverse imaging parameters, they have 
different imaging effects. 

The high-resolution imagery is more suitable for rock type discrimination. 
The detecting ability is better for desert areas when the radar illuminating 
direction is perpendicular to the slipface or sand dunes. The small depression 
angle will bring a lot of shadow and reduce the information. 

We would like to acknowledge Dr. Charles Elachi and Dr. Tom Farr 

of JPL. Also NPO and ALMAZ Corporation, who kindly provided 

SIR-A, SIR-B and ALMAZ SAR data to our study. 
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Figure 1. Jilantai test site of North China 


SIR-A Image (upper left) SIR-B Image (upper right) 

ALMAZ SAR Image (lower left) CAS/SAR Image (lower right) 



Figure 2. Altengaobao area of Inner Mongolia, North China 

SIR-A image showing penetration Landsat MSS image (right), 
phenomena (left). 
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PRIMARY STUDIES OF CHINESE SPACEBORNE SAR 


Wang Zhensong, Wu Guoxiang, Guo Huadong, Wei Zhongquan, Zhu Minhui 



ABSTRACT 

This paper deals with the primary studies on spacebome synthetic aperture 
radar (SAR) in China. The SAR will be launched aboard a Chinese satellite and 
operated at L-band with HH polarization. The purpose of the mission in 
consideration is dedicated to resources and environment uses, especially to 
natural disaster monitoring. The ground resolution is designed as 25 m x 25 m 
for detailed mode and 100 m x 100 m for wide scan-SAR mode. The off-nadir 
angle can be varied from 20 to 40 degrees. The key system concepts are 
introduced here. 


I. INTRODUCTION 

It is well known that SAR can map the Earth surface under poor optical 
conditions. Spacebome SAR data can be widely used for natural resources 
exploration and environment monitoring. SAR is also suggested as the most 
important payload for a natural disaster reduction satellite system, which has 
been considered by many countries and international organizations. To meet 
the requirement for all-weather Earth observation of remote sensing users 
both from China and around the world, China is considering the development 
of a Chinese spacebome SAR system. The primary studies on related 
techniques are conducted in the framework of the National High-Tech 
Research and Development Plan of China, which is to develop techniques to be 
industrialized and practically applied by the end of this century. The primary 
studies on spacebome SAR show that the SAR operated at L-band in a sun- 
synchronized circular orbit can partly fulfill the tasks suggested by Chinese 
users for flood disaster monitoring and other resources and environment 
related fields. 

The SAR will use horizontally polarized L-band (wavelength 23.5 cm) for both 
transmitting and receiving of the SAR signal. Two operating modes, detailed 
mode and wide scan-SAR mode, can be chosen by ground instructions 
according to different working tasks. The satellite will also have a limited 
orbit change capability to meet the need of flood monitoring in the eastern 
part of China, where flood disasters every summer are common. This 
capability means that the satellite can be changed to a lower orbit which is 
frozen to over East China. 

In past years, China has used airborne SAR to monitor flooding areas during 
the flood season. An airborne SAR system developed in China played an 
important role in this kind of monitoring. This system, working in X-band 

with four polarizations, has a resolution of 10 m x 10 m, a swath width of 35 km, 
and an acting range of 150 km. This experience contributed greatly to the 
studies on spacebome SAR in China. 
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The primary results presented here are based on the mission study, on 
domestic experience with airborne SAR systems, and on SAR data processing. 


n. MISSION 

The Chinese spacebome SAR mission will no doubt meet many requirements 
for SAR data in resources and environment uses. Applications of importance 
to China also includes flood and drought monitoring. The SAR system will be 
onboard a Chinese satellite with a circular sun-synchronized orbit. The orbit 
parameters are shown as Table 1. 


Table 1. Parameters of Chine se SAR Satellite 


Orbit height 

Orbit inclination 

Orbit period 

Repetitive period 

Number of revolutions per day 

Total number of revolutions 

Local time of descending node 

Orbit regression 

Space between neighbor orbits 


608 km 

97.8° 

96.83 min. 

31 day 
14+27/31 
461 
10:30 

24.2° westward 
81.7 km (equator) 
75.5 km (20°) 


The planned satellite will be a 3-axis stabilized one. The orbit is chosen 
according to the consideration of Earth coverage and better S/N ratio in the 
SAR operation under the power allocated to the SAR system. 


The SAR system is designed according to user requirements that demand a 
variety of incidence angles from 20 to 40 degrees. A microstrip array antenna 
measuring 12.5 m x 2.5 m is considered with electronic beam steering 
capability. Selected radar characteristics are listed as Table 2. 


Table 2. Chinese Spacebome S AR Characteristics 


Frequency 
Wavelength 
Polarization 
Off-nadir angle 
Spatial resolution 

Swath width 
Pulse bandwidth 
Pulse length 
Pulse repeat frequency 
Noise equivalent 
Quantization bit number 
Down link data rate 


1.275 GHz 
23.5 cm 
HH 

20-40° 

25 m / 100 m 
100 km / 400 km 
20 MHz 
33 micro sec. 
Selectable 
< -20 dB 
3 or 5 
100 Mb/s 
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The SAR system will transmit a linear frequency modulated signal, which is 
generated by the use of a surface acoustic device. A solid state high-power 

amplifier will be used for transmitting high-peak power probing signals. The 

SAR electronics for use onboard the satellite are being developed. The 

primary studies on the spaceborne SAR show the possibility that the 

requirements of spaceborne SAR to the space platform can be realized. 

Techniques for spaceborne SAR image processing are developed in the 

primary studies and SAR images have been reconstructed. Studies on the 

applications of spaceborne SAR data are also being carried out in geology, 
hydrology, agriculture, forestry, oceanology, and surveying and mapping. We 
believe that a mature user group will be the key to the successful realization of 
the future Chinese spaceborne SAR project. 

For peaceful uses of space remote sensing technology, we are looking forward 
to having international collaboration with other countries in the field of 

spaceborne imaging radar development. 


Dr. Wang Zhensong and Prof. Zhu Minhui, Institute of Electronics, Academia 
Sinica; Mr. Wu Guoxiang, National Remote Sensing Center, State Science and 
Technology Commission of China; Prof. Guo Huadong, Institute of Remote 
Sensing Applications, Academia Sinica; Prof. Wai Shongquan, Shanghai 
Institute of Satellite Engineering. 
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A GROUP DISCUSSION AT THE 
3rd SPACEBORNE IMAGING RADAR SYMPOSIUM 
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SAR GROUND DATA SYSTEMS 
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- COMMON FORMATS: DATA STRUCTURES, COMMON READERS 

- COMMON PRESENTATION: PRODUCT DEFINITION 

- DATA PRODUCT QUALITY STANDARDS: CALIBRATION, ALGORITHMS 
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SECURITY OF DATA SYSTEMS 

- CONTROL OF UNAUTHORIZED DISTRIBUTION 

- CONTROL OF UNAUTHORIZED ACCESS 
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GLOBAL ENERGY AND WATER CYCLE EXPERIMENT 
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The Complex Role of Water 
in Climate Processes 









BASELINE SURFACE 
RADIATION NETWORK 
(BSRN) 
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CONTINENTAL-SCALE 
INT'L PROJECT 
(GCIP) 




435 


GCIP STRATEGY: 
DEVELOPING COUPLED 
LAND-SURFACE/ATMOSPHERE MODELS 
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SEA 


GC1P. GCSS. AMIP AND 
SCALE-INTERACTIVE MOIST PROCESSES 
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GCIP FIELD CAMPAIGNS 


o 

Ujjcc 

oP 


CO-J 

m< 


co Si 


sen 

^LU 

IU5 


< ^ 

Scog 

C0-3< 

O^tii 

CO -Jo 

LLI<Q 

*E“ 

->.£0. 

oS| 

(oSs 


Q< 

UIQ 

o 3 

gfc 

2* 

co5 

in Q 

QZ 


so 

X? 
3 CO 

LL§ 

HI CO 

tuj 

<h- 

So 

-js 

<LU 

>CC 


439 




'■'■-Ola- • 

/ZSS73 

N 9 4 " 1 5 93. 9^ b 

AIRSAR South American Deployment 

M. Kobrick 

Jet Propulsion Laboratory 
California Institute of Technology 

Pasadena, California , ; 



PRECEDING PAGE BLANK NOT FILMED 


441 



AIRSAR South American Deployment 

The United States National Aeronautics and Space Administration (NASA) and the 
Brazilian Commission for Space Activities (COBAE) are undertaking a joint experiment 
involving NASA’s DC-8 research aircraft and the Airborne Synthetic Aperture Radar (AIRSAR) 
system during late May and June 1993. The research areas motivating these activities are: 

(1) fundamental research in the role of soils, vegetation, and hydrology in the global carbon 
cycle, and (2) in cooperation with South American scientists, airborne remote sensing research 
for the upcoming NASA Spacebome Imaging Radar (SIR)-C/X-SAR flights on the Space 
Shuttle. 

The following pages show a flight schedule and plans for the deployment that were 
developed at the Ames Research Center and the Jet Propulsion Laboratory (JPL) by using the 
flight requests received and approved by NASA Headquarters in response to a Dear Colleague 
letter dated May 29, 1992. The approved Principal Investigators (Pis), investigation titles, and 
approximate site locations are also listed. 

The plan development included an attempt to accommodate all data requirements in the 
approved requests within the resources available, with the primary limitations being the length of 
the deployment (37 days) and the total number of DC-8 flight hours available within that time 
frame. The only major restriction imposed by those limits was the deletion of the sites south of 
about 24°S latitude. Most of the other approved data requests are accommodated by this plan, 
subject to some caveats listed below. 

Included in the following pages are maps of the site locations and schematic indications of 
flight routes and dates, plots showing swath locations derived from the flight requests and 
generated by flight planning software and, most importantly, a calendar showing which sites will 
be imaged each day. Although the scheduling may not follow any obvious logical pattern, it is 
the result of incorporating a large number of disparate and sometimes conflicting requirements 
such as the desire to acquire multitemporal data at sites, coordination with satellite overflights, 
and so on. It should also be understood that the schedule is not set in concrete and will continue 
to evolve in response to changing PI experiment plans, changes imposed by the acquisition of 
international flight clearances, and any anomalies experienced during the deployment itself. 
Every attempt will be made to coordinate any and all such changes with the Pis involved, and the 
establishment and maintenance of good lines of communication both during the planning and 
operations phases will be of utmost importance. 

The schedule has been developed based upon the following assumptions: 

1. Both the DC-8 and the AIRSAR systems will function normally and any hardware problems 
will be minimal. A small amount of margin has been built into each day’s flight plan, but 
since the deployment cannot be extended beyond the planned number of days, any major 
hardware anomalies will necessarily result in the loss of data. Although historically the 
reliability of both systems has been exceptionally high, and in fact the AIRSAR has not 
failed to acquire planned data due to hardware or system problems in its five years of 
operation, it must be remembered that the deployment is done on a best-effort basis and 
aircraft or radar problems could result in schedule changes or, in the worst case, the deletion 
of some sites. 


442 


2. Clearance will be obtained through the NASA International Affairs Office to overfly and 
acquire radar imaging data from each of the countries involved. At this writing, this has not 
yet been completed and historically it has been difficult to obtain flight clearance for some 
countries. 

3. There will be no significant problems with Air Traffic Control (ATC) or with flying into 
military or other restricted flight zones. The DC-8 flies at the same altitudes as commercial 
air traffic and each day’s flight plan must be cleared and approved by the proper ATC 
authorities. In the past, this has occasionally resulted in changes to flight plans or take off 
times to minimize interference with other air traffic, particularly in foreign countries. Work 
so far indicates that no sites violate any military or restricted flight zones, but a final 
determination will have to wait for the development of more detailed flight plans by the 
DC-8 navigators. 

4. Unless otherwise requested by the PI, all AIRS AR data will be acquired in the standard 
three-frequency quad-polarization mode. Where TOPSAR mode data have been requested, 
C-band polarimetric data will not be acquired but, instead, a CW image will be produced 
along with a digital topographic map registered to the C-band image and L- and P-band 
polarimetric data sets. Note that although JPL has developed operational TOPSAR processor 
software, hardware for the systematic processing of TOPSAR data is not yet available. Thus 
processing of all TOPSAR mode data will be done by arrangement as special products 
outside the normal AIRSAR processing system. 

In addition to AIRSAR, two other interesting and important experiments will also be 
installed on the DC-8 and will acquire data during this deployment: the lightning detection and 
lidar experiments, which will utilize upward-looking optical sensors and will prefer to have data 
acquired during dark conditions, i.e., during early evening or night and with no moon. Since 
AIRSAR flights generally have morning takeoffs so that during the daylight the video and film 
cameras may operate and the Pis collecting surface truth data at ground sites may work, some 
adjustments will be made. 

The plan for the deployment will be to have afternoon takeoffs whenever weather forecasts 
indicate that a thunderstorm imaging opportunity will occur. This will allow the acquisition of: 
SAR data during daylight, and lightning and lidar data during the latter part of the flight in the 
early evening. These opportunities are most likely to occur during mid June 1993, and in the 
central and western portions of Brazil, Bolivia, and Peru. Unfortunately this means that exact 
take-off times cannot yet be determined for each flight, but every attempt will be made to 
coordinate take-off and overflight times with the Pis in the field during the deployment itself. 
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PI 

Affiliation 


Title 


Site 


approx Lon, Lat 


John R. Baker 

British National Space Centre 

AIRSAR Estimation of Tropical Forest Biomass for 
Carbon Cycle Studies 

Manaus-Reserva Du. 
Sena Maduriera 
Rio Tapajos 

-59.57 

-68.67 

-55.00 

-2.57 

-9.25 

-3.50 

Paolo Canuti 

Univ. of Florence 

Contribution of S AR for Estimating Soil Erosion 

Bermejo Basin 
Homahuaca 

-63.00 

-65.50 

-2350 

-23.75 

Luciano Dutra 
INPE 

Rain Forest Type Discrimination, Rood Plain 
Mapping, Calibration Studies, and SAR Model 

Rio Tapajos 
Cachoeira Porteira 

-54.80 

-57.05 

-3.05 

-0.65 


Assessment Using AIRSAR Data: Tapajos Region 


Tony Freeman 
Jet Propulsion Lab 

JERS-l Calibration Underflight (Amazon) 

Manaus 
Rio Tapajos 

-60.10 

-55.20 

-3.20 

-3.70 

Ron Greeley 

The Relationship Between Radar Backscatter and 

Cerro Quisharo, Alti 

-68.75 

-18.75 

Arizona State Univ. 

Aerodynamic Roughness 




Dirk H. Hoekman 

Remote Sensing of Tropical Rain Forests 

Araracuara, Columbia 

-72.07 

0.65 

Wageningen Agricultural Univ. 

Mabura Hill, N, Guyana 

-57.90 

6.70 



Mabura Hill, S, Guyana 

-58.64 

5.08 



San Jose del Guaviare 

-72.61 

2.47 


Bryan Isacks 

S1R-C Analysis of Topography and Climate in the 

Cordillera Blanca, Peru -77.53 

. 

Cornell University 

Central Andes 

Cordillera Real, Bolivia -68.29 

-1 



Quelccaya Ice Cap -70.42 

-1 



Potosi -65.63 

-1 



Hermann Kux Forestry and Land Use in Western Amazonia, Acre Sena Maduriera -68.45 -9.00 

INPE State Brazil 


John Kwiatkowski Applications of Polarimetric Radar in Coherent ^ - t f 

Michigan Technological Univ. Imaging of Earth and Remote Sensing of Precipitation 


Thuy Le Toan 

Centre d’Etude Spatiale des 
Rayonnements 

Study of Tropical Rain Forest by SAR Data 
in French Guiana 

Kaw 

St Elie Paracou 
St Laurent 

-52.10 

-53.03 

-53.97 

4.72 

528 

5.48 

- 

John M. Melack 

Determining the Extent of Inundation on Subtropical 

Anavilhanas 

-60.78 

-2.64 


Univ. of California, Santa 

Roodplains Beneath Vegetation of Varying Types 

Cabal i ana 

-60.90 

-3.40 


Barbara 


Pantanal 

-57.42 

-19.67 

- 

Peter J. Mouginis-Mark 

Radar-Derived Topography of Volcanoes in the 

Galapagos 

-91.00 

-1.00 

ir 

Univ. of Hawaii 

Western Galapagos Islands 
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PI 

Affiliation 


Title 


Site 


approx Lon, Lat 


Waldir Paradella 

Geobotanical Investigation and Lithological 

Curaca Valley 

-39.83 

-9.21 

INPE 

Discrimination by AIRS AR Data in the NE of Brazil 




Jack Paris 

Global Biodiversity: Assessment of Habitat 

Manaus area 

-59.75 

-2.50 

Calstate Fresno 

Change and Species Extinctions with 





Multinanometer SAR 




Kevin 0. Pope 

SIR-C Tropical Wetlands Project 

Merida 

-90.25 

20.75 

GeoEcoArc Research 


Rio Bravo 

-89.03 

17.52 


, 

Rio Hondo 

-88.55 

18.27 

Jeffrey E. Richey 

Land Use, Forest Type and Biomass Assessment 

Manaus-Fazenda Di. 

-59.97 

-2.43 

Univ. of Washington 

Using Polarimetric SAR 

Manaus-Reserva Du. 

-60.10 

-2.93 

Davis D. Sentman 

Optical Imaging of Cloud-tolonosphere Electrical 

Thundercloud sites 



University of Alaska Fairbanks 

Discharges in Brazil 




David Sheres 

SAR Measurements of Rain, Waves, Fronts and 

Gulf of Mexico 

-89.00 

27.00 

Stennis Space Center 

Eddies 




David Skole 

Assessing Secondary Growth and Carbon 

Nova Vida 

-62.79 

-10.19 

University of New Hampshire 

Accumulation in Disturbed Tropical Forests Using 

Paradise 

-63.31 

-9.58 


Airborne Synthetic Aperture Radar 




Joao Vianei Soares 

Surface Hydrology of Semiarid Floodplains 

Bebedouro 

-40.28 

-9.08 

INPE 

(Varzeas) in Northeast Brazil (SIR-C/X-SAR) 




Vincente Paulo Soares 

Relationship between Radar Image Backscatter and 

Minas Gerais 

-42.63 

-19,80 

Federal Univ. of Vicosa 

Eucalyptus Stand Characteristics 




John W. Terborgh 

Vegetation Mapping of a Headwater Catchment in 

Manu National Park 

-12.05 

-7035 

Duke University 

the Amazon Basin: Manu National Park, Peru 




Jakob J. vanZyl 

Use of SAR Imagery in Monitoring the Hydrology of 

Bebedouro 

-40.28 

-9.08 

Jet Propulsion Lab 

Semi-Arid Floodplain Areas in Brazil 




Howard A. Zebker 

Topsar Topographic Mapping Research and Support 

Manaus area 

-59.75 

-2.50 

Jet Propulsion Lab 

for the South American DC-8 Deployment 

Galapagos Islands 

-91.00 

-1.00 


445 


The Schedule 

The following calendars for May and June 1993 show each flight day for the DC-8, indicated by an 
aircraft icon. The notation Data Flight indicates flights where the DC-8 takes off from and returns 
to the same base, and Transit to... indicates takeoff from one base and landing at another, generally 
with data being acquired en route. 

May 25: The DC-8 will depart from its home base at Moffett Field, California and will acquire oceano- 
graphic data over the Gulf of Mexico before landing at Houston International Airport. 

May 27: After one day in Houston the aircraft will transit to San Jose, Costa Rica with data acquired en 
route over the Gulf of Mexico and ecology sites near Merida, Mexico and Rio Hondo and Rio 
Bravo, Belize. 

May 29: After one day in San Jose the next flight will be a transit to Guayaquil, Ecuador via the Galapagos 
Islands. The objective will be to map all of Isla Femandina and Isla Isabela in TOPSAR mode, as 
indicated on the attached map, but as currendy planned this flight is near the duration limit for a 
one day DC-8 flight and a slighdy smaller data set may be acquired. 

May 31: During this flight from Guayaquil the sites at Aracuara and San Jose del Guaviare in Colombia 
will be imaged, as well as the set of lines called the“volcano transect.” Not yet planned in detail, 
these lines will image a number of Andean volcanoes along a track intended to cover as many 
individual sites as possible. 

During the transit to Santa Cruz, Bolivia TOPSAR mode data will be acquired at Cordillera 
Blanca and the Quelccaya Ice Cap in the Peruvian Andes, and the site at Manu National Park will 
also be imaged. 

Data sites for this flight from Santa Cruz include Cordillera Real and Potosi, which will be 
imaged in TOPSAR mode, and Cerro Quisharo near the Chilean border. 

Flying south into Argentina, AIRSAR will collect data at the hydrology sites at Bermejo Basin 
and Homahuaca, then the rest of the flight will be dedicated to Dave Sentman’s lightning imaging 
experiment Since the lightning experiment is weather dependent, the scheduling for this flight 
may be modified in real time to match the opportunities that arise. 

This transit flight to Recife, Brazil will cover the Pantanal site in Brazil and the SIR-C/X-SAR 
“super site” at Bebedouro. 

June 11: This flight will be second of four to overfly Bebedouro, the intent being to acquire data under 
different ground conditions, primarily soil moisture. For that reason the date of this and the 
following flight may be adjusted slightly according to weather conditions. Also covered on this 
day will be the site at Minas Gerais. 

June 12: During this flight data will be acquired at the Curaca Valley site as well as Bebedouro. 

June 15: This transit flight into Manaus provides one last opportunity to image Bebedouro, and in transit 
the site at Rio Tapajos will be covered in both SAR and TOPSAR modes. 

June 17: This will be the first of two flights over the Manaus area, with the plan being to acquire data along 
tracks matching the ground swaths of both the JERS-1 and Shuttle Imaging Radar experiments. 
The dates of this and the following flight may be adjusted to accommodate the lightning imaging 
experiment. 


June 2: 

June 4: 
June 6: 

June 8: 
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June 19: This flight will be divided between the acquisition of data along certain specific tracks in the 
Manaus area for Howard Zebker’s forest investigation, and at targets of opportunity for the 
lightning imaging experiment Due to the weather-dependent nature of the lightning experiment, 
the flight may occur on a different day or be swapped with the preceding flight 


June 21: During this flight data will be acquired at the SIR-C/X-SAR site at Sena Madureira as well as at 
Paradise and Nova Vida. Finally, data lines will be flown at the Rio Tapajos site to match an 
overflight of the JERS-1 satellite which will be imaging the region on that day. 


June 24: The transit to San Juan, Puerto Rico will provide coverage at the Cachoreira Porte ira site northeast 
of Manaus, Kaw, Paracou and St. Laurent along the coast of French Guiana, and Mabura Hill in 
Guyana. 


June 26: This last flight of the deployment will be a transit back to Moffett and no AIRSAR data will be 
acquired. 
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Date 

Latitude 

Longitude 

Transil Data Total 


d m 

d m 

to 



25-May 

37 

24.9 

-122 

2.8 

Moffett 




27 

55.0 

-93 

35.0 

Gulf of Mexico 

3.4 

2.0 



29 

27.0 

-95 

12.6 

Houston 
Block time 

0.3 

0.8 


6.45 

27-May 

29 

27.0 

-95 

12.6 

Houston 




27 

55.0 

-93 

35.0 

Gulf of Mexico 

0.3 

1.0 



20 

45.0 

-90 

15.0 

Merida 

1.0 

0.8 



18 

16.2 

-88 

33.0 

Rio Hondo 

0.4 

0.5 



17 

31.2 

-89 

1.8 

Rio Bravo 

0.1 

0.5 



9 

33.6 

-84 

3.0 

San Jose 

1.2 








Block time 

0.8 


6.57 



31 -May 

-2 

6.0 

-79 30.0 

Guayaquil 





0 

58.8 

-77 54.9 

Mayasquer volcano 0.5 

0.3 



1 

1.8 

-77 48.3 

Cumbal volcano 

0.0 

0.3 



1 

9.0 

-77 30.9 

Azurfal volcano 

0.0 

0.3 



1 

22.2 

-77 6.9 

Galeras volcano 

0.1 

0.3 



1 

56.7 

-76 39.3 

Dona Juana vol. 

0.1 

0.3 



2 

41.1 

-76 10.8 

Huila volcano 

0.1 

0.3 



3 

49.5 

-75 40.5 

Tolima volcano 

0.2 

0.3 



4 

45.9 

-75 22.2 

Ruiz volcano 

0.1 

0.3 



5 

6.9 

-75 34.6 

Herveo volcano 

0.1 

0.3 



2 

28.2 

-72 36.6 

San Jose del Gu. 

0.5 

0.5 



0 

39.0 

-72 4.2 

Araracuara 

0.3 

0.5 



-2 

6.0 

-79 30.0 

Guyaquil 

1.1 







Block time 

0.8 


6.99 

2-Jun 

-2 

6.0 

-79 30.0 

Guayaquil 





-9 

15.0 

-77 32.0 

Cordillera Blanca 

1.0 

1.3 



-12 

2.0 

-71 12.0 

Manu Ntl Park 

0.9 

1.2 



-13 

57.5 

-70 47.5 

Quelccaya Ice Cap 0.3 

1.2 



-17 

28.8 

-63 6.0 

Santa Cruz 

1.1 







Block time 

0.8 


7.59 

4-Jun 

-17 

28.8 

-63 6.0 

Santa Cruz 





-16 

11.5 

-68 16.5 
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